Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,109,511 G→A 6.2% intergenic (+50/‑42) 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,5110GA 6.2% 1013.4 / 35.5 322intergenic (+50/‑42)murA/hisGUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase/ATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (166/136);  new base A (6/14);  total (172/150)
Fisher's exact test for biased strand distribution p-value = 3.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 333 total aligned reads displayed.

TACCCTGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTG  >  NC_002947/1109362‑1109660
                                                                                                                                                     |                                                                                                                                                     
  cccTGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGAt                                                                                                                                                     <  1:1512499/150‑1 (MQ=255)
   ccTGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATg                                                                                                                                                    <  2:6544440/150‑1 (MQ=255)
   ccTGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATg                                                                                                                                                    >  1:2172440/1‑150 (MQ=255)
   ccTGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTGTGAATTGAATGCGGCCAGGTCGATg                                                                                                                                                    >  1:1727363/1‑150 (MQ=255)
    cTGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGa                                                                                                                                                   <  1:1952604/150‑1 (MQ=255)
     tGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGAc                                                                                                                                                  >  1:4972668/1‑150 (MQ=255)
     tGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGATAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGAc                                                                                                                                                  <  1:5080752/150‑1 (MQ=255)
       atcgatcgCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCt                                                                                                                                                <  2:4524085/150‑1 (MQ=255)
        tcgatcgCATCTACCACCTCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCATCGACCTg                                                                                                                                               <  2:2689270/150‑1 (MQ=255)
          gatcgCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGc                                                                                                                                             <  2:4663429/150‑1 (MQ=255)
           atcgCATCTACCACCTCGCCCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTCCCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCATCGACCTGGcc                                                                                                                                            <  2:4293689/150‑1 (MQ=255)
             cgCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCgg                                                                                                                                          <  1:2789514/150‑1 (MQ=255)
               cATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgca                                                                                                                                        <  1:4657295/150‑1 (MQ=255)
               cATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCATCGACCTGGCCgcat                                                                                                                                        <  2:234137/150‑4 (MQ=255)
                aTCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgcag                                                                                                                                       >  1:4246905/1‑150 (MQ=255)
                aTCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgcag                                                                                                                                       >  2:318456/1‑150 (MQ=255)
                aTCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgcag                                                                                                                                       >  2:3357472/1‑150 (MQ=255)
                aTCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgcag                                                                                                                                       >  1:3168100/1‑150 (MQ=255)
                aTCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgcag                                                                                                                                       >  2:6483155/1‑150 (MQ=255)
                  cTACCCCATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCGCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGgcagca                                                                                                                                     <  2:176904/150‑1 (MQ=255)
                       acaTCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGc                                                                                                                                <  1:2705628/150‑1 (MQ=255)
                       acaTCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGc                                                                                                                                <  1:6170767/150‑1 (MQ=255)
                         aTCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCtt                                                                                                                              >  1:6349077/1‑150 (MQ=255)
                         aTCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCATCGACCTGGCCgcat                                                                                                                                        >  2:1116234/1‑137 (MQ=255)
                           cGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTaa                                                                                                                            <  2:4229393/150‑1 (MQ=255)
                            gACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAg                                                                                                                           >  1:5784166/1‑150 (MQ=255)
                            gACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAg                                                                                                                           <  2:1875964/150‑1 (MQ=255)
                             aCCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAgg                                                                                                                          >  1:2462503/1‑150 (MQ=255)
                             aCCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTGTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAgg                                                                                                                          >  2:5582371/1‑150 (MQ=255)
                              ccGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGa                                                                                                                         <  1:854353/150‑1 (MQ=255)
                                   gTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGAc                                                                                                                    >  1:2349914/1‑150 (MQ=255)
                                    ttACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACg                                                                                                                   >  1:665696/1‑150 (MQ=255)
                                          gTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCaa                                                                                                             >  1:5537441/1‑150 (MQ=255)
                                              aTCGAGTAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGtt                                                                                                         <  1:4734398/150‑1 (MQ=255)
                                              aTCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGtt                                                                                                         <  2:5722144/150‑1 (MQ=255)
                                              aTCGAGGAAAAACTGCAGAGGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAAGGCGGCCTGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGtt                                                                                                         <  1:2843784/150‑1 (MQ=255)
                                               tCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTg                                                                                                        >  1:2050976/1‑150 (MQ=255)
                                               tCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTg                                                                                                        >  1:1787782/1‑150 (MQ=255)
                                                cGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGa                                                                                                       >  2:2542563/1‑150 (MQ=255)
                                                 gAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGAc                                                                                                      >  1:4214235/1‑150 (MQ=255)
                                                 gAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGAc                                                                                                      >  1:3642867/1‑150 (MQ=255)
                                                  aGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGAcc                                                                                                     >  2:3699061/1‑150 (MQ=255)
                                                    gAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCAt                                                                                                   >  2:6188908/1‑150 (MQ=255)
                                                      aaaaCTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATcg                                                                                                 <  2:1918884/150‑1 (MQ=255)
                                                      aaaaCTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATcg                                                                                                 >  1:2561594/1‑150 (MQ=255)
                                                      aaaaCTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTGTTGAATTGATTGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATcg                                                                                                 >  1:5311716/1‑150 (MQ=255)
                                                            gCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCttt                                                                                           <  1:3801478/150‑1 (MQ=255)
                                                                    tCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGc                                                                                   >  2:6010365/1‑150 (MQ=255)
                                                                     cGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGcc                                                                                  <  2:3617585/150‑1 (MQ=255)
                                                                     cGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCAACTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGcc                                                                                  <  1:4997169/150‑1 (MQ=255)
                                                                        cgcgAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCa                                                                                                              <  2:4815230/119‑1 (MQ=255)
                                                                          cgAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTAtt                                                                             >  1:536843/1‑150 (MQ=255)
                                                                           gAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATtc                                                                            >  1:3598003/1‑150 (MQ=255)
                                                                           gAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATtc                                                                            >  1:373973/1‑150 (MQ=255)
                                                                               aTCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTcga                                                                        <  2:1277535/150‑1 (MQ=255)
                                                                                  cgccgcGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTcgacga                                                                     <  2:6220158/150‑1 (MQ=255)
                                                                                   gccgcGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGAt                                                                    >  1:1327778/1‑150 (MQ=255)
                                                                                          aCCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGc                                                             <  2:5811455/150‑1 (MQ=255)
                                                                                           ccGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGcc                                                            >  2:4456501/1‑150 (MQ=255)
                                                                                            cGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCg                                                           >  2:3953974/1‑150 (MQ=255)
                                                                                             gggCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGt                                                          >  2:3799876/1‑150 (MQ=255)
                                                                                                    tACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGc                                                   >  1:3447417/1‑150 (MQ=255)
                                                                                                       cccGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGa                                                >  2:4060331/1‑150 (MQ=255)
                                                                                                       cccGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGa                                                >  2:4110379/1‑150 (MQ=255)
                                                                                                       cccGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGa                                                >  2:1599861/1‑150 (MQ=255)
                                                                                                         cgcgGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAgg                                              >  1:5928983/1‑150 (MQ=255)
                                                                                                              gTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGc                                         >  2:436601/1‑150 (MQ=255)
                                                                                                                 cAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATc                                      <  1:3195160/150‑1 (MQ=255)
                                                                                                                  aaGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCg                                     >  1:2655109/1‑150 (MQ=255)
                                                                                                                  aaGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCg                                     >  2:6175550/1‑150 (MQ=255)
                                                                                                                   cggACGCACCTGGTGCCTTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGTCCGACGTTCCAAGGTTGACCATCGCGCTGTCCAAAGGCCGTATTATCGACGAGACCTTGCCGTTGCTCGCCGTGGCCGGCATCGt                                    <  2:5697584/149‑1 (MQ=255)
                                                                                                                     gACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCGGGCCGGCAGCAGCCGCTTAAGGTCCGTCGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGc                                  <  1:6754497/150‑1 (MQ=255)
                                                                                                                        gccccTTGTGAATTGAAGGCGGCCAGGTCGATGACCGGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccga                               <  1:1780632/147‑1 (MQ=255)
                                                                                                                          aCCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccgacc                             >  2:4859637/1‑150 (MQ=255)
                                                                                                                          aCCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccgacc                             >  1:30743/1‑150 (MQ=255)
                                                                                                                          aCCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccgacc                             <  2:1417957/150‑1 (MQ=255)
                                                                                                                           ccTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccgaccg                            <  2:5161091/150‑1 (MQ=255)
                                                                                                                            cTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGccgaccga                           >  1:1591176/1‑150 (MQ=255)
                                                                                                                             ttttGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGaa                          >  1:4454394/1‑150 (MQ=255)
                                                                                                                                tgaattgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAAc                       <  2:2917417/150‑1 (MQ=255)
                                                                                                                                 gaattgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAAcc                      >  1:1034705/1‑150 (MQ=255)
                                                                                                                                    ttgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCgg                   >  2:609351/1‑150 (MQ=255)
                                                                                                                                    ttgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCgg                   >  1:4011885/1‑150 (MQ=255)
                                                                                                                                    ttgaatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCgg                   <  1:3831503/150‑1 (MQ=255)
                                                                                                                                       aatGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACa                >  2:212150/1‑150 (MQ=255)
                                                                                                                                        tgccGGCCAGGTCATCGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACaa               >  2:3449714/4‑150 (MQ=255)
                                                                                                                                          gCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAaga             <  2:4679286/150‑1 (MQ=255)
                                                                                                                                          ccGGCCAGGTCATCGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAaga             >  2:4960670/2‑150 (MQ=255)
                                                                                                                                            ggCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGc           >  2:2187447/1‑150 (MQ=255)
                                                                                                                                            ggCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGc           >  1:6482812/1‑150 (MQ=255)
                                                                                                                                            ggCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGc           >  2:5943481/1‑150 (MQ=255)
                                                                                                                                             gCCAGGTCATCGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGcc          >  2:1317851/1‑150 (MQ=255)
                                                                                                                                                aGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCa       >  1:89186/1‑150 (MQ=255)
                                                                                                                                                  gTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCaaa     <  1:1070215/150‑1 (MQ=255)
                                                                                                                                                    cGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACt   <  1:4839396/150‑1 (MQ=255)
                                                                                                                                                     gATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTg  >  2:1717940/1‑150 (MQ=255)
                                                                                                                                                     gATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTg  >  1:5229933/1‑150 (MQ=255)
                                                                                                                                                     gATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTg  >  1:1183214/1‑150 (MQ=255)
                                                                                                                                                     gATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTg  >  2:6858198/1‑150 (MQ=255)
                                                                                                                                                     gATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCCCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCACTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTg  >  1:6155282/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACCCTGATCGATCGCATCTACCACATCGACCGTGGTTACGAGTGCATCGAGGAAAAACTGCAGATGCTCGGCGCGAAAATCCGCCGCGTACCGGGCTAGTACCCCGCGGGTGCAAGGACGCACCTTTTGAATTGAATGCGGCCAGGTCGATGACCTGGCCGGCAGCAGCCGCTTAAGGACCGACGTTCCAATGTTGACCATCGCGCTTTCCAAAGGCCGTATTCTCGACGATACCTTGCCGTTGCTCGCCGAGGCCGGCATCGTGCCGACCGAAAACCCGGACAAGAGCCGCAAACTG  >  NC_002947/1109362‑1109660

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: