Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,659,886 T→C 5.6% intergenic (+43/‑215) aceA → / → PP_4118 isocitrate lyase/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,659,8860TC 5.6% 1047.5 / 30.2 323intergenic (+43/‑215)aceA/PP_4118isocitrate lyase/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (165/140);  new base C (12/6);  total (177/146)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 332 total aligned reads displayed.

AGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCG  >  NC_002947/4659737‑4660035
                                                                                                                                                     |                                                                                                                                                     
agGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCt                                                                                                                                                       >  1:6063900/1‑150 (MQ=255)
agGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCt                                                                                                                                                       <  2:6729598/150‑1 (MQ=255)
agGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCt                                                                                                                                                       >  2:4207420/1‑150 (MQ=255)
 ggAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTc                                                                                                                                                      <  1:3769438/150‑1 (MQ=255)
  gAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTcc                                                                                                                                                     <  2:1416107/150‑1 (MQ=255)
      tGGGCACTGGCAACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGc                                                                                                                                                 <  2:8225475/150‑1 (MQ=255)
       gggCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGcc                                                                                                                                                >  1:702218/1‑150 (MQ=255)
         gCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCtt                                                                                                                                              >  1:561335/1‑150 (MQ=255)
          cACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCttt                                                                                                                                             <  2:3033336/150‑1 (MQ=39)
           aCTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTg                                                                                                                                            >  2:4229302/1‑150 (MQ=255)
               gCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGtttt                                                                                                                                        >  1:6561650/1‑150 (MQ=255)
                cTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGttttt                                                                                                                                       >  1:7479116/1‑150 (MQ=255)
                 tACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTa                                                                                                                                      >  2:8753294/1‑150 (MQ=255)
                 tACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTa                                                                                                                                      <  2:4949903/150‑1 (MQ=255)
                  aCTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAc                                                                                                                                     <  1:3630084/150‑1 (MQ=255)
                  aCTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTGCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAc                                                                                                                                     <  2:7599835/150‑1 (MQ=255)
                    ttCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAccc                                                                                                                                   <  1:3306891/150‑1 (MQ=255)
                    ttCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTAccc                                                                                                                                   <  1:7589260/150‑1 (MQ=255)
                     tCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCg                                                                                                                                  <  1:669165/150‑1 (MQ=255)
                           acATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAg                                                                                                                            <  1:1774454/150‑1 (MQ=255)
                           acATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCATTGTTTTTACCCGAGAAAg                                                                                                                            <  2:8008337/150‑1 (MQ=255)
                             aTGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCCAGAAAggg                                                                                                                          >  1:8254974/1‑150 (MQ=255)
                                    ccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccgggcc                                                                                                                   >  1:6940774/1‑150 (MQ=255)
                                    ccGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAgggccgggcc                                                                                                                   >  1:854950/1‑150 (MQ=255)
                                       tGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCctgc                                                                                                                <  1:5094449/150‑1 (MQ=255)
                                          tCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTgg                                                                                                             >  1:5295201/1‑150 (MQ=255)
                                          tCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTgg                                                                                                             >  1:7190788/1‑150 (MQ=255)
                                            cAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGcc                                                                                                           >  2:2677461/1‑150 (MQ=255)
                                             aGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCa                                                                                                          >  2:6656324/1‑150 (MQ=255)
                                                  ggcgCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAAtttt                                                                                                     <  2:4292093/150‑1 (MQ=255)
                                                   gcgcTTCGTCGGTGTCTGCACTGAGCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCTAttttt                                                                                                    <  2:8235927/150‑1 (MQ=255)
                                                   gcgcTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAAttttt                                                                                                    >  1:2917784/1‑150 (MQ=255)
                                                     gcTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTAc                                                                                                  >  2:4049168/1‑150 (MQ=255)
                                                     gcTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTAc                                                                                                  >  1:2749632/1‑150 (MQ=255)
                                                       ttCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCa                                                                                                >  1:4551183/1‑150 (MQ=255)
                                                          gtcgGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGcc                                                                                             >  2:5571378/1‑150 (MQ=255)
                                                          gtcgGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGcc                                                                                             <  1:8710958/150‑1 (MQ=255)
                                                          gtcgGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGcc                                                                                             <  2:2084590/150‑1 (MQ=255)
                                                           tcgGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAATGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCg                                                                                            >  2:7664522/1‑150 (MQ=255)
                                                              gTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCAGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAATGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTc                                                                                         >  1:8354756/1‑150 (MQ=255)
                                                                gACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCAc                                                                                       >  2:6174594/1‑150 (MQ=255)
                                                                 aCTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGTGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCAcc                                                                                      <  1:2610773/150‑1 (MQ=255)
                                                                 aCTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCAcc                                                                                      <  2:7808047/150‑1 (MQ=255)
                                                                   tGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGt                                                                                    >  2:6153597/1‑150 (MQ=255)
                                                                        tGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTccg                                                                               >  1:2414473/1‑150 (MQ=255)
                                                                        tGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTccg                                                                               <  2:7176182/150‑1 (MQ=255)
                                                                          aCCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTccgcc                                                                             <  1:3598287/150‑1 (MQ=255)
                                                                           ccGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTccgccg                                                                            >  1:4058447/1‑150 (MQ=255)
                                                                              gTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTa                                                                         >  2:603323/1‑150 (MQ=255)
                                                                                    ccGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTAt                                                                   >  2:8809053/1‑150 (MQ=255)
                                                                                     cGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTAtt                                                                  >  2:5952337/1‑150 (MQ=255)
                                                                                     cGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTAtt                                                                  >  2:166630/1‑150 (MQ=255)
                                                                                      gaagaaGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAATGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTTCCAGCCGTTCACCGTTTCCGCCGCTACCGTATTc                                                                 >  1:5975925/1‑150 (MQ=255)
                                                                                      gaagaaGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTc                                                                 >  2:2182275/1‑150 (MQ=255)
                                                                                       aagaagAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTcg                                                                <  2:4880005/150‑1 (MQ=255)
                                                                                       aagaagAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTcg                                                                >  2:2868901/1‑150 (MQ=255)
                                                                                         gaagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTcgcg                                                              >  1:7366205/1‑150 (MQ=255)
                                                                                          aagaGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCgg                                                             <  2:5329612/150‑1 (MQ=255)
                                                                                            gagCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAGGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGt                                                           >  2:4441197/1‑150 (MQ=255)
                                                                                              gCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGAGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTaa                                                         >  1:5857427/1‑150 (MQ=255)
                                                                                                         gATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTa                                              <  2:2567095/150‑1 (MQ=255)
                                                                                                         gATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTa                                              >  2:3184110/1‑150 (MQ=255)
                                                                                                         gATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTa                                              <  1:3444513/150‑1 (MQ=255)
                                                                                                          aTAGTGGGTTGCCGTTGAACAATGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACGGGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTac                                             >  2:3919698/1‑150 (MQ=255)
                                                                                                           tAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTaca                                            <  1:1469692/150‑1 (MQ=255)
                                                                                                                 gTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATc                                      <  2:8158876/150‑1 (MQ=255)
                                                                                                                 gTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATc                                      <  1:6830942/150‑1 (MQ=255)
                                                                                                                    gCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCggg                                   >  1:4310940/1‑150 (MQ=255)
                                                                                                                     ccGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGgga                                  <  1:1632800/150‑1 (MQ=255)
                                                                                                                     ccGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGgga                                  <  2:5364840/150‑1 (MQ=255)
                                                                                                                     ccGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGgga                                  >  2:1732760/1‑150 (MQ=255)
                                                                                                                     ccGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGgga                                  <  1:5249208/150‑1 (MQ=255)
                                                                                                                      gtaaaatcAATGGCCCGGGGCTTGTGAAAGCTCCGGGGCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAATCCCGCACCTACACTGATCGggag                                 >  2:1595752/8‑150 (MQ=25)
                                                                                                                      cGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                                 <  1:5550304/150‑1 (MQ=255)
                                                                                                                      cGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                                 <  1:4792218/150‑1 (MQ=255)
                                                                                                                      cGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggag                                 <  1:6236808/150‑1 (MQ=255)
                                                                                                                       gTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggagg                                >  2:7430336/1‑150 (MQ=255)
                                                                                                                        ttGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggagga                               <  2:3800518/150‑1 (MQ=255)
                                                                                                                        ttGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGggagga                               <  1:1943384/150‑1 (MQ=255)
                                                                                                                           aaCAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGac                            >  1:8165718/1‑150 (MQ=21)
                                                                                                                           aaCAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGac                            >  1:1317674/1‑150 (MQ=21)
                                                                                                                            aCAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGaca                           <  1:8530865/150‑1 (MQ=255)
                                                                                                                            aCAAAGGCCCGGAGCTTTCACAAGCCCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGaca                           >  1:5730410/1‑150 (MQ=25)
                                                                                                                             cAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGacac                          <  1:1096376/150‑1 (MQ=255)
                                                                                                                              aaaGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGacaca                         >  1:4210045/1‑150 (MQ=255)
                                                                                                                                aGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGAcacaca                       <  1:6611530/150‑1 (MQ=255)
                                                                                                                                 ggCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCAGCAACCGTATTAGCGGGTAAAAAAGCAACTAAACTGAACGGGAGGAGAaacacaa                      >  1:928319/1‑149 (MQ=25)
                                                                                                                                   cccGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCg                    >  2:3351600/1‑150 (MQ=255)
                                                                                                                                    ccGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCgg                   >  1:3167906/1‑150 (MQ=255)
                                                                                                                                    ccGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCgg                   >  1:8010494/1‑150 (MQ=255)
                                                                                                                                    ccGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCgg                   >  2:8437665/1‑150 (MQ=255)
                                                                                                                                    ccGGGGCGTGTGACAGCTCCGGGCCGTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTCAACCCGCCCCGACACTGATCGGGAGGAGACACACAGCgg                   >  2:1859060/1‑150 (MQ=25)
                                                                                                                                       gggCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGa                <  2:6965842/150‑1 (MQ=255)
                                                                                                                                         gCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGAtt              <  2:5762878/150‑1 (MQ=255)
                                                                                                                                           ttGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTgg            <  1:7429554/150‑1 (MQ=255)
                                                                                                                                             tcacaaGCCCCGGGCCTTTGTTTTTACCCGAGAAGGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGgc          <  2:7554083/146‑1 (MQ=255)
                                                                                                                                               gAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGgcgc        >  2:6358759/1‑150 (MQ=255)
                                                                                                                                               gAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGgcgc        >  1:5763947/1‑150 (MQ=255)
                                                                                                                                                aaaGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCa       <  2:1081732/150‑1 (MQ=255)
                                                                                                                                                     tCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCg  >  1:1257816/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGAAGTGGGCACTGGCTACTTCGACGACATGACCACCGTGATCCAGGGCGGCGCTTCGTCGGTGACTGCACTGACCGGTTCGACCGAAGAAGAGCAGTTCCACTGATAGTGGGTTGCCGTTGAACAAAGGCCCGGGGCTTGTGAAAGCTCCGGGCCTTTGTTTTTACCCGAGAAAGGGCCGGGCCTGCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGATCGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCG  >  NC_002947/4659737‑4660035

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