Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,041,467 A→G K31K (AAA→AAG zinT → metal‑binding protein ZinT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,041,4670AG98.6% 246.5 / ‑5.6 71K31K (AAA→AAGzinTmetal‑binding protein ZinT
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (25/45);  minor base T (1/0);  total (26/45)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATT  >  NC_000913/2041222‑2041698
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
cgcggtcttgctgggtacgggactgggtggagaaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCt                                                                                                                                                                                                                                     <  1:162547‑M2/93‑1 (MQ=255)
  cggtcgtgctgggtacgggactgggtggagaaacattgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtacgggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGGCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCTCACGGCAAGCCCTTaa                                                                                                                                                                                                                                  <  2:348952‑M2/96‑1 (MQ=255)
   ggtcttgctgggtacgggactgggtggagaaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTTGCATGGTCATCACTCACACGGCAAGCCCTTAAc                                                                                                                                                                                                                                 >  2:95639‑M2/155‑251 (MQ=255)
            gggtacgggactgggtggagaaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatctCGATTCGTCTTGACAATCGGGCTGTTGCTTTAGGTGTCTTTATTGTGAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTc                                                                                                                                                                                                                          <  1:738533‑M2/104‑1 (MQ=255)
            gggtacgggactgggtggagaaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTc                                                                                                                                                                                                                          <  1:220247‑M2/104‑1 (MQ=255)
                acgggactgggtggagaaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACaa                                                                                                                                                                                                                    <  1:30162‑M2/110‑1 (MQ=255)
                          gtggagaaacactgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCGTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTa                                                                                                                                                                                                          <  2:617688‑M2/120‑1 (MQ=255)
                                  acactgtcgaagattgaatggctggcgctcggcgtaattgtattcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgacagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGtt                                                                                                                                                                                                  >  1:231182‑M2/124‑251 (MQ=255)
                                     ctgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCTTTAACAGAGGACGAACAAAAAGCTGCTAATGGTGttttt                                                                                                                                                                                               >  1:7068‑M2/121‑251 (MQ=255)
                                     ctgtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGttttt                                                                                                                                                                                               >  1:705683‑M2/121‑251 (MQ=255)
                                       gtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTtg                                                                                                                                                                                              <  1:505596‑M2/132‑1 (MQ=255)
                                       gtcgaagattgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTtg                                                                                                                                                                                              <  2:424522‑M2/132‑1 (MQ=255)
                                           aagattgaatggctggcgctcggcgtaattgtcttcgcgttggtactggtcacgttgggaaaatagctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCTATTCGTCTTTACAAACTGGCTGTTGCTTTAGATGTCTTTATTGTTAGCGCTCCTGCCTTTGCGCATGGTCATCACTCACTCGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGGTTTTtatgat                                                                                                                                                                                         <  2:234760‑M2/137‑1 (MQ=255)
                                              attgaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccatttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGcc                                                                                                                                                                                      <  2:596497‑M2/140‑1 (MQ=255)
                                                 gaatggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacacgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAAt                                                                                                                                                                                   >  2:786474‑M2/109‑251 (MQ=255)
                                                    tggctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGt                                                                                                                                                                                 >  1:707061‑M2/106‑250 (MQ=255)
                                                      gctggcgctcggcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACa                                                                                                                                                                              <  1:644167‑M2/148‑1 (MQ=255)
                                                                 gcgtaattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGc                                                                                                                                                                   >  2:297550‑M2/93‑251 (MQ=255)
                                                                  cgtaattgtcttcgcggtggtactggtcacgttgggaaaagatctcttcccggctaaacccgtagtttcgccagttattcaggacgcatcaCGATTCTTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAATAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCt                                                                                                                                                                  <  2:311960‑M2/160‑1 (MQ=255)
                                                                     aattgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGGGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAATCCCTTAACATAGTTCTAACACAACTCTTCTCATTTTTTTTTTTATTATGCCAATGTACAAAACCTAACTCTcat                                                                                                                                                               >  2:507439‑M2/89‑250 (MQ=255)
                                                                      attgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGt                                                                                                                                                              >  1:483804‑M2/88‑251 (MQ=255)
                                                                       ttgtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCAAATGGACAAAACCGAACGCTCAGTg                                                                                                                                                             >  2:480061‑M2/87‑251 (MQ=255)
                                                                         gtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGAc                                                                                                                                                           >  1:498415‑M2/85‑251 (MQ=255)
                                                                         gtcttcgcggtggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGAc                                                                                                                                                           <  1:460444‑M2/167‑1 (MQ=255)
                                                                                 gttggtgctggtcacgtttggaaattatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATgg                                                                                                                                                   <  2:127817‑M2/175‑1 (MQ=255)
                                                                                    ggtactggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGt                                                                                                                                                >  1:102285‑M2/74‑251 (MQ=255)
                                                                                    ggtaatggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGt                                                                                                                                                <  2:737832‑M2/178‑1 (MQ=255)
                                                                                          ggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGc                                                                                                                                           <  1:786474‑M2/183‑1 (MQ=255)
                                                                                          ggtcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCa                                                                                                                                          <  2:131577‑M2/184‑1 (MQ=255)
                                                                                          gggcacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAGCTGGCTGTTGCTTTAGGTGTCTTTATTGGTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCGCAGTGACTGGGATGGAGTCTGGc                                                                                                                                           <  2:695550‑M2/183‑1 (MQ=255)
                                                                                             cacgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATc                                                                                                                                       >  1:46951‑M2/65‑251 (MQ=255)
                                                                                              acgttgggaaaatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATcc                                                                                                                                      >  1:612887‑M2/64‑251 (MQ=255)
                                                                                                    ggaaaatatctcttcccggccaaacccgtagttgctccatttattcaggacgcatcaCGTTTCGTCTTTACAAACTGGCTGTTGTTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTTCAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTAt                                                                                                                                <  2:528962‑M2/194‑1 (MQ=255)
                                                                                                        aatatctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCtt                                                                                                                            <  2:707061‑M2/198‑1 (MQ=255)
                                                                                                           atctcttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTa                                                                                                                          <  1:455335‑M2/200‑1 (MQ=255)
                                                                                                               cttcccggcaaaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACt                                                                                                                        <  1:685678‑M2/202‑1 (MQ=255)
                                                                                                                         aaacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAATCCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCaaa                                                                                                           <  2:825427‑M2/215‑1 (MQ=255)
                                                                                                                          aacccgtagttgcgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACACAAAGCTGCTAATGGTGTTTTTGATGCTGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAAc                                                                                                          >  2:878306‑M2/36‑251 (MQ=255)
                                                                                                                                      cgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAATGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCt                                                                                              <  1:465300‑M2/228‑1 (MQ=255)
                                                                                                                                      cgccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCt                                                                                              <  1:869318‑M2/228‑1 (MQ=255)
                                                                                                                                        ccagttattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCttt                                                                                            <  2:46951‑M2/230‑1 (MQ=255)
                                                                                                                                             tattcaggacgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCaga                                                                                        <  1:480061‑M2/234‑1 (MQ=255)
                                                                                                                                                      cgcatcaCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCgg                                                                               >  1:69405‑M2/8‑250 (MQ=255)
                                                                                                                                                       ccatcaCGATTCGCCTTTACAACCGTGCTGTTGCTTTCGGTGTCTTTATTGTTATCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGATGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGAt                                                                             <  2:514700‑M2/245‑1 (MQ=255)
                                                                                                                                                           tacGATTCGTCTTTACACACTGTCTTTTGCTTTCGTTTACTTTATTGTTAGCGCTCCTGCCTATTCGCATGGTCATCTCTCACGCGTCAAGCCCTTACCAGATGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGTTGCCAATGTACAAAACCGAACGCTCAGTGACTGTGATGCAGTCTGGAAATCCTTTTAAAATTTACTGCAAATTCTCAAACTTGACCCCGTCTTTCAGCAAAAAGCGGATgg                                                                           <  2:138858/247‑2 (MQ=255)
                                                                                                                                                           catgATTCGTCTTTTCAAACTGGCTGTTGCTTTTGGTGTCTTTATTGTTAGCGCTCCTGCCTGTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCTAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTTGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAACCTTGACCCCGTCTTTCAGAAGAAAGCGGTTGCAg                                                                         <  2:698662/248‑1 (MQ=255)
                                                                                                                                                             cGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAt                                                                       <  2:497526/251‑1 (MQ=255)
                                                                                                                                                              gATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAta                                                                      <  1:444013/251‑1 (MQ=255)
                                                                                                                                                                  cGTCTTTACAAACTTGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTTTCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAtaaaac                                                                  <  1:468159/251‑1 (MQ=255)
                                                                                                                                                                  cGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTTTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAtaaaa                                                                   <  1:594746/250‑1 (MQ=255)
                                                                                                                                                                   gTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAtaaaact                                                                 <  1:297550/251‑1 (MQ=255)
                                                                                                                                                                   gTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAtaaaact                                                                 >  1:631088/1‑251 (MQ=255)
                                                                                                                                                                      tttACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAtaaaactaa                                                               >  1:509131/1‑250 (MQ=255)
                                                                                                                                                                       ttACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAtaaaactaaa                                                              <  1:343398/250‑1 (MQ=255)
                                                                                                                                                                           aaaCTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACAtt                                                         >  1:127686/1‑251 (MQ=255)
                                                                                                                                                                           aaaCTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACAtt                                                         <  2:631088/251‑1 (MQ=255)
                                                                                                                                                                                ggCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCt                                                     <  1:479885/250‑1 (MQ=255)
                                                                                                                                                                                ggCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACAGTTGCt                                                     <  1:878306/250‑1 (MQ=255)
                                                                                                                                                                                ggCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCCTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAACGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGAGGCAACCTTGACCCCGTCTTTCAGAAGCAAGCGGATGCAGATAAAACTAAAAAATTTGCTg                                                    >  2:572184/1‑251 (MQ=255)
                                                                                                                                                                                               tgtCTTTATGGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACAGTACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGCTGCAGATAAAACTAAACCATTTGCTGAAATTAAAGattatt                                     >  2:15182/1‑251 (MQ=255)
                                                                                                                                                                                                gtCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGattatt                                     >  1:411422/1‑250 (MQ=255)
                                                                                                                                                                                                    ttATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATcac                                >  2:847585/1‑251 (MQ=255)
                                                                                                                                                                                                    ttATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATca                                 <  1:812850/250‑1 (MQ=255)
                                                                                                                                                                                                        ggTTAGCGCTCCTGCCTTTTCGCATGGTCACCACTCACACGGCTAGCCCGTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGGCAGTGTACAAAACCGAACGCGCAGTGACTGTGATGGAGTCTTGCAATCCGTTTATCCTTTACTGCAAAGTAGCAAACTTGACCCCGTCTTACAGAAAAAATCGGATGCAAATAAAACTAAAACATTTGCTGAAATTAAATATTATTATCACAAAg                            <  2:628766/250‑1 (MQ=255)
                                                                                                                                                                                                               gcTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGc                     <  2:705683/251‑1 (MQ=255)
                                                                                                                                                                                                                         tttCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATc           >  1:659164/1‑251 (MQ=255)
                                                                                                                                                                                                                          ttCGCATGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCg          >  2:331300/1‑251 (MQ=255)
                                                                                                                                                                                                                             gcaagttcatcagccacaGGGCAAGCTCTTAACAGAGTTCGACTAACAAGATGCTCCTGGAGTTTTTCATGATGCCAATGTACAAAACCGAATGCTCAGCGACTGGGATGGAGTCTTGAAATCCGTTTCTCCTTTACTGCCACGCGGCAAACTTGACCCCGCCTTTCAGCAGGAAGCGGATGCAGATAAATCTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGAAACTGATATCgggg       >  2:763133/15‑248 (MQ=255)
                                                                                                                                                                                                                               aTGGTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAgatg     <  2:18211/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gTCATCGCTCACACGGCAAGCGCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGGGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTGTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAATTAAAACATTTGCGGAAATTAAAGATTATTATCACAAAGGTTATGGAACAGATATCGAGATGAtt  <  2:648010/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gTCATCACTCACACGGCAAGCCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGAtt  <  1:387200/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
CCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAAAGATTATTATCACAAAGGTTATGCAACAGATATCGAGATGATT  >  NC_000913/2041222‑2041698

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: