Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,003,508 T→G 21.6% T224P (ACC→CCC)  rarD ← putative chloramphenical resistance permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,003,5080TG21.6% 102.6 / 2.2 60T224P (ACC→CCC) rarDputative chloramphenical resistance permease
Reads supporting (aligned to +/- strand):  ref base T (17/26);  new base G (14/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAA  >  NC_000913/4003275‑4003738
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
tatatTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTg                                                                                                                                                                                                                                        >  1:512100/1‑238 (MQ=255)
             tatttGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAAAAGCGGGACGGGAAGGAGAATACGggcggc                                                                                                                                                                                                                          >  2:722052/1‑239 (MQ=255)
                        tCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGGAAGGTAGTGACAATACCGGCGGCGAACAGCAGTa                                                                                                                                                                                                               >  1:22920/1‑239 (MQ=255)
                                gagTAGAAATCGCAGCCATCACAAAAAGGGCCACCGCGACCCAAATAAAGGCGAAAGTCCCCATCTGATCTGCACCCGGTTTTGCACCATAAAACGTCACAGCCAGCAGGAACAGCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCGGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATGCCGGCGGCGATCAGCAGTAAATTCAg                                                                                                                                                                                                        <  2:221508/236‑1 (MQ=255)
                                 cGTATACATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATACAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAACAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGGAAAAAACAACAGCGGTACGGGAGAGACAAAAGCCGCGGCCGGCAGGAGTAAAATCAACg                                                                                                                                                                                                      >  2:902853/1‑239 (MQ=255)
                                 cGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGCGACGCAAGCGCGGGGCAGCGGCGGCAAAACACAACAGAGGTCAGGTGGTGACAATACCCGCGGCGGACCGCAGAAAAtcacggc                                                                                                                                                                                                      >  2:358979/1‑233 (MQ=255)
                                 cGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTGGGGACAATACCGGCGGCGATCAGCAGGAAATTCAGCg                                                                                                                                                                                                      >  1:826424/1‑239 (MQ=255)
                                 cGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGGACGGTAGTGACAATACCGGCGGCGAGCAGCAGTAAATTCAGCg                                                                                                                                                                                                      >  1:720793/1‑239 (MQ=255)
                                 cGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACAAAAAAACGGGTCGGTAGGGCACATACCGGGGGCGGTCCGCCCGAAAATCAGca                                                                                                                                                                                                      >  2:155679/1‑238 (MQ=255)
                                 cGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAAACCAACAGCGCAACGAGAGAGCCGATACCCGCAGACATCAGAAGAAAAATCAGcc                                                                                                                                                                                                      >  2:302874/1‑238 (MQ=255)
                                 cGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAAAACAACAGCGGGACGGTAGTGACAATACCGGGGGCGATCAGCAGGAAATTCAGCg                                                                                                                                                                                                      >  1:358506/1‑239 (MQ=255)
                                  gTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTGAAAAACAACAGCGGGACGGTAGTGACAATAACCGCGGGGGTCAGGAGGaaaattcacg                                                                                                                                                                                                      >  1:541337/1‑231 (MQ=255)
                                  gTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAAAAACGGGACGGGAGGGACAAAACCGGCGGCGAACAGCAGGAAAATCAGcc                                                                                                                                                                                                      >  1:230837/1‑237 (MQ=255)
                                           gCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCCAGCGCGTGGCAGCGGCGGTAAAAAACCAAAGCGGGACGGTAGAGGCCATACCCGCGGCGGTGAGGCGGAAAATCAGCGAcctaagggt                                                                                                                                                                                            >  2:164876/1‑231 (MQ=255)
                                           gCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGtt                                                                                                                                                                                            <  2:512100/239‑1 (MQ=255)
                                           gCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGACCGGGAGTGACAAAACCGGCGGGGGCCAGCCGAAAATTCAGGCACACCggggt                                                                                                                                                                                            >  1:450422/1‑237 (MQ=255)
                                              tccatcACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCGCCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTg                                                                                                                                                                                         <  1:1053157/239‑1 (MQ=255)
                                                                   gcCACCCAAATAAAGGCGAAAGTCACCAGCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAg                                                                                                                                                                      <  1:718824/237‑1 (MQ=255)
                                                                           cacaaaGGCGAAACGCACCATCTGAGCGGCACCCGGTTTTTCACCATAAAACGTAACAGCCAGCAGGAACATCAGGGTCGGGCCAATGTACTGGAAAAAGCCTAACGTGGAGAGACGCAAGCGCGTGGCCGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACAGCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCaa                                                                                                                                                            <  2:620781/236‑1 (MQ=255)
                                                                                 ggCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAAAAGCGGTACGGGAGTGACCCTTACGGCGGGGGACAACAGAAAATTCAGCGACATCGGGAGTTGCCCCACTAGGGTGGTCGAGATGTCGGCaacaaaca                                                                                                                                                      >  2:106475/1‑233 (MQ=255)
                                                                                     aaaGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGGAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCTGCCATaacaaac‑ac                                                                                                                                                 >  2:948045/1‑238 (MQ=255)
                                                                                            ccATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACCAAACAATGGTACGGGCAGGACAAAAACCGCGGGGGTCAGTAACAAATTCAACGCAACGGGGGTTTGCCCCCACTTGCTGGGTGAGCTGTCGGGAATCACTAAC‑Agcgaacg                                                                                                                                           >  2:240202/1‑232 (MQ=255)
                                                                                                   aTCGACGCCCGGTTTTTCACCATAAAGCGGCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCGGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCa                                                                                                                                   <  2:797611/239‑1 (MQ=255)
                                                                                                     cGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAAAATCGGGAAGGGAGTGACAATACCCTCGGCGATCAGCCGGAAATTCCGCTGCAACGGGGTTTGCCCCATATTGGTGGTTCAGCTGTTGGCAATCACAAAA‑ATGGAAATTGCAGCCACg                                                                                                                                 >  2:134455/1‑239 (MQ=255)
                                                                                                          cccGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGGAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAACAAAA‑AAGGAAATTGCCGCCACGGgca                                                                                                                             >  1:976304/1‑238 (MQ=255)
                                                                                                           ccGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGGAAATTCAGAGACATCCGGTTTTGCCCCATATGGCTGGTTGAACTGTCGGCAATAGCAAAC‑AGGTAAAATGCCGCCACGGgcagc                                                                                                                           >  1:1098458/1‑239 (MQ=255)
                                                                                                                  ttCACCAAAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTCACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCc‑aca                                                                                                                   <  1:554688/239‑1 (MQ=255)
                                                                                                                         tAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGCAGGGACAATATCGGCGGCGGTTAGCAGTAAATTCAGCGACATCGGGCTTTGCCCCaca                                                                                                                                                                                  >  1:875733/1‑169 (MQ=255)
                                                                                                                                  aaagcgagaAGGAACAACAGCGTGGGGCAAAGGAACTGGAAAAAAACTAACGTGGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAATACAACAGCGGTACGGGAGTGACAATTGCGTGGTGGGTCAGCAGTAAATTCAGCGACATCGGGTTTTGCCGCATAAGGCTGGTTGAGCTGTCGGGAATAGCAAAC‑AGTTAAATTTTCGGCACGGGCAGAAGCC‑ACATGGTTTCGATTAAc                                                                                                     <  2:338581/235‑1 (MQ=255)
                                                                                                                                      ccagcAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGAAAAACCAAAAAACCGTACGGGAGAGGCAAATACGGGGGCGGTCAGCAGCAAATTCCGGGACATCGGGGTTTTCCCCATATTGGTGGTTGGGCTGTCCGCAATAGCAAAC‑AGGGAAATTGCCGCCACCGGCAGTCGCC‑ACCGTGTTTCTGTTAACAAG‑ccc                                                                                              >  2:1103836/1‑238 (MQ=255)
                                                                                                                                              aaCATCAGCGTCGGGCCAAGGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTg                                                                                         <  1:164876/236‑1 (MQ=255)
                                                                                                                                                     gCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACCACAGCGGGACGGGAGGGACAATAACGGCGGCGAGCAGCAGGAAATTCAGCGGCATCGGGGTTTGGCCCACATGGCCGGGTGAGATGTCCGCCACAGCAAAC‑AGGGAAATTGCCGCCCCGGGCAGCAGGC‑ACATGGGTTCGCATAAAATG‑CCCGGTTGCCCCTCACCg                                                                               >  2:35987/1‑239 (MQ=255)
                                                                                                                                                     gCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGGAAGGTAGTGACAAAACCGGCGGCGGTCAACAGTAAAATCAGCGACATCGGGGTTTGCCCCATATGGCCGGGTGAGATTTCCGCAATAGCAAAC‑AGGTAAATTACCGCCACGGGGGGCAGCCAACAAGGTTT‑GAATAACaaggccgggtgcgcgtccaca                                                                                >  2:1102995/1‑219 (MQ=255)
                                                                                                                                                       ggcGGTCCAAGGGACTGAAAAAAGCCTAACTGTGAGAGACGCAATCGCGTGGCAGAGGCTGTAAAATACAAGAGAGGGACGTTAGTGACAAGACCGGCGGCGAGCAGCAGAAAATGCAGCGACATCTGGTGTTGCCCCATATGGCGGGTTGAGCGGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGTC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACgt                                                                              <  2:760540/236‑2 (MQ=255)
                                                                                                                                                         cGGGCAAATGTACTGGAAAAAGCCTAACGTTGGGAGGCGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCa                                                                            <  2:976304/238‑1 (MQ=255)
                                                                                                                                                              cAATGTACTGGAAAAAGCCTAACGTAGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAAtcttc                                                                      <  1:948045/239‑1 (MQ=255)
                                                                                                                                                                aTGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATcttctt                                                                    <  1:280155/239‑1 (MQ=255)
                                                                                                                                                                      tGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGcggtacggt                                                                                                                                                                                                                                            >  1:146748/1‑68 (MQ=255)
                                                                                                                                                                      tGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGcggtacggt                                                                                                                                                                                                                                            <  2:146748/68‑1 (MQ=255)
                                                                                                                                                                      tGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAAAACAACAGCGGGACGGGAGAGACAATACCGGCGGCGGGAAGCAGTAAATTCCGCGGCGGCGGGGTTTGCCCCATATGGCGGGTTGGGCGTTCCGCAATAGCAAAA‑CAGTAAATTGCCTCCACGGGGAGACGCA‑CCATGGTTTCGGATAACAAC‑ACGGTTTGGGGGTGAAAGGCAATCTTCTTTCGtc                                                               >  2:1174769/1‑237 (MQ=255)
                                                                                                                                                                      tGGAAAAAGACTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAc                                                              <  2:1085995/239‑1 (MQ=255)
                                                                                                                                                                          aaaaGCCTAACGGTGGGAGGCGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATCCCGGCGGGGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCGTATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGa                                                         <  2:22920/239‑1 (MQ=255)
                                                                                                                                                                             aGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGGACGGGAGTGACAATACCGGCGGCGGTCAGCAGTAAATTCAGCGACATCGGGGTTTGCCCCATATGGCTGGGTGAGCTTTCGGCAATAGCAAAA‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGGTTCGCTTAACAAG‑CCGGTTTGGGCCTCAACGGCCATCTTCTTGGCTAC‑CCGACCg                                                      >  2:977540/1‑239 (MQ=255)
                                                                                                                                                                             aGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAAAAGCGGTACGGTAGTGACAATACCGGCGGCGGTCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCCCGGGGAGCAGCC‑CAATGGTTTCGATTAACATGCCCGGTTTG‑GCCTCAACGGCAATCCTCTTTTGTAC‑Caaaacg                                                      >  2:1008510/1‑234 (MQ=255)
                                                                                                                                                                                       tGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAATCACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCaa                                            <  1:209657/239‑1 (MQ=255)
                                                                                                                                                                                       tGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCa                                             <  1:106475/238‑1 (MQ=255)
                                                                                                                                                                                         agagaCGCAAGCGCGTGGCAGCGGCGGTAAAACACCACAGCGGGACGGTAGTGACAATACCGGCGGCGATCAGCCGTAAATTCAGCGACATCGGGTTTTTCCCCATATGGCTGGGTGAACTGTCGGCAATAACAAACAAGGGAAATT‑TCGCCCCCGGCCGCAGCC‑ACATGGGTTCGCTTAAACTGCCCGGTTTGGGC‑TCAACGGGAAACTTCTTTTGTACACAAACCGAAGAGGGCaaa                                           >  2:656526/1‑239 (MQ=255)
                                                                                                                                                                                         agagaCGCAAGCGCGTGGCAGCGGCGGTAAAAACCAACAGCGGTAAGGGAGTGAAAATAACGGCGGCGATCAGAAATAAAATTAGGGACATAGGGTTTTTCCCCATATGGCTGGTTGTGATGTCGGAGAAAGCAAAA‑AGGTAAATTTTCCCCACGGGAAGAATCC‑CCAACGTTTCCAATTACAAG‑CCGGGTTTCCCGTTAAAGAGAACTTTCTTTCGTAC‑AAAGCCGTAGAAAGCaaa                                           >  2:1073912/1‑238 (MQ=255)
                                                                                                                                                                                          gagaCGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCaaaa                                          <  1:155679/238‑1 (MQ=255)
                                                                                                                                                                                              cGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTaa                                      <  1:722052/238‑1 (MQ=255)
                                                                                                                                                                                               gCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAAt                                    <  1:161928/239‑1 (MQ=255)
                                                                                                                                                                                                 aaGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACGTCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGc                                  <  1:1174769/239‑1 (MQ=255)
                                                                                                                                                                                                           cAGCGGCGGTAAAATACAACAGCGGTACGGGAGTGACAATACCGGCGGCGATCAGCAGTAAATTAAGCGACATCGGGTTTTGCCCCAGATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAg                         <  2:541337/238‑1 (MQ=255)
                                                                                                                                                                                                           cAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAg                         <  1:902853/238‑1 (MQ=255)
                                                                                                                                                                                                           cAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAg                         <  1:358979/238‑1 (MQ=255)
                                                                                                                                                                                                             gcggggGTAAAACACAACGGCGGTACGGTAGTGACAATACCGACGGGGAGCAGCAGGAAATTCAGCGACATCGGGTTTTGCCGCATATGGATGGTTGAGCTGGCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTGCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAgcg                       <  2:12275/238‑1 (MQ=255)
                                                                                                                                                                                                               ggcggTAAAACACAACAGCGGTACGGTAGGGACAATCCCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAgcg                       <  2:937122/236‑1 (MQ=255)
                                                                                                                                                                                                               ggcggTAAAACACAACAGCGGTACGGTAGGGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAgcg                       <  2:230837/236‑1 (MQ=255)
                                                                                                                                                                                                                              caGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGTAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGAt                                                                                                         <  2:204395/141‑1 (MQ=255)
                                                                                                                                                                                                                              caGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGAt                                                                                                         >  1:204395/1‑141 (MQ=255)
                                                                                                                                                                                                                              caGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGaa      >  2:1070579/1‑238 (MQ=255)
                                                                                                                                                                                                                              caGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATA‑CCGGTTTGCGCCTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGCAGGCAAAACTACATGCCAGTCCCAGCGCAATAATAGGTAGCGaa      >  1:559719/1‑238 (MQ=255)
                                                                                                                                                                                                                                 cggtacggtaGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCaa  >  2:290563/1‑239 (MQ=255)
                                                                                                                                                                                                                                 cggtacggtaGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCaa  >  1:1173179/1‑239 (MQ=255)
                                                                                                                                                                                                                                 cggtacggtaGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGCTTTGCCCCATATGGCTGGTTGAACTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAAGTAGCGAACCaa  >  1:674691/1‑239 (MQ=255)
                                                                                                                                                                                                                                  gttactgtaGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCaa  <  1:1020762/236‑1 (MQ=255)
                                                                                                                                                                                                                                  gggacggtcGTGCCAGTCCCGGTGGCGATCAGCAGCAAACTCAGCGACATCGGGTTTTGCACATTATGGCTGGGTGAGCTGTCGGCAATAGCAAAC‑AGGTAAACTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGAGCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTCGCGAAc     <  2:450422/232‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
TATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCC‑ACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAA  >  NC_000913/4003275‑4003738

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