Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,992,420 G→T 7.8% A75A (GCC→GCA) ‡ ycaC‑III ← putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,992,4200GT 7.8% 1183.4 / 63.9 359A75A (GCC→GCA) ‡ycaC‑IIIputative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (157/174);  new base T (17/11);  total (174/185)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 367 total aligned reads displayed.

AGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGACAG  >  NC_002947/5992271‑5992569
                                                                                                                                                     |                                                                                                                                                     
aGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGGCGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAgg                                                                                                                                                       <  2:4833423/150‑1 (MQ=255)
aGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAgg                                                                                                                                                       <  1:1506157/150‑1 (MQ=255)
     cGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCttt                                                                                                                                                  >  1:1138613/1‑150 (MQ=255)
     cGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCttt                                                                                                                                                  <  1:3947672/150‑1 (MQ=255)
     cGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCttt                                                                                                                                                  >  1:2275011/1‑150 (MQ=255)
     cGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCttt                                                                                                                                                  <  2:7229895/150‑1 (MQ=255)
      gggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTc                                                                                                                                                 >  1:9156381/1‑150 (MQ=255)
       ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                                                                                                                                <  2:7759591/150‑1 (MQ=255)
       ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                                                                                                                                <  1:6861779/150‑1 (MQ=255)
            gCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTcc                                                                                                                                           >  2:6725806/1‑150 (MQ=255)
            gCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTcc                                                                                                                                           >  2:2655812/1‑150 (MQ=255)
             ccACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCg                                                                                                                                          >  1:751069/1‑150 (MQ=255)
                 gCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCGTCTTGCCGGTGGCCTTCCCTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCCTCCGGGAATGTTTCTTTCAGCTCCGGCAc                                                                                                                                      >  2:3345316/1‑150 (MQ=9)
                 gCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTATGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCAc                                                                                                                                      >  1:54893/1‑150 (MQ=255)
                   aaaCCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCa                                                                                                                                    >  1:3968986/1‑150 (MQ=255)
                     aCCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGc                                                                                                                                  >  1:263800/1‑150 (MQ=255)
                       cTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCgg                                                                                                                                <  2:2057225/150‑1 (MQ=255)
                           gTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGAGCTGGCCAGTGCGGGCAATGTACGGCGCATCCGGGACCAGGGCTTTCAGTTCCGGCACCAGCGGGCCg                                                                                                                            >  2:7510511/1‑150 (MQ=255)
                           gTCACAACACCTGCAATGATCAGCTGCGTCTTGCCGGTCGCCGTCCCCGCCTTGACGAAATCCTCGGTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCg                                                                                                                            <  2:761485/150‑1 (MQ=255)
                            tCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGccat                                                                                                                           >  1:1608516/1‑148 (MQ=9)
                              acaacaCCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCAttc                                                                                                                         >  1:830551/1‑149 (MQ=2)
                              acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg                                                                                                                         <  2:428245/150‑1 (MQ=255)
                              acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg                                                                                                                         >  2:8822862/1‑150 (MQ=255)
                              acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg                                                                                                                         >  2:4901561/1‑150 (MQ=255)
                              acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg                                                                                                                         <  2:5370899/150‑1 (MQ=255)
                               caacaCCAGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTgg                                                                                                                        <  1:5090160/150‑1 (MQ=255)
                                aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg                                                                                                                       >  2:4317169/1‑150 (MQ=255)
                                     cTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGccgt                                                                                                                  >  2:9074420/1‑148 (MQ=2)
                                      tGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGGTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTg                                                                                                                 <  2:1967914/150‑1 (MQ=255)
                                       gCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGc                                                                                                                <  2:3684844/150‑1 (MQ=255)
                                       gCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGc                                                                                                                <  2:570682/150‑1 (MQ=255)
                                            gATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGaa                                                                                                           <  1:6646911/150‑1 (MQ=255)
                                               cAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCt                                                                                                        <  1:3467338/150‑1 (MQ=255)
                                                 gctgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtgg                                                                                                      >  1:3060364/1‑150 (MQ=255)
                                                  ctgctTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACtggt                                                                                                     <  1:1236478/150‑1 (MQ=2)
                                                  ctgctTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACtggt                                                                                                     <  2:8657704/150‑1 (MQ=2)
                                                   tgctTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACtggtg                                                                                                    >  2:6240309/1‑150 (MQ=2)
                                                   tgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggtg                                                                                                    <  1:3402291/150‑1 (MQ=255)
                                                   tgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggtg                                                                                                    <  2:7114525/150‑1 (MQ=255)
                                                    gctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggtgg                                                                                                   >  2:2320322/1‑150 (MQ=255)
                                                      ttcttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGgtggtg                                                                                                 >  2:226758/1‑150 (MQ=255)
                                                       tcttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGa                                                                                                <  1:8504840/150‑1 (MQ=255)
                                                       tcttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGa                                                                                                <  2:6623920/150‑1 (MQ=255)
                                                        cttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAg                                                                                               >  2:165482/1‑150 (MQ=255)
                                                         ttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAgg                                                                                              >  1:6435098/1‑150 (MQ=255)
                                                          tGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGa                                                                                             <  2:246915/150‑1 (MQ=255)
                                                           gCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGAt                                                                                            <  2:1083690/150‑1 (MQ=255)
                                                            ccGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATg                                                                                           >  2:6472469/1‑150 (MQ=255)
                                                             cGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATgg                                                                                          >  1:6534527/1‑150 (MQ=255)
                                                              ggTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                                                                                         <  1:9089552/150‑1 (MQ=255)
                                                              ggTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                                                                                         <  1:3616052/150‑1 (MQ=255)
                                                                   ccTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCa                                                                                    >  2:4517527/1‑150 (MQ=255)
                                                                   ccTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCa                                                                                    <  2:2704960/150‑1 (MQ=255)
                                                                   ccTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCa                                                                                    <  2:7387101/150‑1 (MQ=255)
                                                                    cTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                                                                                   >  2:3768550/1‑150 (MQ=255)
                                                                     ttCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAgg                                                                                  >  1:714239/1‑150 (MQ=255)
                                                                        aCCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTg                                                                               <  1:7937848/150‑1 (MQ=255)
                                                                        aCCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTg                                                                               >  1:8857446/1‑150 (MQ=255)
                                                                         ccGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTga                                                                              >  1:2797852/1‑150 (MQ=255)
                                                                          cGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTgaa                                                                             <  1:5159536/150‑1 (MQ=255)
                                                                          cGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTgaa                                                                             <  1:6130871/150‑1 (MQ=255)
                                                                            ccTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTgaaga                                                                           >  1:2004539/1‑150 (MQ=255)
                                                                               tGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGTGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGTGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACt                                                                        <  1:4644218/150‑1 (MQ=255)
                                                                               tGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACt                                                                        <  1:1002944/150‑1 (MQ=255)
                                                                                gACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACtt                                                                       <  2:5595812/150‑1 (MQ=255)
                                                                                 aCGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTg                                                                      <  2:8330152/150‑1 (MQ=255)
                                                                                      aTCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAg                                                                 >  1:4856476/1‑150 (MQ=255)
                                                                                          tCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACCGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTca                                                             <  2:8464477/150‑1 (MQ=255)
                                                                                          tCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTca                                                             <  1:3236643/150‑1 (MQ=255)
                                                                                          tCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTca                                                             <  1:5664823/150‑1 (MQ=255)
                                                                                            gTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTcagc                                                           >  1:7619227/1‑150 (MQ=255)
                                                                                            gTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTcagc                                                           <  1:4122773/150‑1 (MQ=255)
                                                                                                 cccAGGCGTTGATCTGGCCACGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAgccagcg                                                      >  2:6918573/1‑150 (MQ=255)
                                                                                                  ccAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAgccagcgc                                                     <  1:1868839/150‑1 (MQ=255)
                                                                                                     ggCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAgccagcgccag                                                  <  1:4057997/150‑1 (MQ=255)
                                                                                                        gTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCAc                                               >  1:4315987/1‑150 (MQ=255)
                                                                                                        gTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCCGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCAc                                               <  2:7029676/150‑1 (MQ=255)
                                                                                                           gATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACgtt                                            <  1:5567330/150‑1 (MQ=255)
                                                                                                           gATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACgtt                                            <  1:7310755/150‑1 (MQ=255)
                                                                                                            aTCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACgttg                                           >  1:1574249/1‑150 (MQ=255)
                                                                                                               tGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTc                                        <  1:6962533/150‑1 (MQ=255)
                                                                                                               tGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTc                                        >  1:8168008/1‑150 (MQ=255)
                                                                                                                ggCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCt                                       <  2:554019/150‑1 (MQ=255)
                                                                                                                 gCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCtt                                      <  1:144116/150‑1 (MQ=255)
                                                                                                                   cAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGa                                    >  1:4509640/1‑150 (MQ=255)
                                                                                                                    aGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGaa                                   <  2:6113810/150‑1 (MQ=255)
                                                                                                                    aGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGaa                                   <  1:1590173/150‑1 (MQ=255)
                                                                                                                     gggcgggcAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAAc                                  >  2:4023633/1‑150 (MQ=255)
                                                                                                                       gcgggcAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACtt                                <  1:610450/150‑1 (MQ=255)
                                                                                                                                  aCGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGa                     >  1:8036808/1‑150 (MQ=255)
                                                                                                                                  aCGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGa                     <  2:1130253/150‑1 (MQ=255)
                                                                                                                                  aCGGCGCATCCGGGAACAGGGCTGTCAGTTCCGGCACCAGCGGGCCGTAGGGGCCGTGCTCGAAGCGGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACATGTCTGGCTCGa                     >  1:625317/1‑150 (MQ=17)
                                                                                                                                     gcgcATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGatat                  >  1:7845423/1‑150 (MQ=255)
                                                                                                                                      cgcATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATc                 <  2:8126750/150‑1 (MQ=255)
                                                                                                                                         aTCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCTCGTTGTTCTTGAACTTGTCTGGCTCGATATcgcg              <  2:4687433/150‑1 (MQ=255)
                                                                                                                                         aTCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATcgcg              >  1:2376643/1‑150 (MQ=255)
                                                                                                                                         aTCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATcgcg              >  2:1480014/1‑150 (MQ=255)
                                                                                                                                              ggAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCa         >  2:2914803/1‑150 (MQ=255)
                                                                                                                                                     ggCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGAcag  >  1:3823223/1‑150 (MQ=255)
                                                                                                                                                     ggCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGAcag  >  2:3174605/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGACAG  >  NC_002947/5992271‑5992569

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: ‑