Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,986,911 G→T 5.4% Q63K (CAA→AAA) ‡ ycaC‑II ← putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,986,9110GT 5.4% 1129.7 / 34.5 334Q63K (CAA→AAA) ‡ycaC‑IIputative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (132/183);  new base T (9/9);  total (142/192)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 335 total aligned reads displayed.

GGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGCG  >  NC_002947/5986762‑5987059
                                                                                                                                                     |                                                                                                                                                    
   gatgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTc                                                                                                                                                   >  1:1933730/1‑150 (MQ=255)
   gatgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTc                                                                                                                                                   >  2:4436954/1‑150 (MQ=255)
     tgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGa                                                                                                                                                 >  1:4338227/1‑150 (MQ=255)
      gatGAGTTGCGTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGTTCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCTGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGaa                                                                                                                                                <  1:6408166/150‑1 (MQ=255)
      gatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGaa                                                                                                                                                <  1:4621951/150‑1 (MQ=255)
         cAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCt                                                                                                                                             >  2:1753487/1‑150 (MQ=255)
          aGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtg                                                                                                                                            <  1:9010169/150‑1 (MQ=255)
            ttgcttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggt                                                                                                                                          >  1:6992440/1‑150 (MQ=255)
             tgcttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtg                                                                                                                                         <  1:8182542/150‑1 (MQ=255)
             tgcttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtg                                                                                                                                         <  1:900472/150‑1 (MQ=255)
                ttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGgtggtg                                                                                                                                      >  2:4085715/1‑150 (MQ=255)
                ttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGgtggtg                                                                                                                                      >  1:4935903/1‑150 (MQ=255)
                  gcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAg                                                                                                                                    >  2:8487287/1‑150 (MQ=255)
                  gcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAg                                                                                                                                    <  1:6007703/150‑1 (MQ=255)
                   cGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAgg                                                                                                                                   <  1:4775771/150‑1 (MQ=255)
                    ggCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGa                                                                                                                                  <  1:3002197/150‑1 (MQ=255)
                    ggCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGa                                                                                                                                  <  1:1012923/150‑1 (MQ=255)
                    ggCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGa                                                                                                                                  >  2:4314046/1‑150 (MQ=255)
                      ccggTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACTGGTGGTCAgaaca                                                                                                                                <  2:3273256/147‑5 (MQ=2)
                        ggTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACTGGTGGTCAgaacagt                                                                                                                              <  1:1514744/149‑7 (MQ=2)
                        cGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                                                                                                                              >  1:1718980/1‑150 (MQ=255)
                         gTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTc                                                                                                                             <  2:4940493/150‑1 (MQ=255)
                            gCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGc                                                                                                                          <  2:1993407/150‑1 (MQ=255)
                              cTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                                                                                                                        >  1:7948562/1‑150 (MQ=255)
                              cTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGAAATGTACGGCGCATCCGGGAACAACTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTAGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                                                                                                                        >  2:5875499/1‑150 (MQ=255)
                               ttGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAgg                                                                                                                       >  2:4760162/1‑150 (MQ=255)
                               ttGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAgg                                                                                                                       <  1:6642685/150‑1 (MQ=255)
                                    tGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCaaa                                                                                                                  <  2:3541212/150‑1 (MQ=255)
                                    tGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCaaa                                                                                                                  <  2:5505722/150‑1 (MQ=255)
                                    tGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCaaa                                                                                                                  <  2:36644/150‑1 (MQ=255)
                                       cTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGaa                                                                                                               <  2:1315379/150‑1 (MQ=255)
                                        ttGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAAc                                                                                                              <  2:6461808/150‑1 (MQ=255)
                                         tGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACt                                                                                                             >  1:330054/1‑150 (MQ=255)
                                          gACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACtt                                                                                                            <  1:7500451/150‑1 (MQ=255)
                                           ccgcaGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg                                                                                                           <  2:6886473/146‑1 (MQ=255)
                                           aCGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg                                                                                                           <  2:4223470/150‑1 (MQ=255)
                                           aCGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg                                                                                                           <  2:6146069/150‑1 (MQ=255)
                                             gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc                                                                                                         >  1:9094824/1‑150 (MQ=255)
                                             gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc                                                                                                         >  1:5092242/1‑150 (MQ=255)
                                             gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc                                                                                                         >  2:2820749/1‑150 (MQ=255)
                                                gTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAg                                                                                                      <  2:4484197/150‑1 (MQ=255)
                                                 tCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAgg                                                                                                     <  2:3170462/150‑1 (MQ=255)
                                                   ttCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTc                                                                                                   >  1:6226597/1‑150 (MQ=255)
                                                   ttCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTc                                                                                                   >  1:7053919/1‑150 (MQ=255)
                                                         gTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCaa                                                                                             <  2:5066801/150‑1 (MQ=255)
                                                            ccAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCCCCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGc                                                                                          >  1:1978544/1‑150 (MQ=255)
                                                            ccAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGc                                                                                          <  2:7614645/150‑1 (MQ=255)
                                                             cAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGcc                                                                                         >  1:8608601/1‑150 (MQ=255)
                                                              aGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCa                                                                                        >  1:8752412/1‑150 (MQ=255)
                                                                  gTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGGCAGGTCGCCCAAGGCCAGCAc                                                                                    >  2:5689867/1‑150 (MQ=255)
                                                                  gTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCAc                                                                                    >  2:7585457/1‑150 (MQ=255)
                                                                       tCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAg                                                                                                                                               >  1:1307859/1‑86 (MQ=255)
                                                                       tCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACgttgt                                                                               >  1:8493549/1‑150 (MQ=255)
                                                                          ggCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCt                                                                            <  1:5505247/150‑1 (MQ=255)
                                                                          ggCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCt                                                                            <  2:3072791/150‑1 (MQ=255)
                                                                           gCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCtt                                                                           >  2:6002322/1‑150 (MQ=255)
                                                                           gCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCtt                                                                           <  2:6696465/150‑1 (MQ=255)
                                                                           gCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCtt                                                                           <  2:7038244/150‑1 (MQ=255)
                                                                              gggACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGaa                                                                        >  1:3513577/1‑150 (MQ=255)
                                                                              gggACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGaa                                                                        <  2:2631803/150‑1 (MQ=255)
                                                                                   gggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGt                                                                   >  1:4528631/1‑150 (MQ=255)
                                                                                   gggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGt                                                                   >  1:5191274/1‑150 (MQ=255)
                                                                                    ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                                                                  <  2:9094290/150‑1 (MQ=255)
                                                                                    ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                                                                  >  1:888499/1‑150 (MQ=255)
                                                                                    ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                                                                  >  1:6966990/1‑150 (MQ=255)
                                                                                    ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                                                                  <  1:6539271/150‑1 (MQ=255)
                                                                                      cAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTgg                                                                >  2:6079051/1‑150 (MQ=255)
                                                                                      cAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTgg                                                                >  2:5122777/1‑150 (MQ=255)
                                                                                        aTGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGc                                                              <  2:391324/150‑1 (MQ=255)
                                                                                        aTGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGc                                                              <  1:5313367/150‑1 (MQ=255)
                                                                                            aCGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTACTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGaaa                                                          >  1:7261933/1‑150 (MQ=255)
                                                                                               gcgcATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAAt                                                       >  1:4249354/1‑150 (MQ=255)
                                                                                                cgcATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATc                                                      <  1:5878592/150‑1 (MQ=255)
                                                                                                cgcATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATc                                                      >  2:5356778/1‑150 (MQ=255)
                                                                                                 gcATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATcc                                                     <  2:1384439/150‑1 (MQ=255)
                                                                                                  cATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCt                                                    <  1:1395731/150‑1 (MQ=255)
                                                                                                  cATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCt                                                    <  1:2422193/150‑1 (MQ=255)
                                                                                                   aTCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTg                                                   <  1:4436954/150‑1 (MQ=255)
                                                                                                     ccGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCa                                                 <  2:8136415/150‑1 (MQ=255)
                                                                                                      cGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCAc                                                <  1:2820749/150‑1 (MQ=255)
                                                                                                        ggAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCa                                              >  2:5639370/1‑150 (MQ=255)
                                                                                                              tctcTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCgaga                                        <  2:6386663/150‑1 (MQ=255)
                                                                                                              tctcTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCgaga                                        <  2:6580975/150‑1 (MQ=255)
                                                                                                               ctctTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGAt                                       >  2:8292783/1‑150 (MQ=255)
                                                                                                                  tttCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAg                                    >  1:2129509/1‑150 (MQ=255)
                                                                                                                  tttCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAg                                    >  1:2650938/1‑150 (MQ=255)
                                                                                                                   ttCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGa                                   <  1:3125624/150‑1 (MQ=255)
                                                                                                                      aGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCg                                >  1:3780015/1‑150 (MQ=255)
                                                                                                                        cTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGt                              <  2:1609871/150‑1 (MQ=255)
                                                                                                                             gCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCtggt                         <  2:215663/150‑1 (MQ=255)
                                                                                                                              cACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCtggtg                        >  1:3020130/1‑150 (MQ=255)
                                                                                                                              cACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCtggtg                        <  1:783113/150‑1 (MQ=255)
                                                                                                                                 cAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTc                     <  2:1463753/150‑1 (MQ=255)
                                                                                                                                  aGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCga                                          <  2:2436160/128‑1 (MQ=255)
                                                                                                                                   gCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGa                   <  2:6226597/150‑1 (MQ=255)
                                                                                                                                          aTTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACcagcag            >  2:7739721/1‑150 (MQ=255)
                                                                                                                                                  cTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACgg    >  2:7177695/1‑150 (MQ=255)
                                                                                                                                                  cTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACgg    <  2:129778/150‑1 (MQ=255)
                                                                                                                                                  cTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACgg    >  1:1594712/1‑150 (MQ=255)
                                                                                                                                                   ttGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGc   <  2:8762505/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGAGATCAGACCGGTCTGGTGGTCGACCAGCAGCACCACGGCG  >  NC_002947/5986762‑5987059

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: