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.C96.4% 187.7 / 1.2 56intergenic (‑25/+24)yjjA/dnaCDUF2501 domain‑containing protein YjjA/DNA replication protein DnaC
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (26/28);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

AAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCG‑C‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTGC  >  NC_000913/4599990‑4600457
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
ccagtttcggttttgtcttctacgtactggaccagatcttcggcccgctaacgttggtttagggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcatcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgctgcgtgcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGAt                                                                                                                                                                                                                             <  1:26150‑M2/32‑1 (MQ=255)
 cagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTa                                                                                                                                                                                                                           <  1:59309‑M2/34‑1 (MQ=255)
  agtttcggttttgtcttcaacgtacgggaccagttcttcggcccgcggacgtgggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTaa                                                                                                                                                                                                                          <  2:183565‑M2/35‑1 (MQ=255)
   gtttcggttttgtcttctacgtcctggcccagatcttcggcccgcggacgttgggttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAAt                                                                                                                                                                                                                         <  2:232002‑M2/36‑1 (MQ=255)
        ggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagatcactttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATActct                                                                                                                                                                                                                    >  2:343426‑M2/211‑251 (MQ=255)
         gttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCtt                                                                                                                                                                                                                   >  1:407276‑M2/210‑251 (MQ=255)
               tcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatccagccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTg                                                                                                                                                                                                             >  2:187933‑M2/204‑251 (MQ=255)
                  tctccgtactggaccagttcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGtt                                                                                                                                                                                                           <  1:306802‑M2/50‑1 (MQ=255)
                    tacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAACCGGGGATTAATACTATTTACCTGTTAcc                                                                                                                                                                                                        >  2:345523‑M2/199‑251 (MQ=255)
                      cgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATGCTCTTTACCTGTTACCCg                                                                                                                                                                                                      <  1:289406‑M2/55‑1 (MQ=255)
                                    tcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTAt                                                                                                                                                                                        >  1:145070‑M2/183‑251 (MQ=255)
                                                   cgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATc                                                                                                                                                                          >  2:98726‑M2/168‑250 (MQ=255)
                                                    gttggtttcgtgcactccgccgataatcggctcactgttgccaagcgacagcttctgcttaagcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCATAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTAGGGTAGCTATCCCAGTTGAAGATCAc                                                                                                                                                                        <  2:228537‑M2/85‑1 (MQ=255)
                                                    gttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCAc                                                                                                                                                                        >  1:11522‑M2/167‑251 (MQ=255)
                                                                 atcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGtt                                                                                                                                                            >  1:340192‑M2/154‑250 (MQ=255)
                                                                   cccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattgctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGAGTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTAcc                                                                                                                                                         <  1:205237‑M2/100‑1 (MQ=255)
                                                                                   gcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCTCCAGAAACGGGGATTAATACTCTTTACCTGTTACCGGTCTACGGTAGCTTTCCCATTTGACGATCACCCACAAACTGTTACCCAGGCGCATACGGTcc                                                                                                                                         <  1:397782‑M2/116‑1 (MQ=255)
                                                                                           gccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCAcgcg                                                                                                                                 <  2:145070‑M2/124‑1 (MQ=255)
                                                                                             caagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATc                                                                                                                                                                                       <  2:72213‑M2/70‑1 (MQ=255)
                                                                                             caagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATc                                                                                                                                                                                       >  1:72213‑M2/126‑195 (MQ=255)
                                                                                                    cagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAg                                                                                                                        <  1:143877‑M2/133‑1 (MQ=255)
                                                                                                        ttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAAGACTCTTTACCTGTTAc                                                                                                                                                                                                         <  1:155658‑M2/52‑1 (MQ=255)
                                                                                                        ttcttcttaatcagactgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTAc                                                                                                                                                                                                         >  2:155658‑M2/115‑166 (MQ=255)
                                                                                                            tcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTgg                                                                                                                >  2:330223‑M2/111‑251 (MQ=255)
                                                                                                                        tgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTcc                                                                                                     >  1:384277‑M2/99‑250 (MQ=255)
                                                                                                                         gtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCat                                                                                                   <  1:216797‑M2/154‑1 (MQ=255)
                                                                                                                              tgttaatgagaaaatcgtaaccgaatgaacgcataaagggagcgagagctcactttattattgttagcaaagtgtgcttcgttcattactgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTAc                                                                                              >  2:429481‑M2/93‑251 (MQ=255)
                                                                                                                               gttaatgagaaaatcgtatccgaatgaacgaatatagggagcaaaagctcccgttattatctctagcaaagattgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTATTCCCCAGGCGCACACGGTCCATCTCGCGTTCGCCCAGCAACTTCGTCATCTCTGCCACATTACt                                                                                             <  2:346676‑M2/160‑1 (MQ=255)
                                                                                                                                ttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTg                                                                                            >  2:378214‑M2/91‑251 (MQ=255)
                                                                                                                                                aaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgctcattcctgaTGTGTCGCC‑AGCATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTCCGGTAGCTATCCCAGTTGAAGATCACCAACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCgg                                                                            >  1:215482‑M2/75‑251 (MQ=255)
                                                                                                                                                  ccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattccttaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTTCCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACTGTCCATCACGCGTTCGCCCAGCAACTTGGTTATCGCTTCCAGATTGCTGTTGGTCAGCATCCCGGtt                                                                          <  2:215482‑M2/179‑1 (MQ=255)
                                                                                                                                                        gaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttccttcctgaTGTGTCGCC‑AGTATCCCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCATTTGCAGATCACCCACAACCTGTTACCCAGGCGCATACGGTCCCTCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGt                                                                    >  2:88410‑M2/67‑251 (MQ=255)
                                                                                                                                                             aataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTTGTCAGCATCCCGGTTGGGCGTTTGga                                                               <  1:371388‑M2/190‑1 (MQ=255)
                                                                                                                                                                                ccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATcga                                              <  1:8075‑M2/207‑1 (MQ=255)
                                                                                                                                                                                    ttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGctct                                        >  2:329498‑M2/39‑251 (MQ=255)
                                                                                                                                                                                    ttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGAt                                          <  1:244976‑M2/211‑1 (MQ=255)
                                                                                                                                                                                      attattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTgg                                      >  2:99503‑M2/37‑251 (MQ=255)
                                                                                                                                                                                      attattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTgg                                      <  1:345523‑M2/215‑1 (MQ=255)
                                                                                                                                                                                        tattgtgagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGTTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCGGGtt                                    >  2:238395‑M2/35‑251 (MQ=255)
                                                                                                                                                                                         attgttagcaaattgttcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCGCCCACAAACGGTTACCCAGGCGCATACGGTCCATCTCGCGTTCGCCCAGCAACTTGGGCATCTCTTCCATATTGCGGTTTGGCAGCATCCCGGTTGGGCGTTTTGAAGAAGAGCGGCGATCGACGATCTGGtt                                    <  2:36665‑M2/217‑1 (MQ=255)
                                                                                                                                                                                         attgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtg                                   >  1:352652‑M2/34‑251 (MQ=255)
                                                                                                                                                                                          ttgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCACCATCCCGGTTGGGCGTTTGGAAGAAGAGAGGCGTTCGCCGCTCTGGttgg                                  >  2:10288‑M2/33‑250 (MQ=255)
                                                                                                                                                                                              tagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTgatgat                              >  1:380889‑M2/29‑251 (MQ=255)
                                                                                                                                                                                               agcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTTCTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATc                             <  2:407276‑M2/224‑1 (MQ=255)
                                                                                                                                                                                                     gtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACtttt                       <  2:85052‑M2/230‑1 (MQ=255)
                                                                                                                                                                                                      tgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACttttt                      >  1:62338‑M2/21‑251 (MQ=255)
                                                                                                                                                                                                          cttcgttccttcctgaTGTGTCG‑CAAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGt                   >  2:54966‑M2/17‑250 (MQ=255)
                                                                                                                                                                                                              gttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtg                                   >  1:100970‑M2/13‑230 (MQ=255)
                                                                                                                                                                                                              gttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtg                                   <  2:100970‑M2/218‑1 (MQ=255)
                                                                                                                                                                                                                  attcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGTTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATcgac                                             <  1:68724‑M2/208‑1 (MQ=255)
                                                                                                                                                                                                                  attcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATcgac                                             >  2:68724‑M2/9‑216 (MQ=255)
                                                                                                                                                                                                                  attcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGAtt          <  1:338384‑M2/243‑1 (MQ=255)
                                                                                                                                                                                                                  attcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGAtt          <  1:187933‑M2/243‑1 (MQ=255)
                                                                                                                                                                                                                     cctgaTGTGTCGCC‑AGTATACCGCCAGCAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGCTCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCgg       >  2:169792‑M2/6‑251 (MQ=255)
                                                                                                                                                                                                                        gaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATAGTACTGTTGGTCAGCATCCCGGTTGGGCGTTGGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTc     <  1:403101‑M2/248‑1 (MQ=255)
                                                                                                                                                                                                                           gtgtCG‑CCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTGc  >  1:180439/1‑250 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
AAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCG‑C‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTGC  >  NC_000913/4599990‑4600457

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