Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,188 A→T 18.8% intergenic (+22/+31) yicH → / ← yicI putative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,1880AT18.8% 259.4 / 50.3 117intergenic (+22/+31)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (42/52);  new base T (8/14);  total (51/66)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-01

Only 100 of 118 total aligned reads displayed.

AGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGA  >  NC_000913/3832090‑3832287
                                                                                                  |                                                                                                   
cttccgatctGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGAc                                                                                                    <  1:1696784/91‑1 (MQ=255)
agcagcAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGAc                                                                                                    <  2:2361323/100‑1 (MQ=255)
agcagcAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGAc                                                                                                    <  2:2342595/100‑1 (MQ=255)
agcagcAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGAc                                                                                                    <  1:1969040/100‑1 (MQ=255)
agcagcAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGAc                                                                                                    <  2:1408979/100‑1 (MQ=255)
  cagcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTc                                                                                                  >  2:1222140/1‑100 (MQ=255)
   agcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGACGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCa                                                                                                 <  2:717521/100‑1 (MQ=255)
   agcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCa                                                                                                 <  2:599949/100‑1 (MQ=255)
       aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                             >  1:242143/1‑100 (MQ=255)
          gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGAcagatcggaag                                                                                          >  2:1696772/1‑90 (MQ=255)
          gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAc                                                                                          >  2:1129230/1‑100 (MQ=255)
           aCTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAcc                                                                                         <  2:1317100/100‑1 (MQ=255)
            cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                        <  1:2016904/100‑1 (MQ=255)
            cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                        >  1:1699843/1‑100 (MQ=255)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                      <  1:2559615/100‑1 (MQ=39)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                      <  2:1355357/100‑1 (MQ=39)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                      <  1:1427392/100‑1 (MQ=39)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                      <  1:1318947/100‑1 (MQ=39)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                      <  1:2062920/100‑1 (MQ=255)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                      >  2:2446020/1‑100 (MQ=255)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                      <  1:956317/100‑1 (MQ=255)
              cAAACCACGAATGCTGGCAGCGTTTCCAGCAGTTAAGTTACATCGACGGAGCCTGTGCGCTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                      <  2:1516406/100‑1 (MQ=37)
              atctccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                      <  1:1920418/96‑1 (MQ=38)
               aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                     <  1:1711116/100‑1 (MQ=39)
               aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                     <  1:1645911/100‑1 (MQ=39)
               aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                     <  2:1463047/100‑1 (MQ=39)
               aaaCCATGAAGGCTGTCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                      <  2:1750266/99‑1 (MQ=255)
                aaCCATGAATGCTGGCATCGTTGCCAGCGGTGTAGTTACGTCGACTGAGCCGGTGCGGTAATCATCAAGCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                      <  2:251220/98‑1 (MQ=37)
                  ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCggagat                                                                                  >  2:1920404/1‑96 (MQ=38)
                  ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcgag                                                                                  <  1:121881/100‑4 (MQ=38)
                  ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGCACGCCTCAGACGTACCCGGCAt                                                                                   >  2:2200672/1‑99 (MQ=255)
                    ttccgatcTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTcc                                                                                <  2:12263/93‑1 (MQ=255)
                     tGAATGGTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                      <  2:1078132/93‑1 (MQ=38)
                       aaTGCTGGCATCGTTTCCAGCTGTGAAGTTACATCGACGGAGCTGGTGCGGTAATCATCAATCTCGCCGTGTACGACTCAGACGTACCCGGCATTCCATc                                                                             <  2:1463848/100‑1 (MQ=255)
                       aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                             <  2:1703964/100‑1 (MQ=255)
                       aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                             <  2:25346/100‑1 (MQ=255)
                       aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                             <  1:1731849/100‑1 (MQ=255)
                          gCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAAt                                                                          <  2:1961581/100‑1 (MQ=255)
                            tGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAg                                                                        <  2:1841315/100‑1 (MQ=255)
                             ggCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCagatcggaa                                                                       >  1:12264/1‑92 (MQ=255)
                             ggCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGa                                                                       >  1:1517831/1‑100 (MQ=255)
                             ggCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGa                                                                       >  2:366220/1‑100 (MQ=255)
                              gCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAt                                                                      <  1:628072/100‑1 (MQ=255)
                                 tCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTACACATCAATCCCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAtt                                                                     <  2:1848272/98‑1 (MQ=255)
                                   gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                 >  1:32289/1‑100 (MQ=255)
                                   gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                 >  2:2320770/1‑100 (MQ=255)
                                   gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAAATTACa                                                                  >  2:504688/1‑99 (MQ=255)
                                    tttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                                <  1:1057009/100‑1 (MQ=255)
                                     ttCGAGCGGTGACGTTACATCGAAGGAGCCGGTGCGGTAAACATCAAGGTCACCGCGATGGACTAAGACCTACCCGTCATTAAATCAACAGATTACAACg                                                               >  2:246152/1‑100 (MQ=255)
                                     ttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACg                                                               <  1:1236387/100‑1 (MQ=255)
                                      tCCAGCGGTGAAGTGACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAAGAATAGATTACAACg                                                               >  2:2103348/1‑99 (MQ=255)
                                            ggTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGt                                                        >  2:1468592/1‑100 (MQ=255)
                                                aaGTTACATCGACGGAGCTGGTGCGGTAAACATCAATATCGCGGGGTACGAGTGAGACGTAACCGGCATTCCATCAATAGATTACAACGTAATTGTCAgc                                                    >  2:2434971/1‑100 (MQ=255)
                                                aaGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgc                                                    <  1:2284275/100‑1 (MQ=255)
                                                 aGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcg                                                   >  1:2221859/1‑100 (MQ=255)
                                                     aCATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAtt                                               >  1:1861569/1‑100 (MQ=255)
                                                      cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTc                                              >  2:864413/1‑100 (MQ=255)
                                                       aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                             >  1:2291731/1‑100 (MQ=255)
                                                       aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                             >  2:1936382/1‑100 (MQ=255)
                                                         cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                           >  2:695731/1‑100 (MQ=255)
                                                          gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCCATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                          <  2:873509/100‑1 (MQ=255)
                                                           aCGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                         >  2:49858/1‑100 (MQ=255)
                                                             ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                       <  1:2477244/100‑1 (MQ=255)
                                                             ggAGCCGGTGCGGTAAACATCAATCTCGCCGGCTACGACTCAGCCGTACCCGGCATTCCATCAATCGATTACAACGTAATTGCCAGCGCATTCCCTTGAg                                       <  1:117073/100‑1 (MQ=255)
                                                              gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                      >  2:2103694/1‑100 (MQ=255)
                                                              gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                      <  1:2375002/100‑1 (MQ=255)
                                                               aGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGc                                     >  1:22261/1‑100 (MQ=255)
                                                               aGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGc                                     >  1:760954/1‑100 (MQ=255)
                                                                  cGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCttn                                  <  2:920606/100‑2 (MQ=255)
                                                                  cGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTc                                  <  1:2446034/100‑1 (MQ=255)
                                                                   ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                 <  1:1400437/100‑1 (MQ=255)
                                                                    gTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcac                                <  2:1911570/100‑1 (MQ=255)
                                                                     tGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacc                               <  1:534456/100‑1 (MQ=255)
                                                                     tGCGGTAAACATAAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACTTAATTGTCCGCGCATTCCCTTGAGGCTTca                                 >  2:289200/1‑98 (MQ=255)
                                                                        ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                            >  2:278217/1‑100 (MQ=255)
                                                                        ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                            >  2:276693/1‑100 (MQ=255)
                                                                        ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                            <  1:2508291/100‑1 (MQ=255)
                                                                         gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCAccacca                           >  2:2140030/1‑100 (MQ=255)
                                                                             aCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGccc                       >  1:2221606/1‑100 (MQ=255)
                                                                              cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                      <  1:136573/100‑1 (MQ=255)
                                                                              cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                      <  1:444805/100‑1 (MQ=255)
                                                                               aTCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCt                     >  2:1943158/1‑100 (MQ=255)
                                                                                  aaTCTCGCTGGGTACGACTCAGACGTACCCGGCATGCCATCGATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                   <  2:1172041/99‑1 (MQ=255)
                                                                                  aaTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCt                  >  1:1207344/1‑100 (MQ=255)
                                                                                    tctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCa                <  1:278221/100‑1 (MQ=255)
                                                                                     ctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCAc               <  2:198130/100‑1 (MQ=255)
                                                                                       tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGCTTACAACGTAATTGTGCGCGCATTCACTTGAGGCTTCACCACCAGCCCCTGCTCACtt             >  2:565064/2‑100 (MQ=37)
                                                                                       tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt             >  1:2395536/2‑100 (MQ=39)
                                                                                       cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtn             <  2:896727/100‑2 (MQ=255)
                                                                                        gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCCCttt            >  2:894022/1‑100 (MQ=38)
                                                                                        gCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt              <  1:2435878/98‑1 (MQ=255)
                                                                                         ccGGGTACGTCTGAGTCGTACCCGGCATTCCCTCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc           >  2:106590/1‑100 (MQ=39)
                                                                                         ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc           >  1:1697378/1‑100 (MQ=39)
                                                                                          cGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTAAGACGGAATTGTCAGCGAATTCGCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa          >  2:1011158/1‑100 (MQ=255)
                                                                                            ggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGc        >  1:1955219/1‑100 (MQ=255)
                                                                                             gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc       >  1:2545760/1‑100 (MQ=255)
                                                                                                cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc    <  1:475725/100‑1 (MQ=255)
                                                                                                cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc    >  2:1624338/1‑100 (MQ=255)
                                                                                                 gACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCg   <  1:818171/100‑1 (MQ=255)
                                                                                                  aCTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGa  >  2:1611083/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
AGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGA  >  NC_000913/3832090‑3832287

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: