Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,600,215 +C intergenic (‑26/+23) yjjA ← / ← dnaC DUF2501 domain‑containing protein YjjA/DNA replication protein DnaC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,600,2141.C97.0% 254.2 / 0.3 67intergenic (‑25/+24)yjjA/dnaCDUF2501 domain‑containing protein YjjA/DNA replication protein DnaC
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (37/28);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

CATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCT  >  NC_000913/4599967‑4600455
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
gggcgtgcgtgtggtcaccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCgcc                                                                                                                                                                                                                                                  <  1:267364‑M2/9‑2 (MQ=255)
      gcgtgtggtcaccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATAc                                                                                                                                                                                                                                           <  1:676869‑M2/16‑1 (MQ=255)
            ggtcaccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAg                                                                                                                                                                                                                                     <  1:473884‑M2/22‑1 (MQ=255)
                  cgttatcagtttcggttttgtcttcttcgtccttgttcagatctgcggccccttgacgttggtttttggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagggaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATCCCGCCAGAAACgg                                                                                                                                                                                                                               <  2:791649‑M2/28‑1 (MQ=255)
                   ggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagatccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACggg                                                                                                                                                                                                                              >  2:199928‑M2/223‑251 (MQ=255)
                           tttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAAt                                                                                                                                                                                                                       >  1:283755‑M2/215‑250 (MQ=255)
                            ttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATa                                                                                                                                                                                                                      >  2:245082‑M2/214‑250 (MQ=255)
                                                     accagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaacgaatgaacgaataaagggagcgcaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGt                                                                                                                                                                                            >  2:685175‑M2/189‑251 (MQ=255)
                                                      ccagttcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTa                                                                                                                                                                                           <  1:487502‑M2/64‑1 (MQ=255)
                                                                ggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttcatcagcctgtggctgttaaagtgaaaatcgtaaccgaatgaacgcataaagggagcgaaagatccctttattattgttagcaaagtgtgattcgttcattcctgaTGTGTCGCGAGGATACCGCCAGAAACGGGGATTAATACTATTTACCTGTTACCCGGCTACGGTAGCTATccc                                                                                                                                                                                   >  2:558697‑M2/178‑249 (MQ=255)
                                                                   ccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTg                                                                                                                                                                              >  1:198694‑M2/175‑251 (MQ=255)
                                                                            ttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCAcc                                                                                                                                                                     >  2:596791‑M2/166‑251 (MQ=255)
                                                                                  tatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACaa                                                                                                                                                                <  1:219902‑M2/91‑1 (MQ=255)
                                                                                        atcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGtt                                                                                                                                                          >  2:11688‑M2/154‑250 (MQ=255)
                                                                                            cgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTAccc                                                                                                                                                      >  1:705992‑M2/150‑250 (MQ=255)
                                                                                                       ggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACgg                                                                                                                                          >  2:39131‑M2/139‑251 (MQ=255)
                                                                                                        gcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTg                                                                                                                                                                                                           >  2:702401‑M2/138‑185 (MQ=255)
                                                                                                        gcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTg                                                                                                                                                                                                           <  1:702401‑M2/48‑1 (MQ=255)
                                                                                                         cgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATa                                                                                                                                             >  2:630402‑M2/137‑246 (MQ=255)
                                                                                                         cgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATa                                                                                                                                             <  1:630402‑M2/110‑1 (MQ=255)
                                                                                                          gcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTcc                                                                                                                                       >  1:295121‑M2/136‑251 (MQ=255)
                                                                                                            actgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCAt                                                                                                                                     >  2:473384‑M2/134‑251 (MQ=255)
                                                                                                            actgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCAt                                                                                                                                     >  1:772886‑M2/134‑251 (MQ=255)
                                                                                                            actgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCAt                                                                                                                                     <  1:312688‑M2/118‑1 (MQ=255)
                                                                                                               gttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCa                                                                                                                                   <  2:123754‑M2/120‑1 (MQ=255)
                                                                                                               gttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCAc                                                                                                                                  <  1:282415‑M2/121‑1 (MQ=255)
                                                                                                               cGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCAcg                                                                                                                                 >  1:657251/1‑251 (MQ=255)
                                                                                                                    caagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGtt                                                                                                                             >  1:422628‑M2/126‑251 (MQ=255)
                                                                                                                            agcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGc                                                                                                                     <  1:685175‑M2/134‑1 (MQ=255)
                                                                                                                                       aatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCa                                                                                                           <  1:11688‑M2/144‑1 (MQ=255)
                                                                                                                                            gcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCtt                                                                                                     >  1:616564‑M2/102‑251 (MQ=255)
                                                                                                                                                gtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTcc                                                                                                   <  2:582800‑M2/152‑1 (MQ=255)
                                                                                                                                                 tggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCat                                                                                                 >  1:749479‑M2/97‑250 (MQ=255)
                                                                                                                                                  ggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCatat                                                                                               >  1:283818‑M2/96‑251 (MQ=255)
                                                                                                                                                    ctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTAcc                                                                                                                                                       <  2:42966‑M2/100‑1 (MQ=255)
                                                                                                                                                    ctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTAcc                                                                                                                                                       >  1:42966‑M2/94‑193 (MQ=255)
                                                                                                                                                      gttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGTGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCGCTTCCATATTa                                                                                             <  2:40084‑M2/158‑1 (MQ=255)
                                                                                                                                                               aaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTc                                                                                    >  1:158429‑M2/83‑249 (MQ=255)
                                                                                                                                                               aaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAg                                                                                  <  1:515182‑M2/169‑1 (MQ=255)
                                                                                                                                                                aaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGc                                                                                 <  2:164142‑M2/170‑1 (MQ=255)
                                                                                                                                                                    cgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTTCCCAGGCGCATACGGTCCATCTCGCGTTCGCCCCGCAACTTGGTCATCTCTTCCATATTACTGGTGGTCAGCATcc                                                                             <  2:90293‑M2/174‑1 (MQ=255)
                                                                                                                                                                     gtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATccc                                                                            <  2:663153‑M2/175‑1 (MQ=255)
                                                                                                                                                                           gaatgaacgaatcaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTTTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTgg                                                                      >  1:563761‑M2/71‑251 (MQ=255)
                                                                                                                                                                                aacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCcac                                                                                                                                                                  >  1:410544‑M2/66‑154 (MQ=255)
                                                                                                                                                                                aacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCcac                                                                                                                                                                  <  2:410544‑M2/89‑1 (MQ=255)
                                                                                                                                                                                   gaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTgg                                                              >  2:581280‑M2/63‑251 (MQ=255)
                                                                                                                                                                                    aataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGTCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTgg                                                              >  1:77537‑M2/62‑250 (MQ=255)
                                                                                                                                                                                           ggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGaagaag                                                        >  1:285897‑M2/55‑249 (MQ=255)
                                                                                                                                                                                            gagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgc                                                     >  1:480579‑M2/54‑251 (MQ=255)
                                                                                                                                                                                               cgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcggc                                                  >  2:416398‑M2/51‑251 (MQ=255)
                                                                                                                                                                                                gaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTc                                                                                                    <  1:233966‑M2/151‑1 (MQ=255)
                                                                                                                                                                                                gaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTc                                                                                                    >  2:233966‑M2/50‑200 (MQ=255)
                                                                                                                                                                                                gaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTCCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcggcg                                                 <  2:283755‑M2/202‑1 (MQ=255)
                                                                                                                                                                                                  aagctccctttattattgttagcaaagggtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGAt                                               <  2:780236‑M2/204‑1 (MQ=255)
                                                                                                                                                                                                   agctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATc                                              >  2:782820‑M2/47‑251 (MQ=255)
                                                                                                                                                                                                     ctcactttattattgttagcaacgtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGACGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGCTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATcg                                             >  2:43445‑M2/45‑250 (MQ=255)
                                                                                                                                                                                                          tttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCTACGATc                                       >  2:187491‑M2/40‑251 (MQ=255)
                                                                                                                                                                                                          tttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATc                                       >  1:785223‑M2/40‑251 (MQ=255)
                                                                                                                                                                                                                attgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGtt                                  <  1:596791‑M2/217‑1 (MQ=255)
                                                                                                                                                                                                                 ttgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtg                                 >  1:661548‑M2/33‑250 (MQ=255)
                                                                                                                                                                                                                              gtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTc                   <  1:416398‑M2/232‑1 (MQ=255)
                                                                                                                                                                                                                                      ttcattcctgaTGTGTTGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCCGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCCGCATCCCGGTTGGGCGTTTGGAAGCAGAGCGGCGATCGACGCACTGGTTGAGGATCACTTTTTCGTATTTCg           >  2:369555‑M2/12‑251 (MQ=255)
                                                                                                                                                                                                                                           tcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTc       <  1:245082‑M2/244‑1 (MQ=255)
                                                                                                                                                                                                                                           tcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCg      >  2:227621‑M2/7‑251 (MQ=255)
                                                                                                                                                                                                                                            cctgaTGTTTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCgg     <  2:563761‑M2/246‑1 (MQ=255)
                                                                                                                                                                                                                                            cctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCgg     <  1:473384‑M2/246‑1 (MQ=255)
                                                                                                                                                                                                                                            cctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCgg     >  1:614981‑M2/6‑251 (MQ=255)
                                                                                                                                                                                                                                               gaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCt  >  1:419440‑M2/3‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
CATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCT  >  NC_000913/4599967‑4600455

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