Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 264,008 C→G 5.2% S253S (TCC→TCG) ‡ gabD‑I → succinate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947264,0080CG 5.2% 581.7 / 15.1 175S253S (TCC→TCG) ‡gabD‑Isuccinate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (85/80);  new base G (3/6);  total (89/86)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 181 total aligned reads displayed.

GCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGA  >  NC_002947/263860‑264157
                                                                                                                                                    |                                                                                                                                                     
 cGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCt                                                                                                                                                     <  2:5349578/150‑1 (MQ=255)
   tGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTgg                                                                                                                                                   >  1:2854545/1‑150 (MQ=255)
    gCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGa                                                                                                                                                  <  2:1039867/150‑1 (MQ=255)
      tGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGc                                                                                                                                                >  2:5398953/1‑150 (MQ=255)
        aGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTg                                                                                                                                              >  2:5181560/1‑150 (MQ=255)
         gtgtGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTgg                                                                                                                                             <  2:1048983/150‑1 (MQ=255)
            tgGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTg                                                                                                                                          >  2:2666721/1‑150 (MQ=255)
             ggTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGACGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTgg                                                                                                                                         >  1:134920/1‑150 (MQ=255)
                 aCCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAAc                                                                                                                                     <  1:4407849/150‑1 (MQ=255)
                      cAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGcccc                                                                                                                                <  1:326869/150‑1 (MQ=255)
                       aGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCg                                                                                                                               >  1:3427251/1‑150 (MQ=255)
                            cGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCAt                                                                                                                          >  2:4319564/1‑150 (MQ=255)
                             ggCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATc                                                                                                                         <  2:2989576/150‑1 (MQ=255)
                             ggCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATc                                                                                                                         >  1:4850149/1‑150 (MQ=255)
                             ggCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATc                                                                                                                         <  2:1612661/150‑1 (MQ=255)
                                gAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTgg                                                                                                                                             >  1:345764/1‑127 (MQ=255)
                                gAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCgtg                                                                                                                      >  1:1800453/1‑150 (MQ=255)
                                gAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCgtg                                                                                                                      >  1:42430/1‑150 (MQ=255)
                                     tGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTcga                                                                                                                 >  2:5275784/1‑150 (MQ=255)
                                        cggcgAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTcgacga                                                                                                              >  2:3591813/1‑150 (MQ=255)
                                         ggcgAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCgacgac                                                                                                             >  2:1238558/1‑150 (MQ=255)
                                            gAACTGACAGGCAACTCCCTGGTACGCAAGATGTCCTTCACCGGCTCGACCGAAATCGGTCGCAAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCAGTTCATCGTGTTCGACGACGcc                                                                                                          >  2:3401791/1‑150 (MQ=255)
                                                tGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGAcc                                                                                                      >  1:4772501/1‑150 (MQ=255)
                                                tGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGAcc                                                                                                      >  2:4089191/1‑150 (MQ=255)
                                                tGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGAcc                                                                                                      >  2:4089199/1‑150 (MQ=255)
                                                tGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGAcc                                                                                                      >  2:4871576/1‑150 (MQ=255)
                                                 gACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGACGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCt                                                                                                     <  2:2883686/150‑1 (MQ=255)
                                                           tCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCg                                                                                           >  1:1090239/1‑150 (MQ=255)
                                                           tCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCg                                                                                           <  1:2380585/150‑1 (MQ=255)
                                                           tCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCg                                                                                           >  2:100638/1‑150 (MQ=255)
                                                           tCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCg                                                                                           >  1:943371/1‑150 (MQ=255)
                                                           tCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCg                                                                                           >  1:1090264/1‑150 (MQ=255)
                                                           tCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCg                                                                                           <  1:2850392/150‑1 (MQ=255)
                                                             ccTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGt                                                                                         <  2:2159488/150‑1 (MQ=255)
                                                              cTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTc                                                                                        >  2:38943/1‑150 (MQ=255)
                                                                ggTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGa                                                                                      <  2:4990115/150‑1 (MQ=255)
                                                                  tACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAgg                                                                                    <  2:1332077/150‑1 (MQ=255)
                                                                  tACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAgg                                                                                    <  1:5343723/150‑1 (MQ=255)
                                                                     gCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGATGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGgcg                                                                                 <  2:3176330/150‑1 (MQ=255)
                                                                         gCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGat                                                                             >  1:2261091/1‑150 (MQ=255)
                                                                          cTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGatc                                                                            <  1:4607466/150‑1 (MQ=255)
                                                                           ggTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGttca                                                                           >  1:398575/2‑150 (MQ=255)
                                                                             tCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGc                                                                                                           <  2:790001/116‑1 (MQ=255)
                                                                             tCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGc                                                                                                           >  1:790001/1‑116 (MQ=255)
                                                                             tCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGatcatc                                                                         >  1:3074263/1‑150 (MQ=255)
                                                                               cTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCAtctc                                                                       >  2:4665093/1‑150 (MQ=255)
                                                                               cTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCAtctc                                                                       >  1:2643152/1‑150 (MQ=255)
                                                                                  ctccGGCTCGACCGAAATCGGTCGCCAGCTGATGGATGTATGCGCCAAGGACTTCAATCAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATGGTGTTCTTCGACGCCGACCTGGACAAGGCGGTAGAGGGCGCGATCTTCTCCaa                                                                    <  1:4904290/148‑1 (MQ=255)
                                                                                  cACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCaa                                                                    <  2:3074263/150‑1 (MQ=255)
                                                                                     cGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTa                                                                 >  1:2914906/1‑150 (MQ=255)
                                                                                     cGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTa                                                                 <  2:684735/150‑1 (MQ=255)
                                                                                            aCCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTaac                                                          >  1:1978527/1‑150 (MQ=255)
                                                                                            aCCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTaac                                                          <  1:45108/150‑1 (MQ=255)
                                                                                              cGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGGCGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTaacaa                                                        >  1:1826817/1‑150 (MQ=255)
                                                                                              cGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTaacaa                                                        >  2:3581730/1‑150 (MQ=255)
                                                                                                 tttCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACaa                                                                                               >  2:481543/3‑108 (MQ=255)
                                                                                                 tttCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACaa                                                                                               <  1:481543/106‑1 (MQ=255)
                                                                                                      gTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGa                                                >  2:5397655/1‑150 (MQ=255)
                                                                                                      gTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGa                                                <  1:3675275/150‑1 (MQ=255)
                                                                                                      gTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGa                                                >  2:5052712/1‑150 (MQ=255)
                                                                                                       tCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACACCgg                                                     >  2:4419445/1‑144 (MQ=255)
                                                                                                       tCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACgg                                                     <  1:4419445/144‑1 (MQ=255)
                                                                                                       tCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGAc                                               >  1:5392053/1‑150 (MQ=255)
                                                                                                         gCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCt                                             <  1:4048770/150‑1 (MQ=255)
                                                                                                         gCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCt                                             <  2:4637124/150‑1 (MQ=255)
                                                                                                           cAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGc                                           <  2:3567857/150‑1 (MQ=255)
                                                                                                           cAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGc                                           <  1:1914069/150‑1 (MQ=255)
                                                                                                           cAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGc                                           >  2:4392832/1‑150 (MQ=255)
                                                                                                           cAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGc                                           <  2:678214/150‑1 (MQ=255)
                                                                                                           cAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGc                                           <  1:5181560/150‑1 (MQ=255)
                                                                                                            aGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCg                                          >  1:4064585/1‑150 (MQ=255)
                                                                                                              cTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTc                                        <  2:3767615/150‑1 (MQ=255)
                                                                                                              cTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTc                                        <  2:4478204/150‑1 (MQ=255)
                                                                                                                gATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTg                                      >  2:2278282/1‑150 (MQ=255)
                                                                                                                 aTGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTgc                                     <  1:3671113/150‑1 (MQ=255)
                                                                                                                   ggAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTgcgc                                   >  1:2365814/1‑150 (MQ=255)
                                                                                                                   ggAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTgcgc                                   >  1:810993/1‑150 (MQ=255)
                                                                                                                       gaaTGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAAc                               <  2:134920/150‑1 (MQ=255)
                                                                                                                            cgcCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTAt                          <  1:5061710/150‑1 (MQ=255)
                                                                                                                              cgccgaccATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTctct                        <  1:2777033/143‑1 (MQ=255)
                                                                                                                              cgccgaccATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTctct                        <  1:2777218/143‑1 (MQ=255)
                                                                                                                               gccgaccATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTa                       <  1:3930155/144‑1 (MQ=255)
                                                                                                                               cAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTa                       >  1:1296172/1‑150 (MQ=255)
                                                                                                                                   gACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGTTCATCTCAAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTg                   >  1:1686754/1‑150 (MQ=255)
                                                                                                                                    aCATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGc                  >  2:186261/1‑150 (MQ=255)
                                                                                                                                     cATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACa                 >  1:1110338/1‑150 (MQ=255)
                                                                                                                                     cATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACGGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACa                 >  1:2092476/1‑150 (MQ=255)
                                                                                                                                     cATCAAGAAGCTTGGGCTTGAGCTTGGCGGCAACGCGCCTTTCATCGTGTTCGACGACGCCGATCTCGATGCCGCCGTCGAAGGTGCTCTTATCGCCAAATTCCGCAACAACAGCCAAACCTGCGTCTGCGCCAACCGTCTCTATGTACa                 <  2:713263/150‑1 (MQ=255)
                                                                                                                                      aTCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGc                                                    >  1:1538259/1‑114 (MQ=255)
                                                                                                                                      aTCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAg                <  2:3523272/150‑1 (MQ=255)
                                                                                                                                       tCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAgg               >  1:671489/1‑150 (MQ=255)
                                                                                                                                         aagaagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGAc             >  1:3197777/1‑150 (MQ=255)
                                                                                                                                         aagaagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGAc             >  1:1757191/1‑150 (MQ=255)
                                                                                                                                          agaagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACg            >  2:1145581/1‑150 (MQ=255)
                                                                                                                                           gaagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACgg           >  1:2267641/1‑150 (MQ=255)
                                                                                                                                           gaagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACgg           >  2:2024866/1‑150 (MQ=255)
                                                                                                                                            aagGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGc          <  2:1436971/150‑1 (MQ=255)
                                                                                                                                                    ccTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTacga  >  2:4617820/1‑150 (MQ=255)
                                                                                                                                                    ccTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTacga  >  1:3823558/1‑150 (MQ=255)
                                                                                                                                                    ccTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTacga  <  1:4495837/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGTGCTGAGTGTGGTTACCGGCAGCGCCGGCGAAGTTGGCGGCGAACTGACCGGCAACTCCCTGGTACGCAAGCTGTCCTTCACCGGCTCGACCGAAATCGGTCGCCAGCTGATGGAAGAATGCGCCAAGGACATCAAGAAGGTTTCCCTGGAGCTGGGTGGCAACGCCCCGTTCATCGTGTTCGACGACGCCGACCTGGACAAGGCGGTCGAGGGCGCGATCATCTCCAAGTACCGTAACAACGGCCAGACCTGCGTCTGCGCCAACCGTATCTACGTGCAGGACGGCGTCTACGA  >  NC_002947/263860‑264157

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: