Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,062,357 T→A 7.6% intergenic (+22/+34) asnV → / ← ldtA tRNA‑Asn/L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,062,3570TA 7.6% 158.7 / 3.6 53intergenic (+22/+34)asnV/ldtAtRNA‑Asn/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base T (27/22);  new base A (2/2);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

CATTTGCTCCTTTAATCAGCTGTCGCGTTCCCCTGCCCTATAAAAGGAGGGTATGCACCACGATGGTTCATTACCCAATAAGATTGAAAGCTCACCACTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCGAGCCATCTGGTTCAGTGGTGTATTTTACTGGCTGGTCGATAATTTGCACACGC  >  NC_000913/2062069‑2062646
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                  
cATTTGCTCCTTGAATCTGCTGTCGCGTGCCCCTGCCCTATAAAAGGATGGTATGCACTACGATGGTTCATTACCCAATAAGATTGATAGCTCACCACTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATGTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCttttt                                                                                                                                                                                                                                                                                       <  1:340276/301‑1 (MQ=37)
      ttCCTTTAATCAGCTGTCGCGTTCCCCTGCCCTATAAAAGGAGGGTATGCACCACGATGGTTCATTACCCAATAAGATTGAAAGCTCACCACTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTc                                                                                                                                                                                                                                                                                 <  2:639799/300‑1 (MQ=255)
                             ccccTGCCCTATAAAAGGAGGGTATGCACCACGATGGTTCATTACCCAATAAGATTGAAAGCTCACCACTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTACTAGCCGGCTTTTTGCTTTCTAATTACCAACGCACTTAAAAcc                                                                                                                                                                                                                                                          >  1:408914/1‑300 (MQ=255)
                                        taAAAGGAGGGTATGCACCACGATGGTTCATTACCCAATAATATTGAAAGCTCACCACTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGCGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTACAACATCTGTCtttac                                                                                                                                                                                                                                               >  1:504707/1‑298 (MQ=255)
                                                                      ttACCCAATAAGATTGAAAGCTCACCACTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGaa                                                                                                                                                                                                                 >  2:535637/1‑300 (MQ=255)
                                                                                    tGAAAGCTCGCCACTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAACAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgc                                                                                                                                                                                                   >  1:231275/1‑301 (MQ=255)
                                                                                     gAAAGCTCACCACTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGATTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTACCGATCGACTTTGCAACgca                                                                                                                                                                                                  >  2:558287/1‑301 (MQ=255)
                                                                                      aaaGCTCACCACTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgcag                                                                                                                                                                                                 <  2:415317/301‑1 (MQ=255)
                                                                                       aaGCTCACCACTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCCGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACCACTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGAACAACGTTGCAACccagc                                                                                                                                                                                                >  1:624314/1‑301 (MQ=255)
                                                                                       aaGCTCACCACTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGAATTCCCGACCGACGTACCAACgcaac                                                                                                                                                                                                >  1:66790/1‑299 (MQ=255)
                                                                                            caccacTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGCCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGAAGCAtttt                                                                                                                                                                                           >  2:348697/1‑300 (MQ=255)
                                                                                            aacctcTTTGTTGAAATTGATAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTg                                                                                                                                                                                           <  1:190150/300‑1 (MQ=255)
                                                                                                       tgAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAACAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACa                                                                                                                                                                                >  1:562977/1‑301 (MQ=255)
                                                                                                                      acaagcaaaAATATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTGATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCg                                                                                                                                                                 <  1:159026/301‑1 (MQ=255)
                                                                                                                           caaaAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGAt                                                                                                                                                            >  2:18638/1‑301 (MQ=255)
                                                                                                                             aaaaaaTGCTTTTCGCCCTTTGACATCACCATGCACGGCCATTAATATGTGCCCCGTTCACACGATTCCGCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATaa                                                                                                                                                          <  1:601213/301‑1 (MQ=255)
                                                                                                                                     cATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGc                                                                                                                                                  <  2:231518/301‑1 (MQ=255)
                                                                                                                                       tttccgatcttGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCaa                                                                                                                                                <  2:223083/292‑1 (MQ=255)
                                                                                                                                             tcTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATgg                                                                                                                                          <  2:625173/300‑1 (MQ=255)
                                                                                                                                               tttGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATgggg                                                                                                                                        <  2:505379/301‑1 (MQ=255)
                                                                                                                                                ttGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGt                                                                                                                                       <  2:563691/301‑1 (MQ=255)
                                                                                                                                                 tGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAAATcggaaa                                                                                                                                      >  1:222840/1‑298 (MQ=255)
                                                                                                                                                  gACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTa                                                                                                                                     >  2:623527/1‑301 (MQ=255)
                                                                                                                                                      tCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAAGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCgg                                                                                                                                 >  2:108566/1‑301 (MQ=255)
                                                                                                                                                        aCCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCa                                                                                                                               >  1:344459/1‑301 (MQ=255)
                                                                                                                                                            tggagttcagacgtgtgctcttccGATCTCACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATgg                                                                                                                           <  2:546301/279‑1 (MQ=255)
                                                                                                                                                                    cATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCACCTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGATACCGACAATGGCACttatt                                                                                                                   >  1:72519/1‑298 (MQ=255)
                                                                                                                                                                            tGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGAc                                                                                                           >  2:253066/1‑301 (MQ=255)
                                                                                                                                                                            tGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCACGCAAAGATGGGGTTACCGGCAATGGCCCTTTTCGGTCAGAc                                                                                                           >  1:489037/1‑301 (MQ=255)
                                                                                                                                                                               gcCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCa                                                                                                        >  2:6503/1‑301 (MQ=255)
                                                                                                                                                                                     ttCACCCGTTTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCt                                                                                                  <  1:252798/301‑1 (MQ=255)
                                                                                                                                                                                       ctcacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGc                                                                                                <  2:409477/299‑1 (MQ=255)
                                                                                                                                                                                         cacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATggagatcggaagagcgtcgtgtagggaaaga                                                                                              >  1:545628/1‑272 (MQ=255)
                                                                                                                                                                                         cacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCAc                                                                                              >  2:238969/1‑301 (MQ=255)
                                                                                                                                                                                            taTTTCTCTGTAGTTCAGTCGTTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGAt                                                                                           <  2:510313/300‑1 (MQ=255)
                                                                                                                                                                                                        cttCAGTCGGTAGAACGGCGGACTTTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATg                                                                               <  1:79609/300‑1 (MQ=255)
                                                                                                                                                                                                          tCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATgaa                                                                             >  2:79699/1‑299 (MQ=255)
                                                                                                                                                                                                                ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGa                                                                       <  2:344922/301‑1 (MQ=255)
                                                                                                                                                                                                                       cgggggACTGTTATTTCGTTGTTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCGTTCTAATTACCAACGCTCTGAAAACATCGGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTcca                                                                <  1:534964/301‑1 (MQ=255)
                                                                                                                                                                                                                                ttttaTCCTTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCg                                                       <  2:353687/297‑1 (MQ=255)
                                                                                                                                                                                                                                gTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAACACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGAATCATATTCTGCACGATTGAGCAACAATGGCCAATGAAATTccacccatctcc                                                       >  1:353224/1‑296 (MQ=255)
                                                                                                                                                                                                                                        gtgtGTCACTGGTTCGAGTCCAGTCAGAGGGGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTGTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCGAGCCATCt                                               <  1:6503/298‑1 (MQ=255)
                                                                                                                                                                                                                                                        agcccatgtagtggtgCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAAGTACCAACGCTCTTAAAACAGCTGTCTTGAACCAGAACTAAGTGGCACAGGCAGTCCCGATCGACGGTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCGGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCGAGCCATCTGGTTCAGTGGTGTAtt                               <  1:348234/286‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTCCAGTCAGAGGAGCCAAATTCAACAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTGTTCGGTCAGCCTCCTATTCTGCACGATTGCGAGCCAATGCCTCATGAACTTCCAGCCATTTCGAGCCATCTGGTTCAGTGGTGTAttt                              >  1:355944/1‑301 (MQ=255)
                                                                                                                                                                                                                                                           ccAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCCGACTCATATTCTGCACGATTGCGCGACAATGACTCATGAACTTCCAGCCAGTTCACGCCATATGGTTCAGTGGTGTAttttt                            >  1:390989/1‑300 (MQ=255)
                                                                                                                                                                                                                                                            cAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCGAGCCATCTGGTTCAGTGGTGTATTTTAc                           <  2:489683/301‑1 (MQ=255)
                                                                                                                                                                                                                                                            cAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCGAGCCATCTGGTTCAGTGGTGTATTTTAc                           <  2:66880/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACCCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCGAGCCATCTGGTTCAGTGGTGTATTTTACTGGCTGGTCGATa              >  2:271763/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                              ccgttcttgCCTGCTTTCTATCAGGCTTTTTGCTTTCTAATTCCCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCGAGCCATCTGGTTCAGTGGTGTATTTTACTGGCTGGTCGATAATTTg         <  1:308508/293‑1 (MQ=255)
                                                                                                                                                                                                                                                                                caaaaaGCCTGCTAGAAAGCAGGCTTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCTTTACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCGAGCCATCTGGTTCAGTGGTGTATTTTACTGGCTGGTCGATAAATTg       >  2:251382/2‑300 (MQ=255)
                                                                                                                                                                                                                                                                                caaaaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCGAGCCATCTGGTTCAGTGGTGTATTTTACTGGCTGGTCGATAATTTGca       >  2:151996/2‑301 (MQ=255)
                                                                                                                                                                                                                                                                                aaaaaaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCGAGCCATCTGGTTCAGTGGTGTATTTTACTGGCTGGTCGATAATTAGca       >  2:29023/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                     aGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCGAGCCATCTGGTTCAGTGGTGTATTTTACTGGCTGGTCGAAAATTTgagaacgg  >  2:308912/1‑294 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                  
CATTTGCTCCTTTAATCAGCTGTCGCGTTCCCCTGCCCTATAAAAGGAGGGTATGCACCACGATGGTTCATTACCCAATAAGATTGAAAGCTCACCACTTTGTTGAAATTGACAGCAAACAAACAAAAAAATGCATTTCACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAAAGATGGGGTTACCGGCAATGGCACTTTTCGGTCAGACTCATATTCTGCACGATTGCGCGACAATGGCTCATGAACTTCCAGCCAGTTCGAGCCATCTGGTTCAGTGGTGTATTTTACTGGCTGGTCGATAATTTGCACACGC  >  NC_000913/2062069‑2062646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: