Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,396,373 C→T T103T (ACC→ACT)  amiB → N‑acetylmuramoyl‑L‑alanine amidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,396,3730CT95.8% 208.7 / ‑0.1 72T103T (ACC→ACT) amiBN‑acetylmuramoyl‑L‑alanine amidase B
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (34/34);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

CTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGC  >  NC_000913/4396139‑4396613
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
ctctGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTGCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAAc                                                                                                                                                                                                                                     <  1:479744/249‑1 (MQ=255)
ctctGATATTCAGGTTTCCAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAAc                                                                                                                                                                                                                                     <  1:131646/249‑1 (MQ=255)
  ctGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCg                                                                                                                                                                                                                                   <  1:375314/249‑1 (MQ=255)
  atGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAAcc                                                                                                                                                                                                                                    <  1:252926/247‑1 (MQ=255)
    gATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGc                                                                                                                                                                                                                               >  1:256363/1‑251 (MQ=255)
     atatTCAGGTTTCTAACGGGAATCATCAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAATATGCACAAACGCTACGGTTTGTGGTCGATCTTACTGAAAACAGTAAAACCGAAGct                                                                                                                                                                                                                              <  1:409040/251‑2 (MQ=255)
     atatTCAGGTTTCTAACGGGAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCAATTAAAGCAAACGCACCTTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTAATGCTGGTGGTCGATCTTTCTGAAAACGGTAAAACCGAAGCg                                                                                                                                                                                                                              <  1:409043/251‑1 (MQ=255)
              tttCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACgg                                                                                                                                                                                                                                           <  1:431737/229‑1 (MQ=255)
              tttCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAAAAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGAAGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACgg                                                                                                                                                                                                                                           >  2:431737/1‑229 (MQ=255)
                      ggTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATgg                                                                                                                                                                                                             >  1:499935/1‑251 (MQ=255)
                         aaTCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGAGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACACACGCTACGGCTGGGGGTCCGTCGTACTGAAAAAGGTAAAACCGaaccgaggaggcggagcaagggagc                                                                                                                                                                                                          >  2:82519/1‑228 (MQ=255)
                                aGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATAAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGAGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTCCTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAAttc                                                                                                                                                                                                    >  1:120849/1‑249 (MQ=255)
                                  gcgtggATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGGTTAGCCATCAAAGCTAACGCTCCGTGGCGTGCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGGTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTacac                                                                                                                                                                                                 <  1:440796/247‑1 (MQ=255)
                                    gcgGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTg                                                                                                                                                                                               <  1:194950/251‑1 (MQ=255)
                                     cgGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTgt                                                                                                                                                                                              <  1:260294/251‑1 (MQ=255)
                                      ggATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTgtc                                                                                                                                                                                             <  1:429939/251‑1 (MQ=255)
                                                  gTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGTGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACTCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACAGTAAAACCGAAGCGATGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGAtt                                                                                                                                                                                 <  2:450500/251‑1 (MQ=255)
                                                  gTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGCAAAAGGTAAAAACGAAACGGTGAAGCGGCAGAATGGACGCAATTACACTGTCGACTTTACGAtt                                                                                                                                                                                 >  2:456089/1‑251 (MQ=255)
                                                          tgtgATCCTGCTTATGCGGTTAGCCCTAAAAGCAAACGCACCGTGGAGCTCGATATCAAACAAACGGGCGTGATTCAGGGATGGCCGTTGTTGTGCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGCTGCACAAACGCTACGGCTGGTGGGCGATCTTACTGAAAACGGTAAACCCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACGGTCGTCTTTACGATTAACGCCg                                                                                                                                                                          <  2:50046/247‑1 (MQ=255)
                                                          ggCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCg                                                                                                                                                                          <  1:169088/250‑1 (MQ=255)
                                                            cgagccTGATTAGGCGTGTTGCCATCAAAGCAAACGCACCTTGGCGCTCGAGATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCTCAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGAt                                                                                                                                                                        <  2:503875/246‑1 (MQ=255)
                                                             gATCCTGATTATGCGTTTAGCCATCCAAGCAAACGCACCGTGGCGCTCGCTATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGACCGCCTAAAGATGCACAAACGCTACGGCTGGCGGTCGATCTTACTGAAAAAGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAAATACAATGACGTCTTTACAATTAACGCCGAtg                                                                                                                                                                       >  2:371654/1‑250 (MQ=255)
                                                                ccTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTgcc                                                                                                                                                                   <  1:382114/251‑1 (MQ=255)
                                                                 cTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGGgccg                                                                                                                                                                  >  1:467396/1‑251 (MQ=255)
                                                                  tGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTAAACTGTCGGCTTTACAGTTAACGCCGCTTTgccgc                                                                                                                                                                 >  2:99062/1‑251 (MQ=255)
                                                                     ttATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCAc                                                                                                                                                              >  2:289567/1‑251 (MQ=255)
                                                                          cGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGcct                                                                                                                                                         <  1:315231/251‑1 (MQ=255)
                                                                           gTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTAAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGACCGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTAATGAAAACGGTAAAAACGAAGCGGTGAAGCGGCCGACTGGCAGCCATTACACTGTCGTCTTTACACTTAACGACGATGTGCCGCCaatgcccc                                                                                                                                                        >  2:55533/1‑244 (MQ=255)
                                                                                     cAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCgccgcc                                                                                                                                              >  1:170664/1‑251 (MQ=255)
                                                                                                       gcgcTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAACACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGGGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGGTGCGAAACGCGt                                                                                                                           >  2:523048/1‑251 (MQ=255)
                                                                                                       gcgcTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGCATAACGCCGATGTGCCGCCACCGCCCCCTCCGCCGCCC‑GGGGGTGCGAAACGCGt                                                                                                                           >  2:481024/1‑251 (MQ=255)
                                                                                                       gcgcTCGAAATCAAAAAAACGGGCGGGATTCAGGGCCTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAAAGCTACGGCTGGTGGTCGATCTTAATGAAAACGGTACAAACGCCGCGGGGCAGCGGCAGAATGGCCGCCATTTCACTGACGTCTTTACGGTTAACGCCGATGTGCCGCCACCGCCTCCTCAGCCGCCCCGTG‑TTGCGAACcgctt                                                                                                                           >  2:523060/1‑249 (MQ=255)
                                                                                                          cTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTg                                                                                                                         >  1:298535/1‑250 (MQ=255)
                                                                                                            cGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAACCCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGaaa                                                                                                                      <  1:433442/251‑1 (MQ=255)
                                                                                                              atatCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGtt                                                                                                                                      >  2:183987/1‑233 (MQ=255)
                                                                                                              atatCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGtt                                                                                                                                      <  1:183987/233‑1 (MQ=255)
                                                                                                                  caaacaaaCGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTACGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTg                                                                                                                <  1:113719/251‑1 (MQ=255)
                                                                                                                      caaacGGGAGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGCACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAACCGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGAAGCAATTACCCTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGCAACGCGTTGCACAGCCTGCggg                                                                                                            >  2:311838/1‑250 (MQ=255)
                                                                                                                       aaacGGGCGGGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCGGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGCCTTTACGATTAACGCAGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGtt                                                                                                           <  1:347580/251‑1 (MQ=255)
                                                                                                                                          aCTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAACCGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAACCCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTTCGATTAAAGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCCAAAAGCGTTGCAACGCCTGCGGTTTTCGCCCCGCGCGttcgcc                                                                                        >  2:503659/1‑246 (MQ=255)
                                                                                                                                              ccGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAAcc                                                                                    <  1:8326/251‑1 (MQ=255)
                                                                                                                                                gttgttgttCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCgg                                                                                  >  1:438514/1‑251 (MQ=255)
                                                                                                                                                  tgttgttCAGCGGCAATAATCTGGTGAAGGCGAGTCGCTCTGGAACGCCTAAAGTTGCACAAACGCTGCGGGTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTTTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGGTGAAACGCCTGCGGTTGTCGCACCGCGCTTCAGCGAACCGgc                                                                                 <  2:234040/250‑1 (MQ=255)
                                                                                                                                                     tgttAAGCGGCAATAATCTGGTGAAGGCTATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACCCTGTCGGCTTTACGATTAACGCCGATGGGCCGCCACCGCCGCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGGCGCAGCGCGCGTCAGCGAACCGgcgcg                                                                              <  2:345454/250‑1 (MQ=255)
                                                                                                                                                         cAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAAt                                                                          <  1:99062/250‑1 (MQ=255)
                                                                                                                                                           gCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGCAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGCACCGGCGCGCAATCCg                                                                       >  1:229287/1‑251 (MQ=255)
                                                                                                                                                              gCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGAGGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCCGCGGTTGTCGCACCCCGCGTCAGCGCACCGGCGCGCaaacccgtt                                                                    >  2:182338/1‑244 (MQ=255)
                                                                                                                                                               cAATAATCTGGTGAAGGCGTTTCGCTCTGGCACGCCTAAAGATGCCCAAACGCTACGGCTGTTGGTCGATCTTACTGAAAACGTTAATACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTa                                                                   <  1:2266/251‑1 (MQ=255)
                                                                                                                                                                aataatCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTa                                                                   >  1:4228/1‑250 (MQ=255)
                                                                                                                                                                 acaatgtggggAAGGCTATTCGCTCGGGATCGCTTAAAGTTGCACAAATGTGAGGGCGGGTGGTCGATCTTACGGAAAACGGTAAAACCGAAGCGGTGAAGCGGCGGAAGGGCCGCAATTACACTGTCGCCTTTAAGATTATCGCCGATGTGCCGCATCCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGAGGTTGGCGCACCGCGCGTCAGCGATCCGGCGCGCAATCCGTTTaaa                                                                 <  2:310400/241‑1 (MQ=255)
                                                                                                                                                                         gTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGaa                                                           >  1:203468/1‑249 (MQ=255)
                                                                                                                                                                           gAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAACCGGAAAGt                                                        >  1:287550/1‑250 (MQ=255)
                                                                                                                                                                               gCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGTTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCg                                                   <  2:438662/251‑1 (MQ=255)
                                                                                                                                                                                   ttCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTAAAAACGGAAAGTCACCGCaac                                               >  2:502064/1‑249 (MQ=255)
                                                                                                                                                                                         ctGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAATCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGAtgt                                                                                                                                                                      <  2:391378/127‑1 (MQ=255)
                                                                                                                                                                                         ctGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAATCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGAtgt                                                                                                                                                                      >  1:391378/1‑127 (MQ=255)
                                                                                                                                                                                         ctGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCTCCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACgg                                           <  1:75796/249‑1 (MQ=255)
                                                                                                                                                                                           ggAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGgt                                         <  1:491783/249‑1 (MQ=255)
                                                                                                                                                                                              aCGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACCCTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTAt                                    <  1:302327/251‑1 (MQ=255)
                                                                                                                                                                                                    aaaGATGCACAACCGCTCCGGCTGGTGGTCGCTCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCCGCAATTAAACTGACGTCTTTACGATTAACGCCGATGTGCCGCCACCACCTCCTCCGCAGCCC‑GTGGTTGCGAAACGCGTTGAAAAGCCTGCGGTTGACGCCCCGCGCGTCCGCGAAACGGCGCGCAATCCGTTTAAAACGGGAAGAAACCGCACCACAGGTGTTATcaaccg                              >  2:518840/1‑247 (MQ=255)
                                                                                                                                                                                                    aaaGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACACGTTGAAACGCCTGCGGTTGACGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATcagcag                              <  1:315318/251‑1 (MQ=255)
                                                                                                                                                                                                     aaGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACCCTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCACCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCCACTCGCGTCCGCCAACCCGCGCGCCATCCGTTTAAAACGGAAAGTAACCGCACTACGGGGGTTTTcaacaat                             >  2:490719/1‑246 (MQ=255)
                                                                                                                                                                                                       gATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGGTTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAACAGCCTGAGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTa                            >  2:187582/1‑250 (MQ=255)
                                                                                                                                                                                                            acaAACGCTACGGCTGGTGGTATATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCACTTCAACTGCCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGAACTACGGGTGTTATCAGCAGTAATACgg                      <  1:82519/251‑1 (MQ=255)
                                                                                                                                                                                                                actcTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCTGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAAc                  <  1:452074/248‑1 (MQ=255)
                                                                                                                                                                                                                aCGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTGACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCTCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAAc                  <  1:55533/251‑1 (MQ=255)
                                                                                                                                                                                                                  gCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGACACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCCGTAATACGGTAAccc                >  2:64465/1‑249 (MQ=255)
                                                                                                                                                                                                                      cGGCTGGTGGTCGATCTTACTGCAAACGGTAAAAACGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTAAGATTAACGACGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGACACGCCTGCGGGTGTCGCAACGCCCGTCAGCGAACCGCAGCGCACTCCGGTTAAAACGGAAAGGAACCGCACACCGGAGGTTATCAGCAGTAAAACGGTAAcgccacc            >  2:205162/1‑247 (MQ=255)
                                                                                                                                                                                                                           ggtggtCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGc         <  1:289567/249‑1 (MQ=255)
                                                                                                                                                                                                                                tCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACCCTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGAGCGCCAGCGAAACGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGGGTTATCAGCAGTAATACGGTAACGCGTCCGGCAgcgcgc  >  1:538037/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cTTACTGAAAACGGTAAAACCGAAGGGGTGAAGCGGCAGAATGGCAGCAATTACCCTTTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGt                     <  2:383844/227‑1 (MQ=255)
                                                                                                                                                                                                                                     cTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGGGTTATCAGCAGTAATACGGt                     >  1:383844/1‑227 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
CTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCC‑GTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGC  >  NC_000913/4396139‑4396613

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