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.C93.5% 160.4 / 4.7 46intergenic (‑25/+24)yjjA/dnaCDUF2501 domain‑containing protein YjjA/DNA replication protein DnaC
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (23/20);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

CGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCG‑C‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTGC  >  NC_000913/4599964‑4600457
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
gatgggcgtgcgtgtggtcaccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCg                                                                                                                                                                                                                                                       >  1:260889‑M2/245‑251 (MQ=255)
                   cctgttatgtgtttcgtttttgtcttctgcgtcctggcccagatcttcgcccctctggcgttgttttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAAc                                                                                                                                                                                                                                   <  2:130336‑M2/26‑1 (MQ=255)
                             gtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAAt                                                                                                                                                                                                                         <  1:544470‑M2/36‑1 (MQ=255)
                                 cggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATActc                                                                                                                                                                                                                     >  1:468859‑M2/212‑251 (MQ=255)
                                   gttgtgtcttctacgtactggacaagaccttcggccccctgacgtttgtttatggcctcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagtgagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCtt                                                                                                                                                                                                                   <  2:317754‑M2/42‑1 (MQ=255)
                                     tttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgcccgcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCttt                                                                                                                                                                                                                  >  1:204015‑M2/208‑250 (MQ=255)
                                         tcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTg                                                                                                                                                                                                             >  1:192540‑M2/204‑251 (MQ=255)
                                                 gtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTTTTACCCgg                                                                                                                                                                                                     >  2:422453‑M2/196‑251 (MQ=255)
                                                                 ccggcccgctgacgttggtttatggcatcccgccgatcagcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTa                                                                                                                                                                                         <  2:233352‑M2/68‑1 (MQ=255)
                                                                           gacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAg                                                                                                                                                                             <  2:204015‑M2/80‑1 (MQ=255)
                                                                                  gtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCcac                                                                                                                                                                    <  1:571323‑M2/89‑1 (MQ=255)
                                                                                             cccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTAcc                                                                                                                                                         <  1:410675‑M2/100‑1 (MQ=255)
                                                                                               cgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCa                                                                                                                                                       <  2:260889‑M2/102‑1 (MQ=255)
                                                                                                 ccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAgg                                                                                                                                                     >  1:524317‑M2/148‑251 (MQ=255)
                                                                                                  cgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGgc                                                                                                                                                    >  1:476296‑M2/147‑251 (MQ=255)
                                                                                                               nctgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCAt                                                                                                                                       >  2:287465‑M2/134‑251 (MQ=255)
                                                                                                                cggttgccaagcgccagcttgttcttaatcagcctgtggctgtgaatgagaaaatcgtaaccgaatgaacgaataaagggatcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTTCCCGGTTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCTCATACGGTCAATc                                                                                                                                      <  2:325858‑M2/119‑1 (MQ=255)
                                                                                                                            gccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGccc                                                                                                                          >  2:368841‑M2/121‑251 (MQ=255)
                                                                                                                                      tcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACt                                                                                                                   <  1:336655‑M2/138‑1 (MQ=255)
                                                                                                                                            tcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCAtct                                                                                                          >  2:183949‑M2/105‑251 (MQ=255)
                                                                                                                                               gcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCtt                                                                                                       >  2:373744‑M2/102‑251 (MQ=255)
                                                                                                                                                      gctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagatacctttattattgttagcaaagtgagcttcgttcattcctgaTGTGTCGCC‑AGTAGACCGCCAGCAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATAtt                                                                                                >  2:394756‑M2/95‑251 (MQ=255)
                                                                                                                                                           taatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcgtgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTg                                                                                            <  2:503316‑M2/161‑1 (MQ=255)
                                                                                                                                                               gagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGt                                                                                       <  2:524317‑M2/166‑1 (MQ=255)
                                                                                                                                                                          aaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgagcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCgg                                                                            >  1:75470‑M2/75‑251 (MQ=255)
                                                                                                                                                                           accgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGt                                                                           <  1:287465‑M2/178‑1 (MQ=255)
                                                                                                                                                                                     cgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTg                                                                 >  1:533208‑M2/64‑251 (MQ=255)
                                                                                                                                                                                           aagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGgaaga                                                            >  2:623490‑M2/58‑250 (MQ=255)
                                                                                                                                                                                                 gcgcaaactcccttcattattgttagcacagagtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGGAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCAAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCAGGTTGGGCGTTTGGAAGAAGAgcgg                                                     >  2:96368‑M2/52‑251 (MQ=255)
                                                                                                                                                                                                    aaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGgaag                                                             <  1:188409‑M2/192‑1 (MQ=255)
                                                                                                                                                                                                    aaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCAAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGgaag                                                             >  2:188409‑M2/49‑240 (MQ=255)
                                                                                                                                                                                                              ttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATc                                         <  1:623490‑M2/212‑1 (MQ=255)
                                                                                                                                                                                                                attattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAGCGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTAGCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATTCGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCTGCATCCCTGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTg                                       <  2:9598‑M2/214‑1 (MQ=255)
                                                                                                                                                                                                                attattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTgg                                      >  1:590530‑M2/37‑251 (MQ=255)
                                                                                                                                                                                                                  tattgttagcaaagtgtgctgctttcattcctgaTGTCTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGGTGAAGATCACCCACAAACTGTTACCCAGTCGCATACGGTCCATCGCGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTCCTGTTTGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGt                                     <  2:47007‑M2/216‑1 (MQ=255)
                                                                                                                                                                                                                   attgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGa                                                  <  2:506786‑M2/203‑1 (MQ=255)
                                                                                                                                                                                                                   attgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGa                                                  >  1:506786‑M2/34‑236 (MQ=255)
                                                                                                                                                                                                                   attgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtg                                   >  2:265734‑M2/34‑251 (MQ=255)
                                                                                                                                                                                                                        tagcaaagtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTgatgat                              >  2:464949‑M2/29‑251 (MQ=255)
                                                                                                                                                                                                                              agtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACttt                        <  1:183949‑M2/229‑1 (MQ=255)
                                                                                                                                                                                                                               gtgtgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATATGGTTGATGATCACtttt                       >  1:220597‑M2/22‑251 (MQ=255)
                                                                                                                                                                                                                                  tgcttcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATCCTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCg                    >  1:33062‑M2/19‑251 (MQ=255)
                                                                                                                                                                                                                                      tcgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTAtt                >  1:573890‑M2/15‑251 (MQ=255)
                                                                                                                                                                                                                                       cgttcattcctgaTGTGTCGCC‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTAttt               <  2:220597‑M2/238‑1 (MQ=255)
                                                                                                                                                                                                                                                     gtgtCG‑CCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTGc  >  2:208560/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      tgtCG‑CCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCgg       <  1:373012/244‑1 (MQ=255)
                                                                                                                                                                                                                                                      tgtCG‑CCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCgg       >  2:373012/1‑244 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
CGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCG‑C‑AGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTGC  >  NC_000913/4599964‑4600457

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