Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,109,513 T→C 5.9% intergenic (+52/‑40) murA → / → hisG UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase/ATP phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,109,5130TC 5.9% 1088.7 / 34.2 340intergenic (+52/‑40)murA/hisGUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase/ATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (177/143);  new base C (6/14);  total (183/157)
Fisher's exact test for biased strand distribution p-value = 3.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

Only 100 of 347 total aligned reads displayed.

CCCTGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTGAT  >  NC_002947/1109364‑1109662
                                                                                                                                                     |                                                                                                                                                     
 ccTGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATg                                                                                                                                                      >  1:2177684/1‑150 (MQ=255)
 ccTGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATg                                                                                                                                                      <  2:6544283/150‑1 (MQ=255)
 ccTGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTGTGAATTGAATGCGGCCAGGTCGATg                                                                                                                                                      >  1:1727363/1‑150 (MQ=255)
   tGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTCCCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGAc                                                                                                                                                    >  1:1151076/1‑150 (MQ=255)
   tGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGTAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGAc                                                                                                                                                    >  2:3458209/1‑150 (MQ=255)
   tGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGAc                                                                                                                                                    >  1:6843116/1‑150 (MQ=255)
   tGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGAc                                                                                                                                                    >  1:4972668/1‑150 (MQ=255)
   tGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTGTTGAATTGAATGCGGCCAGGTCGATGAc                                                                                                                                                    >  1:2481601/1‑150 (MQ=255)
   tGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGATAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGAc                                                                                                                                                    <  1:5080752/150‑1 (MQ=255)
    gatcgatcGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGAcc                                                                                                                                                   >  2:2260984/1‑150 (MQ=255)
         atcgCATCTACCACCTCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCATCGACCTGGcc                                                                                                                                              <  2:4473093/150‑1 (MQ=255)
         atcgCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCATCGACCTGGcc                                                                                                                                              <  1:524642/150‑1 (MQ=255)
         atcgCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCATCGACCTGGcc                                                                                                                                              <  1:2724410/150‑1 (MQ=255)
          tcgCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCg                                                                                                                                             >  1:5511209/1‑150 (MQ=255)
          tcgCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCg                                                                                                                                             >  2:1429960/1‑150 (MQ=255)
           cgCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCgg                                                                                                                                            <  2:340997/150‑1 (MQ=255)
           cgCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCgg                                                                                                                                            <  1:6254120/150‑1 (MQ=255)
            gCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgc                                                                                                                                           >  1:307659/1‑150 (MQ=255)
             cATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgca                                                                                                                                          <  1:6451676/150‑1 (MQ=255)
             cATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgca                                                                                                                                          <  1:4657295/150‑1 (MQ=255)
              aTCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgcag                                                                                                                                         >  1:4246905/1‑150 (MQ=255)
                cTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgcagca                                                                                                                                       >  1:6677447/1‑150 (MQ=255)
                cTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgcagca                                                                                                                                       <  1:6525712/150‑1 (MQ=255)
                 tACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGcagcag                                                                                                                                      >  1:2010275/1‑150 (MQ=255)
                      caTCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCt                                                                                                                                 >  2:1625025/1‑150 (MQ=255)
                       aTCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCATCGACCTGGCCgcat                                                                                                                                          >  2:1116234/1‑137 (MQ=255)
                        tCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTa                                                                                                                               >  1:5899996/1‑150 (MQ=255)
                         cGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTaa                                                                                                                              <  2:4229393/150‑1 (MQ=255)
                           aCCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTGTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAgg                                                                                                                            >  2:5582371/1‑150 (MQ=255)
                             cGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGAc                                                                                                                          >  1:3325262/1‑150 (MQ=255)
                              gTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGAcc                                                                                                                         >  2:5919652/1‑150 (MQ=255)
                                    aCGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTCCCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTGTGGAAGGGAATGCGGCCAGGTCGCTGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGtt                                                                                                                   >  2:3282334/1‑150 (MQ=255)
                                             tCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTg                                                                                                          >  1:1787782/1‑150 (MQ=255)
                                             tCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTg                                                                                                          >  1:2050976/1‑150 (MQ=255)
                                               gAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAATGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGAc                                                                                                        <  2:1755246/150‑1 (MQ=255)
                                               gAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGAc                                                                                                        <  2:4638361/150‑1 (MQ=255)
                                               gAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGAc                                                                                                        >  1:4214235/1‑150 (MQ=255)
                                                aGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGAcc                                                                                                       >  2:3699061/1‑150 (MQ=255)
                                                   aaaaaCTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATc                                                                                                    <  2:3621051/150‑1 (MQ=255)
                                                        cTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCt                                                                                               >  2:6178399/1‑150 (MQ=255)
                                                        cTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCt                                                                                               >  1:5312338/1‑150 (MQ=255)
                                                         tGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCtt                                                                                              >  2:4290344/1‑150 (MQ=255)
                                                           cAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTc                                                                                            <  2:6205937/150‑1 (MQ=255)
                                                               tGCTCGGCGAGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCTAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCGGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCaaa                                                                                        <  1:5022104/150‑1 (MQ=255)
                                                                gCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAg                                                                                       <  1:6170325/150‑1 (MQ=255)
                                                                       gcgAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTAt                                                                                >  2:2570402/1‑150 (MQ=255)
                                                                                 gccgcGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTCTTCTCGACGAt                                                                      >  1:3773866/1‑150 (MQ=255)
                                                                                 gccgcGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGAt                                                                      >  1:1327778/1‑150 (MQ=255)
                                                                                 gccgcGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGAt                                                                      <  1:3987520/150‑1 (MQ=255)
                                                                                          cGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCg                                                             >  2:3952380/1‑150 (MQ=255)
                                                                                           gggCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGt                                                            >  2:3799876/1‑150 (MQ=255)
                                                                                             gCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTg                                                          <  2:1914731/150‑1 (MQ=255)
                                                                                              cTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGc                                                         >  1:1499848/1‑150 (MQ=255)
                                                                                                     cccGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGa                                                  >  2:4110379/1‑150 (MQ=255)
                                                                                                     cccGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGa                                                  >  2:4060331/1‑150 (MQ=255)
                                                                                                     cccGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGa                                                  >  1:3030301/1‑150 (MQ=255)
                                                                                                     cccGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGa                                                  <  1:608204/150‑1 (MQ=255)
                                                                                                     cccGCGGGTGCAAGGACGCACCTTTTGAATTGAAAGCGGCCAGGTCGATGACCAGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAAAGTTGACCAACGCGCTTTCCAAAGGCCGTATTCTCGACGAAACATTGCCGTAGCTCGCCGa                                                  >  2:5105890/1‑150 (MQ=255)
                                                                                                            gTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGc                                           >  2:436601/1‑150 (MQ=255)
                                                                                                                  ggACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTg                                     <  2:1280528/150‑1 (MQ=255)
                                                                                                                   gACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGc                                    <  2:6351612/150‑1 (MQ=255)
                                                                                                                      gccccTTGTGAATTGAAGGCGGCCAGGTCGATGACCGGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccga                                 <  1:1780632/147‑1 (MQ=255)
                                                                                                                      gCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccga                                 >  2:4674457/1‑150 (MQ=255)
                                                                                                                      gCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccga                                 <  1:2574039/150‑1 (MQ=255)
                                                                                                                        aCCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccgacc                               >  2:4859637/1‑150 (MQ=255)
                                                                                                                        aCCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTCGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccgacc                               >  1:888024/1‑150 (MQ=255)
                                                                                                                         ccTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccgaccg                              >  2:1356252/1‑150 (MQ=255)
                                                                                                                         ccTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccgaccg                              <  2:5701423/150‑1 (MQ=255)
                                                                                                                          cTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccgaccga                             >  1:1591176/1‑150 (MQ=255)
                                                                                                                           ttttGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGaa                            >  2:1803347/1‑150 (MQ=255)
                                                                                                                           ttttGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGaa                            >  1:280059/1‑150 (MQ=255)
                                                                                                                           ttttGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGaa                            >  1:4454394/1‑150 (MQ=255)
                                                                                                                               gaattgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAAcc                        >  1:1034705/1‑150 (MQ=255)
                                                                                                                                atttgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAAccc                       <  1:3187194/148‑1 (MQ=255)
                                                                                                                                aattgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAAccc                       <  2:2172440/150‑1 (MQ=255)
                                                                                                                                aattgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAAccc                       <  2:5784166/150‑1 (MQ=255)
                                                                                                                                 attgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCg                      <  2:3030301/150‑1 (MQ=255)
                                                                                                                                 attgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCg                      <  2:6107864/150‑1 (MQ=255)
                                                                                                                                  ttgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCgg                     >  2:4790300/1‑150 (MQ=255)
                                                                                                                                  ttgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCgg                     <  1:3831503/150‑1 (MQ=255)
                                                                                                                                  ttgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCgg                     >  1:3727958/1‑150 (MQ=255)
                                                                                                                                        gCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAaga               <  1:4734657/150‑1 (MQ=255)
                                                                                                                                         cGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAagag              >  1:2662623/1‑150 (MQ=255)
                                                                                                                                          ggCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGc             <  2:2010275/150‑1 (MQ=255)
                                                                                                                                          ggCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGc             >  2:5943750/1‑150 (MQ=255)
                                                                                                                                          ggCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGc             >  1:3641590/1‑150 (MQ=255)
                                                                                                                                          ggCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCGTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGGTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGc             <  1:5582371/150‑1 (MQ=255)
                                                                                                                                            ccAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCg           >  2:6056281/1‑150 (MQ=255)
                                                                                                                                               ggTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCaa        >  1:998903/1‑150 (MQ=255)
                                                                                                                                                gTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACCAGAGCCGCaaa       >  1:3576236/1‑150 (MQ=255)
                                                                                                                                                gTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCaaa       <  1:1056072/150‑1 (MQ=255)
                                                                                                                                                  cGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACt     <  2:4933549/150‑1 (MQ=255)
                                                                                                                                                  cGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACt     <  1:4839396/150‑1 (MQ=255)
                                                                                                                                                   gATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTg    >  1:2420242/1‑150 (MQ=255)
                                                                                                                                                   gATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTg    >  2:6858198/1‑150 (MQ=255)
                                                                                                                                                    aTGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTGa   <  2:4395494/150‑1 (MQ=255)
                                                                                                                                                    aTGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTGa   <  2:2763881/150‑1 (MQ=255)
                                                                                                                                                    aTGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTGa   <  2:6664846/150‑1 (MQ=255)
                                                                                                                                                     tGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTGat  >  2:5667739/1‑150 (MQ=255)
                                                                                                                                                     tGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTGat  <  2:2462503/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCTGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTGAT  >  NC_002947/1109364‑1109662

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: