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→ACTamiBN‑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

CTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGC  >  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)
                                                                                                       gcgcTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAACACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGGGCCGCCACCGCCTCCTCCGCCGCCCGTGGGTGCGAAACGCGt                                                                                                                           >  2:523048/1‑251 (MQ=255)
                                                                                                       gcgcTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGCATAACGCCGATGTGCCGCCACCGCCCCCTCCGCCGCCCGGGGGTGCGAAACGCGt                                                                                                                           >  2:481024/1‑251 (MQ=255)
                                                                                                       gcgcTCGAAATCAAAAAAACGGGCGGGATTCAGGGCCTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAAAGCTACGGCTGGTGGTCGATCTTAATGAAAACGGTACAAACGCCGCGGGGCAGCGGCAGAATGGCCGCCATTTCACTGACGTCTTTACGGTTAACGCCGATGTGCCGCCACCGCCTCCTCAGCCGCCCCGTGTTGCGAACcgctt                                                                                                                           >  2:523060/1‑249 (MQ=255)
                                                                                                          cTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTg                                                                                                                         >  1:298535/1‑250 (MQ=255)
                                                                                                            cGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAACCCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGaaa                                                                                                                      <  1:433442/251‑1 (MQ=255)
                                                                                                              atatCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGtt                                                                                                                                      >  2:183987/1‑233 (MQ=255)
                                                                                                              atatCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGtt                                                                                                                                      <  1:183987/233‑1 (MQ=255)
                                                                                                                  caaacaaaCGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTACGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTg                                                                                                                <  1:113719/251‑1 (MQ=255)
                                                                                                                      caaacGGGAGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGCACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAACCGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGAAGCAATTACCCTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGCAACGCGTTGCACAGCCTGCggg                                                                                                            >  2:311838/1‑250 (MQ=255)
                                                                                                                       aaacGGGCGGGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCGGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGCCTTTACGATTAACGCAGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGtt                                                                                                           <  1:347580/251‑1 (MQ=255)
                                                                                                                                          aCTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAACCGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAACCCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTTCGATTAAAGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCCAAAAGCGTTGCAACGCCTGCGGTTTTCGCCCCGCGCGttcgcc                                                                                        >  2:503659/1‑246 (MQ=255)
                                                                                                                                              ccGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAAcc                                                                                    <  1:8326/251‑1 (MQ=255)
                                                                                                                                                gttgttgttCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCgg                                                                                  >  1:438514/1‑251 (MQ=255)
                                                                                                                                                  tgttgttCAGCGGCAATAATCTGGTGAAGGCGAGTCGCTCTGGAACGCCTAAAGTTGCACAAACGCTGCGGGTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTTTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGGTGAAACGCCTGCGGTTGTCGCACCGCGCTTCAGCGAACCGgc                                                                                 <  2:234040/250‑1 (MQ=255)
                                                                                                                                                     tgttAAGCGGCAATAATCTGGTGAAGGCTATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACCCTGTCGGCTTTACGATTAACGCCGATGGGCCGCCACCGCCGCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGGCGCAGCGCGCGTCAGCGAACCGgcgcg                                                                              <  2:345454/250‑1 (MQ=255)
                                                                                                                                                         cAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAAt                                                                          <  1:99062/250‑1 (MQ=255)
                                                                                                                                                           gCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGCAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGCACCGGCGCGCAATCCg                                                                       >  1:229287/1‑251 (MQ=255)
                                                                                                                                                              gCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGAGGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCCGCGGTTGTCGCACCCCGCGTCAGCGCACCGGCGCGCaaacccgtt                                                                    >  2:182338/1‑244 (MQ=255)
                                                                                                                                                               cAATAATCTGGTGAAGGCGTTTCGCTCTGGCACGCCTAAAGATGCCCAAACGCTACGGCTGTTGGTCGATCTTACTGAAAACGTTAATACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTa                                                                   <  1:2266/251‑1 (MQ=255)
                                                                                                                                                                aataatCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTa                                                                   >  1:4228/1‑250 (MQ=255)
                                                                                                                                                                 acaatgtggggAAGGCTATTCGCTCGGGATCGCTTAAAGTTGCACAAATGTGAGGGCGGGTGGTCGATCTTACGGAAAACGGTAAAACCGAAGCGGTGAAGCGGCGGAAGGGCCGCAATTACACTGTCGCCTTTAAGATTATCGCCGATGTGCCGCATCCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGAGGTTGGCGCACCGCGCGTCAGCGATCCGGCGCGCAATCCGTTTaaa                                                                 <  2:310400/241‑1 (MQ=255)
                                                                                                                                                                         gTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGaa                                                           >  1:203468/1‑249 (MQ=255)
                                                                                                                                                                           gAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAACCGGAAAGt                                                        >  1:287550/1‑250 (MQ=255)
                                                                                                                                                                               gCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGTTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCg                                                   <  2:438662/251‑1 (MQ=255)
                                                                                                                                                                                   ttCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTAAAAACGGAAAGTCACCGCaac                                               >  2:502064/1‑249 (MQ=255)
                                                                                                                                                                                         ctGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAATCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGAtgt                                                                                                                                                                      <  2:391378/127‑1 (MQ=255)
                                                                                                                                                                                         ctGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAATCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGAtgt                                                                                                                                                                      >  1:391378/1‑127 (MQ=255)
                                                                                                                                                                                         ctGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCTCCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACgg                                           <  1:75796/249‑1 (MQ=255)
                                                                                                                                                                                           ggAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGgt                                         <  1:491783/249‑1 (MQ=255)
                                                                                                                                                                                              aCGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACCCTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTAt                                    <  1:302327/251‑1 (MQ=255)
                                                                                                                                                                                                    aaaGATGCACAACCGCTCCGGCTGGTGGTCGCTCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCCGCAATTAAACTGACGTCTTTACGATTAACGCCGATGTGCCGCCACCACCTCCTCCGCAGCCCGTGGTTGCGAAACGCGTTGAAAAGCCTGCGGTTGACGCCCCGCGCGTCCGCGAAACGGCGCGCAATCCGTTTAAAACGGGAAGAAACCGCACCACAGGTGTTATcaaccg                              >  2:518840/1‑247 (MQ=255)
                                                                                                                                                                                                    aaaGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACACGTTGAAACGCCTGCGGTTGACGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATcagcag                              <  1:315318/251‑1 (MQ=255)
                                                                                                                                                                                                     aaGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACCCTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCACCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCCACTCGCGTCCGCCAACCCGCGCGCCATCCGTTTAAAACGGAAAGTAACCGCACTACGGGGGTTTTcaacaat                             >  2:490719/1‑246 (MQ=255)
                                                                                                                                                                                                       gATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGGTTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAACAGCCTGAGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTa                            >  2:187582/1‑250 (MQ=255)
                                                                                                                                                                                                            acaAACGCTACGGCTGGTGGTATATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCACTTCAACTGCCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGAACTACGGGTGTTATCAGCAGTAATACgg                      <  1:82519/251‑1 (MQ=255)
                                                                                                                                                                                                                actcTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCTGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAAc                  <  1:452074/248‑1 (MQ=255)
                                                                                                                                                                                                                aCGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTGACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCTCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAAc                  <  1:55533/251‑1 (MQ=255)
                                                                                                                                                                                                                  gCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGACACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCCGTAATACGGTAAccc                >  2:64465/1‑249 (MQ=255)
                                                                                                                                                                                                                      cGGCTGGTGGTCGATCTTACTGCAAACGGTAAAAACGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTAAGATTAACGACGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGACACGCCTGCGGGTGTCGCAACGCCCGTCAGCGAACCGCAGCGCACTCCGGTTAAAACGGAAAGGAACCGCACACCGGAGGTTATCAGCAGTAAAACGGTAAcgccacc            >  2:205162/1‑247 (MQ=255)
                                                                                                                                                                                                                           ggtggtCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGc         <  1:289567/249‑1 (MQ=255)
                                                                                                                                                                                                                                tCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACCCTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGAGCGCCAGCGAAACGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGGGTTATCAGCAGTAATACGGTAACGCGTCCGGCAgcgcgc  >  1:538037/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cTTACTGAAAACGGTAAAACCGAAGGGGTGAAGCGGCAGAATGGCAGCAATTACCCTTTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGt                     <  2:383844/227‑1 (MQ=255)
                                                                                                                                                                                                                                     cTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGGGTTATCAGCAGTAATACGGt                     >  1:383844/1‑227 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 
CTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGC  >  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: