Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 566,627 A→G 6.9% intergenic (+36/+37) 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,6270AG 6.9% 761.1 / 29.7 246intergenic (+36/+37)bfr‑I/uvrAbacterioferritin/uvrABC excinuclease protein A
Reads supporting (aligned to +/- strand):  ref base A (120/109);  new base G (6/11);  total (126/120)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

Only 100 of 251 total aligned reads displayed.

TCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTTGGAAC  >  NC_002947/566479‑566776
                                                                                                                                                    |                                                                                                                                                     
tCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCAt                                                                                                                                                      >  2:5698423/1‑150 (MQ=255)
  cTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCAGCg                                                                                                                                                    >  2:3300319/1‑150 (MQ=255)
       tgcgCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTAcc                                                                                                                                               >  2:5172179/1‑150 (MQ=255)
       tgcgCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTAcc                                                                                                                                               >  1:3236253/1‑150 (MQ=255)
          gcAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAgg                                                                                                                                            <  2:4406319/150‑1 (MQ=255)
          gcAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAgg                                                                                                                                            >  1:9144823/1‑150 (MQ=255)
          gcAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAgg                                                                                                                                            <  2:2636981/150‑1 (MQ=255)
            aGCTGGCCGACACCGAAGAAGATCACACCGACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAgggg                                                                                                                                          <  2:5819299/150‑1 (MQ=255)
             gCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGTACCGATGGTGCCAGGGGc                                                                                                                                         <  2:6418603/150‑1 (MQ=255)
               tGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCtt                                                                                                                                       >  2:6271067/1‑150 (MQ=255)
                  ccGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGTACCGATGGTGCCAGGGGCtttta                                                                                                                                    <  1:2996163/150‑2 (MQ=39)
                    gACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCa                                                                                                                                  <  1:1193739/150‑1 (MQ=255)
                     acacCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAt                                                                                                                                 <  1:1047358/150‑1 (MQ=255)
                      cacCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAtt                                                                                                                                <  1:8145008/150‑1 (MQ=255)
                       acCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTg                                                                                                                               <  2:7070059/150‑1 (MQ=255)
                        ccGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgc                                                                                                                              <  2:2257277/150‑1 (MQ=255)
                         cGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgcg                                                                                                                             >  1:7696971/1‑150 (MQ=255)
                          gaagaaGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgcgc                                                                                                                            >  2:1631345/1‑150 (MQ=255)
                          gaagaaGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgcgc                                                                                                                            <  2:3125076/150‑1 (MQ=255)
                          gaagaaGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgcgc                                                                                                                            <  2:8273790/150‑1 (MQ=255)
                            agaagaTCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCg                                                                                                                          >  2:4313519/1‑150 (MQ=255)
                                      acCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCa                                                                                                                >  2:4812066/1‑150 (MQ=255)
                                      acCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCa                                                                                                                >  1:7574496/1‑150 (MQ=255)
                                      acCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCa                                                                                                                >  2:7603505/1‑150 (MQ=255)
                                      acCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCAACGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCa                                                                                                                >  2:6806104/1‑150 (MQ=255)
                                       ccTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAg                                                                                                               <  1:8481805/150‑1 (MQ=255)
                                         tACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGc                                                                                                             >  1:2624049/1‑150 (MQ=255)
                                            tGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCATATTTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCt                                                                                                          <  1:6556192/150‑1 (MQ=255)
                                            tGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCt                                                                                                          >  2:3600189/1‑150 (MQ=255)
                                                ttGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATcc                                                                                                      >  2:1720986/1‑150 (MQ=255)
                                                          gcagGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCAACAGACTCAAGTAAAAGCCACTGGCACCATCGGTACCAGGGGCTTTTTCATTGAGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCa                                                                                            >  1:4806582/1‑150 (MQ=255)
                                                            agGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAAt                                                                                          <  1:2763366/150‑1 (MQ=255)
                                                             gggCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATg                                                                                         <  2:6575656/150‑1 (MQ=255)
                                                              ggCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATgg                                                                                        >  1:4613908/1‑150 (MQ=255)
                                                              ggCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATgg                                                                                        >  1:1406322/1‑150 (MQ=255)
                                                                ccTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGtt                                                                                      >  1:8270173/1‑150 (MQ=255)
                                                                ccTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGtt                                                                                      >  1:981378/1‑150 (MQ=255)
                                                                  tGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTc                                                                                    <  1:2356719/150‑1 (MQ=255)
                                                                  tGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTc                                                                                    <  1:3302920/150‑1 (MQ=255)
                                                                  tGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTc                                                                                    <  2:4186619/150‑1 (MQ=255)
                                                                    aTCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAg                                                                                  <  1:8764969/150‑1 (MQ=255)
                                                                    aTCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAg                                                                                  <  1:86479/150‑1 (MQ=255)
                                                                    aTCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAg                                                                                  <  2:1058825/150‑1 (MQ=255)
                                                                    aTCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAg                                                                                  <  1:1328967/150‑1 (MQ=255)
                                                                       aaGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTa                                                                               >  2:1360366/1‑150 (MQ=255)
                                                                       aaGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTa                                                                               >  1:1712626/1‑150 (MQ=255)
                                                                            cATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGc                                                                          >  2:389680/1‑150 (MQ=255)
                                                                             aTCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGcc                                                                         >  1:4120469/1‑150 (MQ=255)
                                                                              tCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCg                                                                        <  2:4656122/150‑1 (MQ=255)
                                                                              tCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCg                                                                        <  1:3482537/150‑1 (MQ=255)
                                                                                gggCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGt                                                                      <  2:5944503/150‑1 (MQ=255)
                                                                                 ggCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGta                                                                     <  1:8893202/150‑1 (MQ=255)
                                                                                    tGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTatat                                                                  <  2:8822315/150‑1 (MQ=255)
                                                                                    tGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTatat                                                                  >  2:5333861/1‑150 (MQ=255)
                                                                                    tGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTatat                                                                  <  1:8209095/150‑1 (MQ=255)
                                                                                    tGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTatat                                                                  <  1:7880678/150‑1 (MQ=255)
                                                                                      gataaCTACCTGGAGTCGCTGATGTAAGTTATCCACATACTCAAGTAACAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGTCCCTATCCCGCTCCATCAATTGTTTCAGGTAATGGCCGGTATATGa                                                                <  2:1418771/147‑1 (MQ=255)
                                                                                                   aGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCg                                                   <  1:2472031/150‑1 (MQ=255)
                                                                                                    gTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGATTTTTAATTGCGCCTGCAATACTCGGGCCCTATCCCGCTCCAGCAATGGTTTCTGGTAATGGCCGGTATATGACTGCTTCATCTCGc                                                  >  1:7506310/1‑150 (MQ=255)
                                                                                                       gCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATg                                                                            >  1:2069089/1‑121 (MQ=255)
                                                                                                       gCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCa                                               >  1:4002019/1‑150 (MQ=255)
                                                                                                        cAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAg                                              <  2:9136876/150‑1 (MQ=255)
                                                                                                        cAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAg                                              <  2:1921835/150‑1 (MQ=255)
                                                                                                        cAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAg                                              <  1:7057829/150‑1 (MQ=255)
                                                                                                              tAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCtcttc                                        >  2:2684899/1‑150 (MQ=255)
                                                                                                                aGTTATCCACAGACTCAAGTAACAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTg                                      <  2:8669433/150‑1 (MQ=255)
                                                                                                                aGTTATCCACAGACTCAAGCAAAAGCCCCTGGCACCATCGGTATCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTg                                      <  2:8994770/150‑1 (MQ=255)
                                                                                                                   gaTCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCCGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGgtg                                   <  2:3693649/149‑1 (MQ=255)
                                                                                                                     tCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTa                                 <  2:7983803/150‑1 (MQ=255)
                                                                                                                     tCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCGCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTGCATCTCGCTCAGCTCTTCTGGTGTa                                 >  1:6716675/1‑150 (MQ=255)
                                                                                                                        acaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCg                              >  2:8589258/1‑150 (MQ=255)
                                                                                                                            aCTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAgg                          <  1:7355967/150‑1 (MQ=255)
                                                                                                                              tCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCACTGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCa                        <  1:1742182/150‑1 (MQ=255)
                                                                                                                               cAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCaa                       <  2:1238466/150‑1 (MQ=255)
                                                                                                                                 aGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATg                     >  1:565415/1‑150 (MQ=255)
                                                                                                                                 aGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATg                     <  1:4466653/150‑1 (MQ=255)
                                                                                                                                     aaaGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCt                 <  2:7256334/150‑1 (MQ=255)
                                                                                                                                      aaGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTg                <  1:5698423/150‑1 (MQ=255)
                                                                                                                                       aGCCCCTGGTACCGATGGTGCCAGGGGCTTTTACATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCAtct                                                     >  1:4856730/1‑112 (MQ=38)
                                                                                                                                       aGCCCCTGGTACCGATGGTGCCAGGGGCTTTTACATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCAtct                                                     <  2:4856730/112‑1 (MQ=38)
                                                                                                                                       aGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGc               <  1:7201963/150‑1 (MQ=255)
                                                                                                                                        gCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGcc              >  1:3112392/1‑150 (MQ=255)
                                                                                                                                        gCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGcc              <  2:7696971/150‑1 (MQ=255)
                                                                                                                                         ccccTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGccc             <  1:6327152/150‑1 (MQ=255)
                                                                                                                                          cccTGGTACCGATGGTGCCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCa            >  1:6828041/1‑150 (MQ=255)
                                                                                                                                          cccTGGTACCGATGGTGCCAGGGGCTTTTACATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCa            >  2:8062766/1‑150 (MQ=255)
                                                                                                                                          cccTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCa            >  1:7245330/1‑150 (MQ=255)
                                                                                                                                          cccTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCa            >  2:8377764/1‑150 (MQ=255)
                                                                                                                                           ccTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAg                                              <  1:5911665/115‑1 (MQ=255)
                                                                                                                                           ccTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAc           >  1:603437/1‑150 (MQ=255)
                                                                                                                                            cTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTACTGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAcc          >  2:8722263/1‑150 (MQ=255)
                                                                                                                                            cTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAcc          >  2:7774665/1‑150 (MQ=255)
                                                                                                                                            cTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAcc          >  2:7516122/1‑150 (MQ=255)
                                                                                                                                            cTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAcc          >  2:2086966/1‑150 (MQ=255)
                                                                                                                                             tGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGccc             <  2:8640242/146‑1 (MQ=255)
                                                                                                                                             tGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGccc             >  1:8640242/1‑146 (MQ=255)
                                                                                                                                              ggCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCt        >  2:3454612/1‑150 (MQ=255)
                                                                                                                                               gCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCtt       <  1:5172179/150‑1 (MQ=255)
                                                                                                                                                 aCCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTTgg     >  1:5389039/1‑150 (MQ=255)
                                                                                                                                                    aTCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCtcttc                                        <  2:5983957/112‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGGAAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCACCCTTGGAAC  >  NC_002947/566479‑566776

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: