Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,594,108 G→T 5.1% intergenic (‑264/+37) PP_1399 ← / ← kgtP hypothetical protein/alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,594,1080GT 5.1% 593.1 / 15.0 177intergenic (‑264/+37)PP_1399/kgtPhypothetical protein/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base G (83/85);  new base T (2/7);  total (85/92)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 177 total aligned reads displayed.

TGGGCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCG  >  NC_002947/1593961‑1594256
                                                                                                                                                   |                                                                                                                                                    
         aTGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCATTGCATGGAgcgc                                                                                                                                           <  2:3629271/150‑1 (MQ=255)
           gTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGccc                                                                                                                                         >  2:1370694/1‑150 (MQ=255)
           gTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGccc                                                                                                                                         >  2:4406932/1‑150 (MQ=255)
             cGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCt                                                                                                                                       <  1:1608515/150‑1 (MQ=255)
             cGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCt                                                                                                                                       <  1:1608668/150‑1 (MQ=255)
               tGCTTCAGGCCGCACGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCttt                                                                                                                                     <  1:3511288/150‑1 (MQ=255)
               tGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCttt                                                                                                                                     >  1:1205565/1‑150 (MQ=255)
                 cTTCAGGCCGCACGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTa                                                                                                                                   <  2:3142210/150‑1 (MQ=255)
                  ttCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAAGTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCACTGGGGCGCCCCtgtgag                                                                                                                                  >  2:3779632/1‑145 (MQ=255)
                      ggCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCCCCATTGCAGGGAGCGCCCCttttctattc                                                                                                                              <  2:4404566/150‑7 (MQ=38)
                        ccGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTgc                                                                                                                            <  1:2275000/150‑1 (MQ=255)
                         cGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTgcc                                                                                                                           >  1:237577/1‑150 (MQ=255)
                           cAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTgccag                                                                                                                         >  2:78356/1‑150 (MQ=255)
                               gctggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAATGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTgccaggcca                                                                                                                     >  1:3162544/1‑150 (MQ=255)
                                 tggctTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGc                                                                                                                   >  2:3389044/1‑150 (MQ=255)
                                            cGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                        <  1:3482/150‑1 (MQ=255)
                                            cGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTg                                                                                                        <  2:4949827/150‑1 (MQ=255)
                                                  cgcgGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATc                                                                                                  <  1:2256221/150‑1 (MQ=255)
                                                         gcCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGt                                                                                           >  2:3818119/1‑150 (MQ=255)
                                                           cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATgg                                                                                                                                                >  2:700818/1‑95 (MQ=255)
                                                           cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATgg                                                                                                                                                <  1:700818/95‑1 (MQ=255)
                                                           cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATcc                                                                                                 >  2:4065647/1‑142 (MQ=255)
                                                           cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATcc                                                                                                 <  1:4065647/142‑1 (MQ=255)
                                                           cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgc                                                                                         >  1:4725389/1‑150 (MQ=255)
                                                           cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgc                                                                                         >  1:2733512/1‑150 (MQ=255)
                                                           cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgc                                                                                         >  2:3508995/1‑150 (MQ=255)
                                                           cAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGtgc                                                                                         >  2:3438106/1‑150 (MQ=255)
                                                             aCGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTa                                                                                                               <  2:940981/126‑1 (MQ=255)
                                                             aCGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTa                                                                                                               >  1:940981/1‑126 (MQ=255)
                                                               gTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTgctgct                                                                                     >  2:2824029/1‑150 (MQ=255)
                                                                 tCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTg                                                                                   <  2:1817059/150‑1 (MQ=255)
                                                                   tGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgt                                                                                 <  1:1701605/150‑1 (MQ=255)
                                                                   tGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgt                                                                                 >  1:425786/1‑150 (MQ=255)
                                                                     ccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgtgt                                                                               >  2:4072873/1‑150 (MQ=255)
                                                                     ccagccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGgtgt                                                                               >  2:1672105/1‑150 (MQ=255)
                                                                       agccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCt                                                                             >  2:4356166/1‑150 (MQ=255)
                                                                        gccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCtt                                                                            <  2:5331103/150‑1 (MQ=255)
                                                                        gccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCtt                                                                            >  2:3010656/1‑150 (MQ=255)
                                                                         ccaTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCttt                                                                           <  2:4881771/150‑1 (MQ=255)
                                                                           aTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTACGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCa                                                                         <  2:2522091/150‑1 (MQ=255)
                                                                           aTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTcagatc                                                                         >  1:4923066/1‑145 (MQ=255)
                                                                           aTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCa                                                                         <  1:1094859/150‑1 (MQ=255)
                                                                           aTTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCa                                                                         <  1:5240888/150‑1 (MQ=255)
                                                                                aGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTc                                                                    <  1:4979395/150‑1 (MQ=255)
                                                                                  cTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCtt                                                                                                                                      >  1:2922830/1‑82 (MQ=255)
                                                                                  cTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCtt                                                                                                                                      <  2:2922830/82‑1 (MQ=255)
                                                                                  cTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTcgc                                                                  <  2:1002721/150‑1 (MQ=255)
                                                                                  cTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTcgc                                                                  >  2:1620559/1‑150 (MQ=255)
                                                                                     gATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTGTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTa                                                               >  2:580755/1‑150 (MQ=255)
                                                                                      aTAAGGTGGAATTAATAGGGCGTTAAGAAGGAAATGTCTGTATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTaa                                                              >  2:475447/1‑150 (MQ=255)
                                                                                       tAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTaaa                                                             <  1:1660492/150‑1 (MQ=255)
                                                                                           gTGGAATTAATGGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTa                                                                                                                                   <  1:4334645/76‑1 (MQ=255)
                                                                                           gTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcag                                                         >  2:1999146/1‑150 (MQ=255)
                                                                                             ggAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCGTGGGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagca                                                       <  2:5308207/150‑1 (MQ=255)
                                                                                             ggAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagca                                                       <  2:4210021/150‑1 (MQ=255)
                                                                                              gAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACcagcag                                                      >  2:2372355/1‑150 (MQ=255)
                                                                                               aaTTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCagcagc                                                     <  2:3741912/150‑1 (MQ=255)
                                                                                                aTTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCg                                                    >  1:3164866/1‑150 (MQ=255)
                                                                                                 ttAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGa                                                   <  2:1711599/150‑1 (MQ=255)
                                                                                                       gggCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGc                                             <  1:626858/150‑1 (MQ=255)
                                                                                                         gCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGa                                           <  2:4674391/150‑1 (MQ=255)
                                                                                                         gCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGa                                           >  2:5312249/1‑150 (MQ=255)
                                                                                                         gCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGa                                           >  2:328984/1‑150 (MQ=255)
                                                                                                           gTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCTATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATa                                         <  2:1315564/150‑1 (MQ=255)
                                                                                                           gTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATa                                         <  2:3058273/150‑1 (MQ=255)
                                                                                                           gTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATa                                         >  1:3043704/1‑150 (MQ=255)
                                                                                                           gTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATa                                         <  2:575416/150‑1 (MQ=255)
                                                                                                             tAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACa                                       <  1:378466/150‑1 (MQ=255)
                                                                                                                gaagGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGc                                    >  1:3112572/1‑150 (MQ=255)
                                                                                                                  agGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCgg                                  <  2:1344056/150‑1 (MQ=255)
                                                                                                                   ggAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt                                 >  2:282718/1‑150 (MQ=255)
                                                                                                                   ggAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt                                 >  1:1449144/1‑150 (MQ=255)
                                                                                                                   ggAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGt                                 >  1:2006362/1‑150 (MQ=255)
                                                                                                                    gAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTg                                <  1:105541/150‑1 (MQ=255)
                                                                                                                    gAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTg                                >  2:2202893/1‑150 (MQ=255)
                                                                                                                    gAAAAGTCTGAAGAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTg                                <  2:2406893/150‑1 (MQ=255)
                                                                                                                         gTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTa                           <  2:2073065/150‑1 (MQ=255)
                                                                                                                          tCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                          >  2:5203984/1‑150 (MQ=255)
                                                                                                                          tCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTAc                          <  1:1246770/150‑1 (MQ=255)
                                                                                                                            tGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCa                        >  2:5377324/1‑150 (MQ=255)
                                                                                                                            tGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCa                        <  1:1922809/150‑1 (MQ=255)
                                                                                                                            tGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCa                        <  1:3389044/150‑1 (MQ=255)
                                                                                                                             gAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAg                       >  1:1206607/1‑150 (MQ=255)
                                                                                                                             gAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAg                       <  1:2951719/150‑1 (MQ=255)
                                                                                                                              aaTAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGt                      >  1:2755416/1‑150 (MQ=255)
                                                                                                                              aaTAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGt                      >  1:2282485/1‑150 (MQ=255)
                                                                                                                                   aaaaGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAg                 >  1:5009861/1‑150 (MQ=255)
                                                                                                                                   aaaaGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAg                 >  2:5044078/1‑150 (MQ=255)
                                                                                                                                    aaaGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGc                <  1:4356166/150‑1 (MQ=255)
                                                                                                                                        gggCGCCCCATTGCATGGAGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGa                                           >  1:1714009/1‑119 (MQ=39)
                                                                                                                                          gcgcTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGctct          >  2:3612627/1‑150 (MQ=255)
                                                                                                                                          gcgcTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGctct          >  1:2028118/1‑150 (MQ=255)
                                                                                                                                              tCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCAt      >  2:629186/1‑150 (MQ=255)
                                                                                                                                               ccATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATc     <  2:1205912/150‑1 (MQ=255)
                                                                                                                                               ccATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATc     <  2:1205565/150‑1 (MQ=255)
                                                                                                                                                cATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATcc    >  1:2083186/1‑150 (MQ=255)
                                                                                                                                                cATGCAATGGGGAGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCAtgg    <  2:3162544/150‑3 (MQ=255)
                                                                                                                                                 aTGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCAt                                                                        >  1:553104/1‑81 (MQ=255)
                                                                                                                                                  tGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCg  <  1:1370694/150‑1 (MQ=255)
                                                                                                                                                  tGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCg  >  1:2496624/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGGCGCGTATGTCGTGCTTCAGGCCGCAAGGCTGGCTTGGGAGCGTTCGCGCGGGCGCCAACGTTCTGCCAGCCATTGCAGCTTGATAAGGTGGAATTAATAGGGCGTTAAGAAGGAAAAGTCTGAATAGAAAAGGGGCGCTCCATGCAATGGGGCGCCCCTTTTAGTACTGCCAGGCCAGCGTTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCG  >  NC_002947/1593961‑1594256

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