Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,456,828 G→A 7.1% intergenic (+36/‑145) gapB → / → PP_2150 glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,456,8280GA 7.1% 289.0 / 8.6 99intergenic (+36/‑145)gapB/PP_2150glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (49/42);  new base A (2/5);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 104 total aligned reads displayed.

GACAACCGCGTTGTTCTGTACGTCTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGC  >  NC_002947/2456679‑2456976
                                                                                                                                                     |                                                                                                                                                     
gACAACCGCGTTGTTCTGTACGTCTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCg                                                                                                                                                       <  2:84978/150‑1 (MQ=255)
 aCAACCGCGTTGTTCTGTACGTCTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCgg                                                                                                                                                      <  1:1945875/150‑1 (MQ=255)
 aCAACCGCGTTGTTCTGTACGTCTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCgg                                                                                                                                                      <  2:1459989/150‑1 (MQ=255)
   aaCCGCGTTGTTCTGTACGTCTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGt                                                                                                                                                    >  2:1494691/1‑150 (MQ=255)
          ttgtgCTGTACGTCGGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTTGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGtt                                                                                                                                            <  2:77944/150‑1 (MQ=255)
           tgttCTGTACGTCTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGttt                                                                                                                                           <  1:918832/150‑1 (MQ=255)
           tgttCTGTACGTCTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCCGAAGGTTCCCGttt                                                                                                                                           <  2:1008234/150‑1 (MQ=255)
            gttCTGTACGTCTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGtttt                                                                                                                                          >  1:468300/1‑150 (MQ=255)
            gttCTGTACGTCTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGtttt                                                                                                                                          <  1:1388692/150‑1 (MQ=255)
            gttCTGTACGTCTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCCGAAGGTTCCCGtttt                                                                                                                                          <  2:1543766/150‑1 (MQ=255)
             ttCTGTACGTCTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGttttt                                                                                                                                         >  1:2315406/1‑150 (MQ=255)
                       cTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTg                                                                                                                               <  1:2243499/150‑1 (MQ=255)
                          gTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAAt                                                                                                                            >  2:1266068/1‑150 (MQ=255)
                             cGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAAtgtg                                                                                                                         >  2:694189/1‑150 (MQ=255)
                              gACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAAtgtgt                                                                                                                        <  2:1385825/150‑1 (MQ=255)
                               aCAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAAtgtgtg                                                                                                                       <  1:1494691/150‑1 (MQ=255)
                                 aaCGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGtgtt                                                                                                                     >  2:1472280/1‑150 (MQ=255)
                                  aCGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCg                                                                                                                                                       <  2:1357486/116‑1 (MQ=255)
                                  aCGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCg                                                                                                                                                       >  1:1357486/1‑116 (MQ=255)
                                    gAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTc                                                                                                                  >  1:646951/1‑150 (MQ=255)
                                     aaTTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCa                                                                                                                 <  2:2310554/150‑1 (MQ=255)
                                       ttCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGt                                                                                                               <  1:1472280/150‑1 (MQ=255)
                                               taGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGccc                                                                                                       <  2:1637196/150‑1 (MQ=255)
                                                     ccAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCACAATTCCCACCAAGATTGGCATGTCCTGCTCCAGAACTTTTTCCACCTCGTCGTTTCGTTTTGTTTTTATTCTTAATGTGTGTTGTCAGTGCCGGCCCTTTcgc                                                                                                 <  2:1433422/150‑1 (MQ=255)
                                                       aGGTGGTTCGGGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGG                                                                                               <  2:2315406‑M1/150‑1 (MQ=255)
                                                        ggtggtTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGG                                                                                              <  2:359694‑M1/150‑1 (MQ=255)
                                                         gtggtTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGT                                                                                             >  2:1754157‑M1/1‑150 (MQ=255)
                                                         gtggtTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGT                                                                                             >  2:1424221‑M1/1‑150 (MQ=255)
                                                         gtggtTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGT                                                                                             >  2:422068‑M1/1‑150 (MQ=255)
                                                          tggtTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTA                                                                                            <  2:47621‑M1/150‑1 (MQ=255)
                                                              tCGTGTGATGGAAGACATGGCTGGTGTGAACTCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACC                                                                                        <  1:1679741‑M1/150‑1 (MQ=255)
                                                                    gATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTT                                                                                                    >  2:2284798‑M1/1‑132 (MQ=255)
                                                                    gATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTT                                                                                                    <  1:2284798‑M1/132‑1 (MQ=255)
                                                                     aTGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGTCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGtgttgt                                                                                                                   <  1:2244601/116‑1 (MQ=255)
                                                                     aTGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGTCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGtgttgt                                                                                                                   >  2:2244601/1‑116 (MQ=255)
                                                                         aaGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAG                                                                             <  1:694189‑M1/150‑1 (MQ=255)
                                                                              aTGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGtt                                                                                                                                                   <  1:261563/76‑1 (MQ=255)
                                                                              aTGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGtt                                                                                                                                                   >  2:261563/1‑76 (MQ=255)
                                                                                  cTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAAtgtgt                                                                                                                        <  2:52649/98‑1 (MQ=255)
                                                                                  cTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAAtgtgt                                                                                                                        >  1:52649/1‑98 (MQ=255)
                                                                                      tgtgAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGcc                                                                >  1:702373/1‑150 (MQ=255)
                                                                                       gtgAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCATTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCc                                                               <  2:1658840‑M2/150‑2 (MQ=255)
                                                                                         gAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTg                                                                  >  2:1821485/1‑145 (MQ=255)
                                                                                         gAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTg                                                                  <  1:1821485/145‑1 (MQ=255)
                                                                                            cccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAACACGTGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgta                                                          >  1:1424221‑M2/1‑144 (MQ=255)
                                                                                             ccGCCAGCGTTTCCACGTTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtag                                                         >  1:2055375‑M2/1‑143 (MQ=255)
                                                                                             ccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtag                                                         <  1:2233110‑M2/150‑8 (MQ=255)
                                                                                                 cAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGggcagcgataatcctgtaggagc                                                     >  2:1149621‑M1/1‑127 (MQ=255)
                                                                                                 cAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGggcagcgataatcctgtaggagc                                                     <  1:422068‑M1/150‑24 (MQ=255)
                                                                                                  aGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGG                                                                            >  2:1976894‑M1/1‑126 (MQ=255)
                                                                                                  aGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGG                                                                            <  1:1976894‑M1/126‑1 (MQ=255)
                                                                                                   gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcgg                                                   >  1:1040565‑M2/1‑137 (MQ=255)
                                                                                                    cGTTTCCACGCTGATTCGCTTGTAAGGTAATGAAAAAACGGGAACGTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGGCCTTTCGTGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcggg                                                  >  1:620096‑M2/1‑136 (MQ=255)
                                                                                                    cGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcggg                                                  >  2:702373‑M2/1‑136 (MQ=255)
                                                                                                    cGTTTCCACGCGGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcggg                                                  >  1:1266068‑M2/1‑136 (MQ=255)
                                                                                                     gTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcgggt                                                 >  2:1826697‑M2/1‑135 (MQ=255)
                                                                                                         ccACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCATTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcgggtttac                                             >  1:1658840‑M2/1‑131 (MQ=255)
                                                                                                           aCGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcgggtttaccc                                           >  2:2233110‑M2/1‑129 (MQ=255)
                                                                                                            cGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTC                                                                                   >  1:1689297‑M1/1‑109 (MQ=255)
                                                                                                            cGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTC                                                                                   <  2:1689297‑M1/109‑1 (MQ=255)
                                                                                                             gCTGATTCGCTTGTAAGGCAATGAAAAACCGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTA                                                                                <  1:1898538‑M1/112‑1 (MQ=255)
                                                                                                             gCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTA                                                                                >  2:1898538‑M1/1‑112 (MQ=255)
                                                                                                             gCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATcctgtaggagcgggtttacccgc                                         <  1:1530592‑M2/150‑24 (MQ=255)
                                                                                                             gCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcgggtttacccgc                                         >  2:468300‑M2/1‑127 (MQ=255)
                                                                                                                 aTTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATcctgtaggagcgggtttacccgcgaaa                                     >  2:1530592‑M2/1‑123 (MQ=255)
                                                                                                                 aTTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGG                                                                                               >  2:261251‑M1/1‑92 (MQ=255)
                                                                                                                 aTTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGG                                                                                               <  1:261251‑M1/92‑1 (MQ=255)
                                                                                                                  ttCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcg                                                    <  2:2279092‑M2/134‑13 (MQ=255)
                                                                                                                  ttCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcg                                                    >  1:2279092‑M2/1‑122 (MQ=255)
                                                                                                                     gCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCCGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                 <  2:1802985/150‑1 (MQ=255)
                                                                                                                      cTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGcc                                <  1:1754157/150‑1 (MQ=255)
                                                                                                                       ttGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATcctgtaggagcgggtttacccgcgaaagggccg                               >  1:1149621‑M2/1‑117 (MQ=255)
                                                                                                                         gTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcgggtttac                                             >  2:1436017‑M2/1‑115 (MQ=255)
                                                                                                                         gTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcgggtttac                                             <  1:1436017‑M2/134‑20 (MQ=255)
                                                                                                                          tagggCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCa                            >  1:2329793/4‑150 (MQ=255)
                                                                                                                            aGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGTCCTTTCGCGGGTAAACCATCTCCTACAAGTTTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGGATAGGGCAGCAcc                          >  2:1088156/1‑150 (MQ=255)
                                                                                                                             ggCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCg                         <  2:800949/150‑1 (MQ=255)
                                                                                                                               cAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGc                       >  1:1899586/1‑150 (MQ=255)
                                                                                                                                   aaaaaaacGGGAACCCGAAGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcgggtttacccgcgaa                                      <  1:1158626‑M2/130‑27 (MQ=255)
                                                                                                                                   aaaaaaacGGGAACCCGAAGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCcctgtaggagcgggtttacccgcgaa                                      >  2:1158626‑M2/2‑105 (MQ=255)
                                                                                                                                    aaaaaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCg                  <  2:477390/150‑1 (MQ=255)
                                                                                                                                    aaaaaaCGGGAACCCGAAGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGTTATCGCTGCCcctgtaggagcggg                                                  <  2:249123‑M2/118‑15 (MQ=255)
                                                                                                                                    aaaaaaCGGGAACCCGAAGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGTTATCGCTGCCcctgtaggagcggg                                                  >  1:249123‑M2/1‑104 (MQ=255)
                                                                                                                                      aaaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGggcagcgataatcctgtaggagc                                                     <  1:323481‑M1/113‑24 (MQ=255)
                                                                                                                                      aaaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGggcagcgataatcctgtaggagc                                                     >  2:323481‑M1/1‑90 (MQ=255)
                                                                                                                                      aaaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAG                                                                             <  1:478704‑M1/89‑1 (MQ=255)
                                                                                                                                      aaaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAG                                                                             >  2:478704‑M1/1‑89 (MQ=255)
                                                                                                                                       gaaCTGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCgaga               >  1:451003/2‑150 (MQ=255)
                                                                                                                                       aaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGTAGCGGGTTTACCCGCGAAGAGGACAGAAACGGCAACCgaga               >  1:295467/1‑150 (MQ=255)
                                                                                                                                       aaaCGGGAACCCGAAGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCg                  <  1:2033561/147‑1 (MQ=255)
                                                                                                                                         aCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGaaa             >  1:976491/1‑150 (MQ=255)
                                                                                                                                          cGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAAc            >  2:726743/1‑150 (MQ=255)
                                                                                                                                          cGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACA                                                                              <  2:1806652‑M1/84‑1 (MQ=255)
                                                                                                                                          cGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGACCTTTCGCGGGTAAACCCGCTCCTACA                                                                              >  1:1806652‑M1/1‑84 (MQ=255)
                                                                                                                                            ggAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAAccc          <  2:1899586/150‑1 (MQ=255)
                                                                                                                                             gAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAaga                                    <  1:1945333/123‑1 (MQ=255)
                                                                                                                                                 cTTCGGGTTCCCGTTTTTTTTTTTTCTTAATTTGTTTTGTCAGTGCCGGCCCTTTCCCGGGTAAACCCGCCCCTACAGGATTATCCCTCCCCTTGTAGGAGGGGGTTTACCGCCGAAGGGGCCACCCCCGGCAACAGGAAAAcccccctt     >  2:1418541/1‑145 (MQ=255)
                                                                                                                                                 cTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCgg                        >  2:591656/1‑131 (MQ=255)
                                                                                                                                                 cTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCgg                        <  1:591656/131‑1 (MQ=255)
                                                                                                                                                    cGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTCGTCAGTGCCGGCCCTTTCGCGGGCAAGCCCGCTCCTACAGTAGTAACGCTGCCCCTGTAGGAGCGGGGTTAACCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTgc  <  1:1616322/150‑1 (MQ=38)
                                                                                                                                                     |                                                                                                                                                     
GACAACCGCGTTGTTCTGTACGTCTGGTACGACAACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGC  >  NC_002947/2456679‑2456976

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: