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.C98.2% 203.2 / NA 57intergenic (‑25/+24)yjjA/dnaCDUF2501 domain‑containing protein YjjA/DNA replication protein DnaC
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (29/27);  minor base A (1/0);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTG  >  NC_000913/4599974‑4600456
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
cgtgtggtcaccggtatcagtttcggttttgtcttctacggactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagctaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATAcc                                                                                                                                                                                                                                           <  1:549986‑M2/17‑1 (MQ=255)
      gtcaccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGcc                                                                                                                                                                                                                                        >  1:320932‑M2/229‑248 (MQ=255)
         cccggtatcagtttcggttttgtcttctacgtactggaccagatcgtcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAAc                                                                                                                                                                                                                                  <  1:250223‑M2/26‑1 (MQ=255)
           cggtatcggtttcggttttttcttctccgttctggaccagatcttcggcccgctgacgttgtttaatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctcccttgattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACgg                                                                                                                                                                                                                                <  2:320040‑M2/28‑1 (MQ=255)
             gtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGACACgggg                                                                                                                                                                                                                              >  1:211147‑M2/222‑251 (MQ=255)
               atcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgtgggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattggtagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGAt                                                                                                                                                                                                                            <  1:342341‑M2/32‑1 (MQ=255)
               agcagtttcggttttgtcttctacgtactggaccagatcgtcggcccgctgacgttggtgtatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgtgcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGAt                                                                                                                                                                                                                            <  1:341350‑M2/32‑1 (MQ=255)
                  agtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatccgcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcatttctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTa                                                                                                                                                                                                                          >  2:455143‑M2/217‑250 (MQ=255)
                   gtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAAt                                                                                                                                                                                                                        <  1:223559‑M2/36‑1 (MQ=255)
                               tcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTg                                                                                                                                                                                                            >  1:149770‑M2/204‑251 (MQ=255)
                                   ctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTAc                                                                                                                                                                                                        >  2:123235‑M2/200‑251 (MQ=255)
                                             gaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattcttgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCCGAAACGGGGCTTAATACTCTTTACCTGTTACCCGGCTACgg                                                                                                                                                                                              >  1:353550‑M2/190‑251 (MQ=255)
                                             gaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgacaagcgccagcttcttcttaatcagcctgtggctgttaatgagcaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACgg                                                                                                                                                                                              >  1:378211‑M2/190‑251 (MQ=255)
                                                    tcttcggcccgctgactttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTa                                                                                                                                                                                        <  2:343656‑M2/68‑1 (MQ=255)
                                                     cttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATc                                                                                                                                                                                      >  1:415806‑M2/182‑251 (MQ=255)
                                                       tctgccagctgacgttgttttatggcagcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgtgaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGTGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATc                                                                                                                                                                                      <  1:402283‑M2/70‑1 (MQ=255)
                                                         ggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTAACTGTTACCCGGCTACGGTAGCTATccc                                                                                                                                                                                    >  1:537559‑M2/178‑249 (MQ=255)
                                                              gctgacgttggtttatggcataccgccgatcatcggcgcactgttgccaagcgacagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagaaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGa                                                                                                                                                                              >  1:540730‑M2/173‑250 (MQ=255)
                                                                 gacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGAt                                                                                                                                                                          <  1:539089‑M2/82‑1 (MQ=255)
                                                                      tggtttatggcaccccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaaggagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCAcc                                                                                                                                                                      <  2:529244‑M2/86‑1 (MQ=255)
                                                                            atggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAAc                                                                                                                                                               <  2:428898‑M2/93‑1 (MQ=255)
                                                                              ggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGTTAGCTATCCCAGTTGAAGATCACCCACACACTg                                                                                                                                                             >  1:517662‑M2/157‑251 (MQ=255)
                                                                                    ccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTAccc                                                                                                                                                       <  2:415806‑M2/101‑1 (MQ=255)
                                                                                             atcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATa                                                                                                                                              <  2:178508‑M2/110‑1 (MQ=255)
                                                                                               cggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTGCCTGTTACCCGGCTACGGTAGCTATCCCAGTGGAAGATCACCCACAATCTGTTACCCAGGCGCATACg                                                                                                                                            <  2:550922‑M2/112‑1 (MQ=255)
                                                                                                ggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGt                                                                                                                                                                                                           <  1:153956‑M2/49‑1 (MQ=255)
                                                                                                ggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGt                                                                                                                                                                                                           >  2:153956‑M2/139‑187 (MQ=255)
                                                                                                              aatcgccagcttcttattaagcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgattaaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTCCCTGTTACCCGGCTACGGTAGCTGTCCCAGTTGAAGATCACCCACAAACGGTTTCCCAGGCGCATACGGTACATCACGCGTTc                                                                                                                             <  2:253815‑M2/127‑1 (MQ=255)
                                                                                                              aagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTc                                                                                                                             >  1:103432‑M2/125‑251 (MQ=255)
                                                                                                                               taatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGt                                                                                                              <  1:168735‑M2/142‑1 (MQ=255)
                                                                                                                                aatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCCACTTGgtcct                                                                                                           >  2:475113‑M2/107‑249 (MQ=255)
                                                                                                                                   cagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCAtct                                                                                                         <  1:440107‑M2/147‑1 (MQ=255)
                                                                                                                                             ctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATAtt                                                                                               >  2:421710‑M2/94‑250 (MQ=255)
                                                                                                                                                  aatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGtt                                                                                         >  2:166950‑M2/89‑251 (MQ=255)
                                                                                                                                                         aaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGc                                                                                  >  1:154409‑M2/82‑251 (MQ=255)
                                                                                                                                                          aatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGc                                                                                  >  2:428318‑M2/81‑250 (MQ=255)
                                                                                                                                                              gtaaccgaatgaacgaataaagggagctaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATccc                                                                             <  2:434528‑M2/175‑1 (MQ=255)
                                                                                                                                                                                 aagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGgaaga                                                           >  2:295209‑M2/58‑250 (MQ=255)
                                                                                                                                                                                   gggagcgaaagctccctttattattgttagcaaagtgtgcttcggtcattcctgaTGTGTCGCCAGTATACCGCCATAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTAGCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTGCGCCCAGCAACTTGGTCAGCTCGTCCATATGACTTTTGGGCATCATCCCGTTTGGGCGTTTGGAagaaga                                                        <  2:540730‑M2/196‑1 (MQ=255)
                                                                                                                                                                                     gagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCCTATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgc                                                      >  1:567257‑M2/54‑251 (MQ=255)
                                                                                                                                                                                          aaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcggcg                                                  >  1:558879‑M2/49‑250 (MQ=255)
                                                                                                                                                                                          aaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCAAGACACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCACAGTTGGAGATCACCAACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCCTCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGCAGAgcggcg                                                  >  2:88621‑M2/49‑250 (MQ=255)
                                                                                                                                                                                          aaagctccctttattattgttagcaaaaagtgcttcgcttattcctgaTGTGTCGACAGTATACCGCCAGAAACGGGGATTAATACTCTTTCACTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACACACAAACTGTTACCCAGGCGCATACGGTCCATCACGCATTCGCCCAGCAAATTGGTCATCTCTTCCATATTAATTTTGGTCAGCATCCAGGTTGGGCGTTTGGAAGAAGtgcggcgt                                                 >  2:53228‑M2/49‑250 (MQ=255)
                                                                                                                                                                                           aagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGAt                                                >  1:461529‑M2/48‑251 (MQ=255)
                                                                                                                                                                                             gctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATcg                                              <  1:492688‑M2/206‑1 (MQ=255)
                                                                                                                                                                                                    ttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGAt                                         <  1:198737‑M2/211‑1 (MQ=255)
                                                                                                                                                                                                    ttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCt                                       >  1:353977‑M2/39‑251 (MQ=255)
                                                                                                                                                                                                            gttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtgat                                <  2:149770‑M2/220‑1 (MQ=255)
                                                                                                                                                                                                             ttagcaaagtgtgcttcgttcattccttaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCGTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtgatg                               <  2:456629‑M2/221‑1 (MQ=255)
                                                                                                                                                                                                                 caaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATGCCGCCAGAAACGGGGATTAATTCTCTTTACCTGTTTCCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCGTCCATATTACTGTTGGTCAGCATCCCGGTTGGTCGGTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCAc                          <  1:405140‑M2/226‑1 (MQ=255)
                                                                                                                                                                                                                     gtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACtttt                      >  1:172487‑M2/22‑251 (MQ=255)
                                                                                                                                                                                                                                 cattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGAt          >  1:94970‑M2/10‑251 (MQ=255)
                                                                                                                                                                                                                                      ctgaTGTGTCGCCAGTATACCGCCAGCAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATc                                        >  2:92252‑M2/5‑216 (MQ=255)
                                                                                                                                                                                                                                      ctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCcaca                                                                                                                                                                  >  2:288315‑M2/5‑94 (MQ=255)
                                                                                                                                                                                                                                      ctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATc                                        <  1:92252‑M2/212‑1 (MQ=255)
                                                                                                                                                                                                                                      ctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACGCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCcaca                                                                                                                                                                  <  1:288315‑M2/90‑1 (MQ=255)
                                                                                                                                                                                                                                         aTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTctc  >  1:183401‑M2/2‑250 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
CATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTG  >  NC_000913/4599974‑4600456

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: