Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 566,623 C→T 5.5% intergenic (+32/+41) bfr‑I → / ← uvrA bacterioferritin/uvrABC excinuclease protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947566,6230CT 5.5% 816.9 / 20.4 254intergenic (+32/+41)bfr‑I/uvrAbacterioferritin/uvrABC excinuclease protein A
Reads supporting (aligned to +/- strand):  ref base C (124/116);  new base T (6/8);  total (130/124)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

TGACATCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTT  >  NC_002947/566474‑566771
                                                                                                                                                     |                                                                                                                                                    
 gACATCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCa                                                                                                                                                     <  1:3440537/150‑1 (MQ=255)
     tCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCAt                                                                                                                                                 >  2:1226944/1‑150 (MQ=255)
     tCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCAt                                                                                                                                                 <  2:5238002/150‑1 (MQ=255)
        tgcgtgcgCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCATATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCgg                                                                                                                                              <  1:2890894/150‑1 (MQ=255)
         gcgtgcgCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGt                                                                                                                                             <  1:5341029/150‑1 (MQ=255)
             gcgcAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCa                                                                                                                                         <  2:2168140/150‑1 (MQ=255)
               gcAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAgg                                                                                                                                       >  1:2270179/1‑150 (MQ=255)
                  gCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGc                                                                                                                                    <  1:1765652/150‑1 (MQ=255)
                  gCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCAAAGACTCAAGTAAAAGCCCCTGGTACCGATGGTGCCAGGGGc                                                                                                                                    <  1:68770/150‑1 (MQ=255)
                     ggCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCttt                                                                                                                                 >  2:2931224/1‑150 (MQ=255)
                      gCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCtttt                                                                                                                                >  2:890775/1‑150 (MQ=255)
                        cGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTATTTc                                                                                                                              <  2:358196/150‑1 (MQ=255)
                         gACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCACAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCa                                                                                                                             >  1:1281071/1‑150 (MQ=255)
                          acacCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAt                                                                                                                            <  2:759754/150‑1 (MQ=255)
                           cacCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAtt                                                                                                                           <  1:3632864/150‑1 (MQ=255)
                           cacCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAtt                                                                                                                           <  2:3679967/150‑1 (MQ=255)
                            acCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTg                                                                                                                          >  2:5000124/1‑150 (MQ=255)
                            acCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTg                                                                                                                          >  1:4647545/1‑150 (MQ=255)
                            acCGAAGAAGAACACACCTACTCGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTg                                                                                                                          <  2:3203669/150‑1 (MQ=255)
                              cGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgcg                                                                                                                        >  2:53628/1‑150 (MQ=255)
                               gaagaaGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgcgc                                                                                                                       <  2:4837859/150‑1 (MQ=255)
                                aagtagATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAATTGCCCCTGGCACCATCGGTAGCAGGGGGTTTTTCTTTgcgcg                                                                                                                      >  2:4261450/1‑149 (MQ=255)
                                aagaagATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGcc                                                                                                                      >  1:810658/1‑150 (MQ=255)
                                 agaagaTCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCg                                                                                                                     <  2:1731219/150‑1 (MQ=255)
                                  gaagaTCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCgg                                                                                                                    <  2:5510815/150‑1 (MQ=255)
                                    agaTCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCa                                                                                                                  >  1:5295132/1‑150 (MQ=255)
                                     gaTCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCaa                                                                                                                 <  1:4137398/150‑1 (MQ=255)
                                       tCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATa                                                                                                               >  1:5070748/1‑150 (MQ=255)
                                       tCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATa                                                                                                               >  2:329846/1‑150 (MQ=255)
                                        cacaCCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATAc                                                                                                              <  1:4489900/150‑1 (MQ=255)
                                           acCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCa                                                                                                           >  1:1570330/1‑150 (MQ=255)
                                              tACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGc                                                                                                        <  2:1363123/150‑1 (MQ=255)
                                               aCTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGcc                                                                                                       >  1:371184/1‑150 (MQ=255)
                                               aCTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGcc                                                                                                       >  1:1627827/1‑150 (MQ=255)
                                               aCTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGACAATACTCAGGcc                                                                                                       <  2:4604622/150‑1 (MQ=255)
                                                 tGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCt                                                                                                     >  1:3321298/1‑150 (MQ=255)
                                                   gCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTAt                                                                                                   >  1:2641178/1‑150 (MQ=255)
                                                     ttGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATcc                                                                                                 >  2:2995824/1‑150 (MQ=255)
                                                       gagaAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCg                                                                                               >  1:1875413/1‑150 (MQ=255)
                                                        agaAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGc                                                                                              >  1:872524/1‑150 (MQ=255)
                                                           agcagcAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGTCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTcc                                                                                           <  1:1231128/150‑1 (MQ=255)
                                                           agcagcAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTcc                                                                                           >  1:4845254/1‑150 (MQ=255)
                                                             cagcagGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAg                                                                                         <  1:3427704/150‑1 (MQ=255)
                                                               gcagGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCa                                                                                       >  2:2061302/1‑150 (MQ=255)
                                                                   ggCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATgg                                                                                   >  1:4293061/1‑150 (MQ=255)
                                                                       tGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTc                                                                               >  1:3631838/1‑150 (MQ=255)
                                                                       tGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTc                                                                               <  1:455995/150‑1 (MQ=255)
                                                                        gATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCa                                                                              <  2:363276/150‑1 (MQ=255)
                                                                         aTCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCACTGGTTTCAg                                                                             >  2:2710120/1‑150 (MQ=255)
                                                                         aTCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAg                                                                             <  2:5142659/150‑1 (MQ=255)
                                                                         aTCAAGGCCATCGCGCTGGAAAACTACCAGCAGTCGCAGATGAAAGTAATCCACAGACACAACAAAAAGCCCCTCGCACCATCGGTACCAGGGGCAAATTCATTGCGCCGGCAATACTCAGGCCCTATCCCCCTCCAGCAATGGTTTCAg                                                                             <  2:4775596/150‑1 (MQ=255)
                                                                          tCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTCCCAGGGGCTTTGTCATTGCGCCGGCAATACTCAGGCCCTCTCCCGCGCCAGCAATGGTTTCAgg                                                                            >  1:3888648/1‑150 (MQ=255)
                                                                          tCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAgg                                                                            >  2:4981297/1‑150 (MQ=255)
                                                                             aGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCCATGGTTTCAGGTaa                                                                         >  2:2869424/1‑150 (MQ=255)
                                                                             aGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTaa                                                                         <  2:2661507/150‑1 (MQ=255)
                                                                              ggCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAAt                                                                        <  2:4344933/150‑1 (MQ=255)
                                                                               gCCATCGGGCTGGAAATCTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAATAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATg                                                                       >  2:3843665/1‑150 (MQ=255)
                                                                               gCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATg                                                                       >  2:3740625/1‑150 (MQ=255)
                                                                                  aTCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTc                                                                               >  1:4891447/1‑139 (MQ=255)
                                                                                   tCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCg                                                                   <  1:3098177/150‑1 (MQ=255)
                                                                                   tCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCg                                                                   <  2:30302/150‑1 (MQ=255)
                                                                                    cGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCgg                                                                  >  1:5543930/1‑150 (MQ=255)
                                                                                     gggCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGt                                                                 <  2:1141075/150‑1 (MQ=255)
                                                                                      ggCTGGAAAACTACCTGCATTCGCAGATGTAAGTTATCAACAGACTAAAGTAAAAGCCCCTGGCACCATCTGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGta                                                                <  1:4204289/150‑1 (MQ=255)
                                                                                         tGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTatat                                                             <  1:2177898/150‑1 (MQ=255)
                                                                                          ggAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGcc                                                                    <  1:1077041/142‑1 (MQ=255)
                                                                                          ggAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGcc                                                                    >  2:1077041/1‑142 (MQ=255)
                                                                                           gAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGa                                                           >  2:2681811/1‑150 (MQ=255)
                                                                                           gAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGa                                                           <  1:566749/150‑1 (MQ=255)
                                                                                             aaaCTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACt                                                         <  1:329846/150‑1 (MQ=255)
                                                                                               aCTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGc                                                       <  2:151039/150‑1 (MQ=255)
                                                                                                 tACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCtt                                                     >  1:4616458/1‑150 (MQ=255)
                                                                                                   ccTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCa                                                   >  1:1769916/1‑150 (MQ=255)
                                                                                                      gCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCAtct                                                >  2:4710450/1‑150 (MQ=255)
                                                                                                          tCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCgg                                                                  <  1:2672468/128‑1 (MQ=255)
                                                                                                          tCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCt                                            >  1:5211964/1‑150 (MQ=255)
                                                                                                            gCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCa                                          >  1:3080394/1‑150 (MQ=255)
                                                                                                            gCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCa                                          >  2:539211/1‑150 (MQ=255)
                                                                                                                    aaGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCtcttct                                  <  2:5124618/150‑1 (MQ=255)
                                                                                                                    aaGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAGTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCtcttct                                  <  2:5525029/150‑1 (MQ=255)
                                                                                                                     aGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTg                                 <  2:4467260/150‑1 (MQ=255)
                                                                                                                     aGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTg                                 <  2:4119771/150‑1 (MQ=255)
                                                                                                                      gTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTgg                                >  2:2991616/1‑150 (MQ=255)
                                                                                                                       ttATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTATTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGgt                               <  1:3241295/150‑1 (MQ=255)
                                                                                                                            cacaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTAcc                          >  2:687124/1‑150 (MQ=255)
                                                                                                                             acaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCCATGGTTTCAGGTCATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCg                         >  2:2413101/1‑150 (MQ=255)
                                                                                                                              caGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGc                        <  2:2489981/150‑1 (MQ=255)
                                                                                                                               cgACTCACGTAAAAGCCCCGGGCACCTTCGGTACCAGGGGCTTTTTCAGTGCGCCGGCAATTCTCAGGCCCGTTGGCGCTCCCGCTATGGTTTGAGGTAATGGCCGGTATATGTCTGCTTCATCTAGCTCCGCTCTTCTGGTGTACCGct                       <  2:2195424/149‑2 (MQ=255)
                                                                                                                               aGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCa                       <  2:1854942/150‑1 (MQ=255)
                                                                                                                                   tCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCa                   >  1:5415542/1‑150 (MQ=255)
                                                                                                                                        tAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCAGCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGAt              >  1:933230/1‑150 (MQ=255)
                                                                                                                                        aaaaaGCCCCTGGTACCGATGGTGCCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCtctt                                    <  2:1851024/127‑1 (MQ=38)
                                                                                                                                        aaaaaGCCCCTGGTACCGATGGTGCCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCtctt                                    >  1:1851024/2‑128 (MQ=38)
                                                                                                                                         aaaaGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATc             >  1:1086089/1‑150 (MQ=255)
                                                                                                                                          aaaGCCCCGGGTACCGATGGTGCCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCAAGCAATGGTTTAAGGTAATGGCCGGTATATGACTGCTACATCTCGCTCAGCTCTTCTGTTGTACCGCAGGCAATGAtat            >  1:196458/1‑148 (MQ=255)
                                                                                                                                               cccTGGTACCGATGGTGCCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGTCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCa       >  1:4018845/1‑150 (MQ=255)
                                                                                                                                               cccTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCa       >  1:603203/1‑150 (MQ=255)
                                                                                                                                                ccTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAc      >  1:5158069/1‑150 (MQ=255)
                                                                                                                                                    gCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCtt  >  1:2537313/1‑150 (MQ=255)
                                                                                                                                                    gCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCtt  <  1:3740625/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGACATCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTT  >  NC_002947/566474‑566771

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: