Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,0390TG45.0% 93.8 / 16.6 61N108T (AAC→ACC) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (8/25);  new base G (27/0);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 4.48e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  NC_000913/2677892‑2678187
                                                                                                                                                   |                                                                                                                                                      
aGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGtttt                                                                                                                                                     <  1:77532/151‑1 (MQ=255)
 gctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTc                                                                                                                                                    <  2:391075/151‑1 (MQ=255)
 gctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCCGGCAGAGGgttt                                                                                                                                                     >  1:385489/1‑150 (MQ=255)
    gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCg                                                                                                                                                  <  2:915968/150‑1 (MQ=255)
    gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGc                                                                                                                                                 <  2:76333/151‑1 (MQ=255)
        ttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACGTACCCCGGCAGAGGGTTTCCGCCAc                                                                                                                                              >  1:654514/1‑150 (MQ=255)
        ttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGAACGAGGCAGAGGGTTTTCGcc                                                                                                                                                >  2:665643/1‑148 (MQ=255)
          gcgcCGGGCTTAGCGTTGAAACAAGTGTTTCAGGACGAACGTCAACGCCCAGCGCCTGTAGGCAACGTTGGGCAGCCTGCGCCATTTGCAGGTAATCAACCATGCCGTTGCGGCGGGGCCATTGCCCGAGGCAGAGTTTTTCCGCCACTG‑t                                                                                                                                          <  1:18438/151‑1 (MQ=255)
           cgcCGGGCTTGGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACt                                                                                                                                             <  1:985421/148‑1 (MQ=255)
                 gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCCACCGAGGCAGAGGGTTTTCGCCCCTG‑TCCGccc                                                                                                                                    >  1:778318/1‑150 (MQ=255)
                 gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTCCCCCACTG‑TCAGccc                                                                                                                                    >  2:367013/1‑150 (MQ=255)
                   ttAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCAGGCAGGGGGTTTTCGCCCCCG‑TCAGcccct                                                                                                                                  >  1:120656/1‑148 (MQ=255)
                      gCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCCGAGGGTTTCCGCCCCTG‑GCACCCCTTcctc                                                                                                                              >  1:769698/1‑149 (MQ=255)
                           gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACACACCCAGGCAGGGGATTTCTGCCACTG‑TCAGCCCCTCCA‑CCCGAc                                                                                                                        >  1:544982/1‑151 (MQ=255)
                             aaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACtt                                                                                                                      <  2:759122/151‑1 (MQ=255)
                             aaCAAGTTGTTCAGGACTAACGCCAACGCCCAGCGCCTGTAAGCAACGTTGGGAATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACAGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACtt                                                                                                                      <  1:664182/151‑1 (MQ=255)
                               cAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCCAGGCAGAGGGTTTTCGCCACCG‑TCAGCCCCTccc‑ccccacct                                                                                                                      >  1:671445/1‑143 (MQ=255)
                               cAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTaa                                                                                                                    <  2:261263/151‑1 (MQ=255)
                                aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAAc                                                                                                                   <  2:1042286/151‑1 (MQ=255)
                                    tgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGCGGCGGAGGGTTTTCCCCACCG‑TCAACCCTTCCC‑CCAGACTTAAATcc                                                                                                                >  1:534878/1‑150 (MQ=255)
                                       tCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCCGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGa                                                                                                             <  2:376452/150‑1 (MQ=255)
                                       tCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACACAACCAGGGAGAGGGTTTTCGCCCCCC‑TCAACCCCTCCA‑CCACACTTTACttcc                                                                                                               >  1:413027/1‑145 (MQ=255)
                                        cAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGAt                                                                                                            <  2:415741/150‑1 (MQ=255)
                                        cAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATa                                                                                                           <  1:165382/151‑1 (MQ=255)
                                         aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATaa                                                                                                          <  2:460108/151‑1 (MQ=255)
                                             cTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAAcc                                                                                                       <  2:175226/150‑1 (MQ=255)
                                                aCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACcgcgc                                                                                                   <  2:1011965/151‑1 (MQ=255)
                                                  gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCAGAGGGTTTTCGCCACTG‑TCAGCCCTTCCC‑CCCGACTTAAATCCTGATAAAAcgc                                                                                                     >  2:130770/1‑147 (MQ=255)
                                                   tCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGAGGGTTTTCGCCACTG‑TCAGCCCCTCCC‑CCAGAATTAACTCCTGCAAAACCccccgc                                                                                                 >  2:935116/1‑150 (MQ=255)
                                                    cAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGAACGAGGGAGAGGGTTTTCGCCCCTG‑TCAGCCCCTCCC‑CCCGACCTAACTCCTTATaaaaccc                                                                                                     >  2:679903/1‑141 (MQ=255)
                                                    cAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCAc                                                                                               <  2:594258/151‑1 (MQ=255)
                                                     aaCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCCGAGGGTTTTCGCCCCTG‑TCAGCCCCTTCC‑CCAGAATTAAATCCCGGTAAACCCCGCGCAcc                                                                                              >  1:1075462/1‑151 (MQ=255)
                                                     aaCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCAGGCAGAGGGTTTTCCCCCCTG‑TCACCCCCTCCC‑CCCGACTTTAATCCTGATAAAACCCGCGcccc                                                                                              >  2:171731/1‑151 (MQ=255)
                                                          ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAg                                                                                         <  2:596140/151‑1 (MQ=255)
                                                            aGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCAGGGAGAGGGTTTTCGCCCCCG‑TCAGCCCCTCCC‑CCCGACTTAACTCCCGATAAAACCCGCGCACCCCCCGtt                                                                                       >  2:937648/1‑151 (MQ=255)
                                                            aGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGtt                                                                                       <  2:915062/151‑1 (MQ=255)
                                                             gcgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTc                                                                                      <  2:804027/151‑1 (MQ=255)
                                                               gcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAg                                                                                    <  1:286698/151‑1 (MQ=255)
                                                               gcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTCCGCCCCTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCACGCGCAccccccattc                                                                                      >  2:278169/1‑144 (MQ=255)
                                                                 cTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGGCCGAGGCAGAGGGTTTTCCCCCCAG‑TCAGCCCCTCCACCCAAACTTAA‑TCCCGATAAAACGCCCCCAccccccggtccacg                                                                                  >  1:793648/1‑142 (MQ=255)
                                                                  tGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTTACCGCGGGAGAGGGTTTTCGCCCCTG‑TTAGCCCCTCCC‑CCAGACCTAAATCCTGATAAAACGCGCGCAccccccccttagccg                                                                                 >  1:900402/1‑141 (MQ=255)
                                                                         acctttGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACgg                                                                            <  1:605836/145‑1 (MQ=255)
                                                                           cGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCTGCCCCGGGAGAGGGTTTTCCCCACTG‑TCAGCCCCTTCT‑CCCCACCTTACTCCCTATAAAACCCCCGCCCCCCCCGTTCCGCGGCACcgcgg                                                                         >  2:289761/1‑149 (MQ=255)
                                                                                 gCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCAGAGGGTTTTCCCCCCCG‑TCAGCCCCTCCA‑CCCGACTTTACTCCTTATAAAACCCGCGCCCCCCCAGTTCAGCGGCACCGCGAGTCAGc                                                                   >  2:1081400/1‑150 (MQ=255)
                                                                                   aTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTa                                                                <  2:90617/151‑1 (MQ=255)
                                                                                       tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGAACGAAGCAGAGGGTTTCCGCCACTG‑TCAGCCCCTCCA‑CCCAACCTAACTCCTGATAAACCGCGCGCCCCCCCCATTCAGCGGCACGGCGAGTCAGCGTAgact                                                            >  2:506888/1‑148 (MQ=255)
                                                                                         cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGAGGCAGAGGGTTTTCCCCCCCG‑TCAGCCCCTCCC‑CCCGACTTAAATCCTGGTAAACCCCGCGCCCCCCCCCCTCAGCGGCCCCGCGGGGCAGCGTAGcct                                                            >  1:736654/1‑147 (MQ=255)
                                                                                          gcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGGCAGCGTAGCtt                                                            >  2:75733/1‑148 (MQ=255)
                                                                                                 gCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTCAACCGGGGGGAGGGTTTCCCCCACTGTTCAGCCC‑TCCT‑CCAGCCTTAACTCCTTTAAAACCCCGCGCCCCCCCACTTCCGCGGCCCGGCGAGTCAGCGTTGCGTCGTCAcc                                                     >  1:740139/1‑148 (MQ=255)
                                                                                                   aGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGGAGAGGGTTTTCGCCCCTG‑TCAGCCCCTTCA‑CCCGACCTTACTCCTTATAAACCGCGCGCACCCCCCGGTCAGCGGCACGGCGAGGCAGCCGAGCCTCGTCACCCTcc                                                 >  1:410254/1‑150 (MQ=255)
                                                                                                    ggTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCa                                                <  2:892598/150‑1 (MQ=255)
                                                                                                          caatcaTGCCGTTGCGGCGGGGCCACTTACCCAGGCAGGGGTTTTTCTCCCCTG‑TCAGCCCTTCCC‑CCAGACCTAACTCCTTATAAAACGCCCGCACCCCCCGTTCAGCGGCACGGCGAGACAGCGTAGCTTCTTCACCCTCCAGTcgc                                           >  2:359142/1‑149 (MQ=255)
                                                                                                                      tGCGGCGGGGGCCCTGACCGCGGCAGAGGGTTTTCCCCCCTG‑TCAGCCCCTCCA‑CCCCCCCTAACTCCCGATAAAACCCGCGCCCCCCCCATTCAGCGGCCCCGCGCGGCAGCGTAGCTTCGTCACCTTCCACTCGCGCCTCCCCAATcc                              >  1:1075428/1‑150 (MQ=255)
                                                                                                                        cggcggGGCCACTGACCCAGGCAGAGGGTTTCCGCCACTG‑TCAGCCCTTCCA‑CCCGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGCCACCCTCCAGTCGCGTCTCGCCAATcc                              >  1:749918/1‑148 (MQ=255)
                                                                                                                         ggcggGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGat                          <  1:94386/151‑1 (MQ=255)
                                                                                                                           cggGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat                        <  1:531225/151‑1 (MQ=255)
                                                                                                                               gCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCAcc                    >  2:139107/1‑151 (MQ=255)
                                                                                                                                 cACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGc                  >  2:1050481/1‑151 (MQ=255)
                                                                                                                                   cTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCt                 >  2:1023676/1‑150 (MQ=255)
                                                                                                                                      aCCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCCCGTCTCGCCAATCCAGATATCACCGCTATcc             >  2:709142/1‑151 (MQ=255)
                                                                                                                                             cAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGtt      >  2:162111/1‑151 (MQ=255)
                                                                                                                                                agGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGc   >  1:87062/1‑151 (MQ=255)
                                                                                                                                                 ggTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  >  2:804617/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
AGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTG‑TCAGCCCTTCCA‑CCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  NC_000913/2677892‑2678187

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑