Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 507,782 T→G 23.1% H100P (CAC→CCC)  ybaP ← TraB family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913507,7820TG23.1% 94.9 / 2.1 57H100P (CAC→CCC) ybaPTraB family protein
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base G (6/10);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCAGGGAGCTGGAGCAACATGGCAATCTGGTTTTCAGCCCCTTCCAGTTCAATCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCG  >  NC_000913/507561‑508009
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                      
tggtaGGGTGGTGGATCAACTTGGCATTCTGGTTTTCAGCCCCTTCGTGGTCAACCGCGGGTTCATTTTTTTGACGCGCCGCCTGCAATAGCGGGTAATAGACACAGGATTCTGCCCGGAGCCCCAGTTTTTGTGCCTGAGTAGTCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGtgctgc                                                                                                                                                                                                                                   <  1:704144/223‑1 (MQ=255)
tGTCAGGGAGCTGGAGCAACATGGCAATCTGGTTTTCAGCCCCTTCCAGTTCAATCACGGGTTTATGTTGTTGCTTCGCCACCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTGGAAAAGGGGGAGGGAGAAATGCCCATCCCCTGGCTAATGTGCTGCa                                                                                                                                                                                                                                  >  2:657233/1‑228 (MQ=255)
           tGGAGCAACATGGCAATCTGGTTTTCAGCCCCTTCCAGTTCAATCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATGTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGt                                                                                                                                                                                                                       <  2:556956/228‑1 (MQ=255)
            ggAGCAACATGGCAATCTGGTTTTCAGCCCCTTCCAGTTCAATCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGGAAAGAGGGAGGGAGAAATGCCCACAACCTGGTTAATGGGCTGCAGGTTTTGGAGtt                                                                                                                                                                                                                      >  2:670352/1‑228 (MQ=255)
                  aataggcaagtgtgTTTCCAGCCCCTTACAGTGTAATCACGGGTTTATGGTTTTGCTACGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTGTTGGAGTTGctc                                                                                                                                                                                                                 <  2:467828/213‑1 (MQ=255)
                                         ccTTCCAGTTCAATCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGGGCTGAAGGGTTTGGAGTTGCTCCTCGCTAATCGCGTCTTCCAgcgc                                                                                                                                                                                        >  2:528857/1‑227 (MQ=255)
                                                      tCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGTTAGAAAAGAGTGAGGGAGAAAGGCCCATCTCCTGGCTAATGGTCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCGCTTCCAGCGCCTCGACggcgggca                                                                                                                                                                          >  2:328146/1‑228 (MQ=255)
                                                      tCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCTTCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAggcaggca                                                                                                                                                                          >  2:130243/1‑228 (MQ=255)
                                                       cagggttttagtttgttgCTTCGCCGCGTGCAATAGCGGGTAAGCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCaa                                                                                                                                                                         <  2:276634/216‑1 (MQ=255)
                                                           ggTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGTGGAAATGCCCATCTCCTGGCTAATGGGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCCCCTCGCAGGCGGGCAAATTa                                                                                                                                                                     >  1:129151/1‑228 (MQ=255)
                                                              ttttggtgttgCGTCGCCGCCTGGAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCa                                                                                                                                                                  <  2:671095/222‑1 (MQ=255)
                                                               tttgttttTGGTGCGCCCCCGGCAAGAGCTGGTAAGCGATCCCGGATTCTGCCCGCAGCCCCAGTTTTTGTTCCTGCGTAGGCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCaa                                                                                                                                                                 <  1:541768/226‑1 (MQ=255)
                                                                      ttgCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACCGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGt                                                                                                                                                          >  2:121399/1‑228 (MQ=255)
                                                                      ttgCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGCAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGGTCCTCGCTAATGCGCTCTTCCAGCGCCTCGACGGCAGGCAAGTTAGCAAAAGGCGt                                                                                                                                                          >  2:677512/1‑228 (MQ=255)
                                                                      ttgCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTATGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGACTGAAGAACAATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCCGGCTAATGGGCTGCAGGTTTTGGAGGTGCGCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCCAATTAGCCAAAGACGt                                                                                                                                                          >  2:483645/1‑228 (MQ=255)
                                                                           tCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGGGATTTGCCACAGCGGTTGGGTAGAAAAGAGTAGGGGGGAAATGCCCATCCCCCGGCTAATGGGCTGCAGGGTTTGGAGTTGGTCCTCGCTAATGCGCCCTTCCACCGCCCGGCAGGCGGGGAAATTAGCACAAgagggatcg                                                                                                                                                      >  1:234475/1‑219 (MQ=255)
                                                                           tCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGAATCGc                                                                                                                                                     >  2:566059/1‑228 (MQ=255)
                                                                                   tGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGtgc                                                                                                                                                                                                                                      >  1:143421/1‑141 (MQ=255)
                                                                                    gCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGtgct                                                                                                                                                                                                                                     <  2:143421/141‑1 (MQ=255)
                                                                                             ggTAATCGATACCGTATTCTGCACGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAGGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCCGGGTTTGGAGTTGCTCCTCCCTAATTCGCTCTTCCAGCTCCTGCCAGGCAGGCAAATTACCAAAAGGCGTATCGCTCGGGGACAAATcgcccc                                                                                                                                  >  1:240131/1‑223 (MQ=255)
                                                                                               taacggATACCTTATTCTGCCCGCAGCCCCAGTTGTTGGGCCTGGGGCGCCTGAAGAACCATCGCGAATTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTcg                                                                                                                                <  2:212789/223‑1 (MQ=255)
                                                                                                        ccGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTACCAAAAGTCGGATCGCTGGTGGAACACTCCGCCTCGACGATAAgc                                                                                                                       >  2:672477/1‑228 (MQ=255)
                                                                                                                tGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGAGATTTGCCACAGCGGTTGGGGAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCCGGCAGGCAAATTAGCAATAGGAGTATCGCTGGTGGAAACATCGGCCTCGACGCTCAGCGCATCGGc                                                                                                               >  2:527711/1‑228 (MQ=255)
                                                                                                                  cccGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGTAAACATCCGCCTCGAGGATCAGAGCATCGGCGt                                                                                                             >  2:480314/1‑228 (MQ=255)
                                                                                                                   ccGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGGGTTGCTCCTCGCTAATGGGCCCTTCCAGCGCCTCGCAGGCCGGCAAATTAGCAAAAGGCGTATCGTCGGGGGAAACCTCCCCCCCGACGATCAGCGAACCAGCGtt                                                                                                            >  1:686695/1‑228 (MQ=255)
                                                                                                                   ccGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGCTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGCCGATCAGCGCATCGGCGtt                                                                                                            >  1:482964/1‑228 (MQ=255)
                                                                                                                       aGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGCCGATCAGCGCATCGGCGtttttg                                                                                                        >  2:336881/1‑228 (MQ=255)
                                                                                                                               tgtttGTGCCTGGCGCGCCCGAAGAACCATCGCGATTTGCCACAGGGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATTACCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGtttttgagcttttt                                                                                                <  2:522067/226‑1 (MQ=255)
                                                                                                                                         tGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGGGCAAACgc                                                                                      >  2:160110/1‑227 (MQ=255)
                                                                                                                                           gcaggCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACggg                                                                                     <  2:234475/223‑1 (MQ=255)
                                                                                                                                           cGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTACTGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAGAAGCCGTATCGCCGGTGGAAGCATCAGCCTCGACGATCAGCGCATCGGCGCTTTTGAGCTTTTTTAGCAAACggc                                                                                     >  2:245207/1‑226 (MQ=255)
                                                                                                                                           cGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTAGCAGGCAGGCAAATTAGCAAAAGGCGGATCGCTGGGTGGAAACATCCGCCTCGCCGATCAGCGCATCGGCGGTTTTGAGCTTTTTGAGCAAACgg                                                                                      >  2:517336/1‑227 (MQ=255)
                                                                                                                                                   taaaaaCCTTCGCGATTTTCCACAGCGGTTGGGTAGAAAAGTGTGGGGGATAAAGCCCCATCTGCGGGCTAATGTGCGGCAGGTTTTGGAGTTGCTCCTCGCAAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTTTCGCTGGTGGAAACACCCGCCTCGGCGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAggg                                                                            <  1:421487/224‑1 (MQ=255)
                                                                                                                                                   gaagtaCCATCGCGTTTTGCCGCTGGGGTGGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAAGGTGCTGCAGGGTTTGGAGTTGCTCCTCGCTAATGCGCCCTGCCAGCGCATCGCAGGCAGGCAAATTAGCAAAAGGGGTATCGTTGGTTGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAggg                                                                            <  1:355854/228‑1 (MQ=255)
                                                                                                                                                    aagcaccaacctgtttttCCACGGGGGTTGGGTGAAAAAGAGGTGGGGAGAAAGGCCCATCTCCGGGCTAAGGGGCGGCAGGGTTTGGAGTTGCTCCTCGATAATGGGCTCTTCCGGCGCCTCGCAGGCAGGGAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGGCGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAggg                                                                            <  1:657233/213‑1 (MQ=255)
                                                                                                                                                     cgaaCCATCGCGATTTGCCACAGAGGTTGGGAAGAAAAGAGGGAGGGAGAAATGCCCATCTCCTGGCTAAGGTGCTGCAGGCTTTGGAGTTGCACCTCGCTAATGCGCTCTTCAAGCGCCTAGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGGAAAATCGGCCTCACCGATCAGAGAAACGGGGTTTTTGAGCTTTGTGATCAAACGGTGGGGGAggggg                                                                          >  1:450043/2‑227 (MQ=255)
                                                                                                                                                       aaCCATCGGGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAg                                                                         <  2:150171/227‑1 (MQ=255)
                                                                                                                                                             cgcgATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGCTCAGCGCATAGGCGTTTTTGAGCTTTTTGTGCCAACGGGTGGGCAGGGGAGCCATATc                                                                  >  2:615558/1‑228 (MQ=255)
                                                                                                                                                              gcgATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGCCAGGGGGGCCATATCg                                                                 >  2:481996/1‑228 (MQ=255)
                                                                                                                                                                   ttGCCACAGCGTTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTgg                                                              <  2:473315/226‑1 (MQ=255)
                                                                                                                                                                   tggCCACAGCTGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCGTCTCCTGGCTAATGGGCTGCAGGTTATGGAGTTGGTCCTCGCTAATGCGCTCTTCCAGCGCCGCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTgg                                                              <  2:30880/224‑1 (MQ=255)
                                                                                                                                                                    tGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAACCATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCACACGGGTgg                                                                                  >  2:411820/1‑205 (MQ=255)
                                                                                                                                                                    tGCCACAGCGGTGGGGTAGAAAAGGGTGAGGGGGAAATGCCCATCTCCTGGCTAATGTGCGGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTa                                                           <  2:198466/228‑1 (MQ=255)
                                                                                                                                                                     tataagaggaaggggTAGAAAAGAGGGGGGGATAAATCCCCCTCGCCGGGGTAATGTGCTGCAGGTTTTGGAGTTGTCCCGCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGGGGAACCATCCGCCTCGACGAGCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTa                                                           <  1:454721/215‑1 (MQ=255)
                                                                                                                                                                     gCCACGGGGGTTGGGTAAGAAAGGGGTAGGGAGAAAGCCCCATCTCCTGGCTAATGGGCTGCAGGTTTGGGAGTTGCTCCCCGGTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGGTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTggg                                                                                 <  1:411820/205‑1 (MQ=255)
                                                                                                                                                                      tccgaGCGGTTGGTTGGAAAAGGGTGGGGGGGAAATGCACCTTGCCGGGTTAATAGGCTGGAGGTGTTGGAGTTGCTCATGGCTAATGCAGTCTTCCGGCGCCTCGCAGGCAGGAAAATTAGCAAAAGGGGTATCGGTGGTGGAACCAGCCGCCTGGACGCTCAGCGCAGCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTAcc                                                         <  1:566059/224‑1 (MQ=255)
                                                                                                                                                                      ctctgggggTTGGGTGGGAATGAGGGGGGGAGAATCGCCCATCTCGTGCTTATGGTGCTGCAGGTTGTGGAGTTGCTCCTGGATAAGGCGCTCTTCCAGCGCCTCGCAGGCTGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAACCATCCGCCTCGACGATCAGCGCAGCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTAcc                                                         <  1:483645/221‑1 (MQ=255)
                                                                                                                                                                      ccACAGCGGTTGGGTAGAAAGGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTAcc                                                         <  2:570641/228‑1 (MQ=255)
                                                                                                                                                                         ctggggTGGGGTAGATAAGGGTGAGGGAAAAATGCCCATTCCCTGGTTAAGGGGCGGCAGGTTTGGGATGTGCTCCTCGCTAATGCGCTCTTCCGGCGCCTCGCAGGCAGGTAAATTAGCAAAAGGCGTATCGCTGGTGGAACCATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCat                                                      <  1:458240/224‑1 (MQ=255)
                                                                                                                                                                         caGCGGTGGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGATAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCat                                                      <  2:83250/228‑1 (MQ=255)
                                                                                                                                                                          aGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTCGGGGAACCATCCCCCTCGCCGATCACGGAATAGGCGTTTTTTAGGGTGTTGAGCAACCGGGTGGGCAGGGGGGcactatgatggggcccaaaa                                                     >  1:246626/1‑209 (MQ=255)
                                                                                                                                                                              gTTGGGTAGAAAAGAGTGAGGGAGAAAGGCCCATCTCCTGGCTAATGGGCTGCAGGTTTTGGAGTTGATCCTCGCTAAGGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGaa                                                 <  2:380137/228‑1 (MQ=255)
                                                                                                                                                                                        aaaGGGTGGGGGAAAAATCCCCATCTCCTGGCTAATGGGCGGCAGGTTTTGGGTTTGCTCCTCGATAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACAGCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAAt                                       <  1:481996/228‑1 (MQ=255)
                                                                                                                                                                                         aaGAGTGAGGGAGAAATGCCCATCTCTGGGCGAATGTGCTGCAGGTTTTGGAGTGGCTCCTAGCTAATGCGATCTGCCAGCGCCTCGCAGGCAGGTAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATAGCGTGGCTACCCATATGAATACTGCCAAt                                       <  2:527939/227‑1 (MQ=255)
                                                                                                                                                                                            attgggggAGAAATGCCCATCTCCGGGCTAATGGGCTGCAGGTTTTGGAGTTGCTCCGCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGa                                   <  2:392091/223‑1 (MQ=255)
                                                                                                                                                                                            agTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAAtctgtctcttatc                                    >  2:454721/1‑215 (MQ=255)
                                                                                                                                                                                             gTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGAGCATAGGCGTTTTTGAGCTTTTTGAGCAACCGGGTGGGCAGGGGAGCCCGATCGGGGCTACCCATATGAATACCGcaatacacag                                  >  1:657888/1‑219 (MQ=255)
                                                                                                                                                                                                        atccccATCCCCTGGCTAATGTGCTGCAGGTTTGGGAGTTGCCCCTCGCAAAGGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCTTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGc                         <  2:540639/223‑1 (MQ=255)
                                                                                                                                                                                                        acgCCCCTCTCCGTGGTAATGGGCGTCGGGTTTTGGATTTGCTCCGAGCTAACGTGTTTGCCCAGCGCCGCGCCGGCAGGGAAATTAGCAAAAGGCGTATCGCGGGTGGAAACATCCGCCTCGACGACCGGCGCAGCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGa                       <  1:677512/226‑1 (MQ=255)
                                                                                                                                                                                                         cccccAGCTCCTGGCTGAGTGGCTGCGGGTTGTAGGATGCGTCCTAGCTAATGCGCTCTTACAGCGCCTCGCAGGCAGGCAAATTAGCAAAGGGCGTGTCGGTGGTGTAAACATCCGCCTCGACGATCAGCGCATCGGAGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGa                       <  1:506594/225‑1 (MQ=255)
                                                                                                                                                                                                                           tgtgCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCCAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGAGGAAAATGGCGAATGCAGGGGGGGGGGaaaa    >  1:12232/1‑225 (MQ=255)
                                                                                                                                                                                                                             tgctgcAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCCTATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGCGGGGGataccc  >  2:487854/1‑225 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                      
TGTCAGGGAGCTGGAGCAACATGGCAATCTGGTTTTCAGCCCCTTCCAGTTCAATCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGCTCCTCGCTAATGCGCTCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCTGGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCG  >  NC_000913/507561‑508009

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: