Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,003,509 T→G 22.3% 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,5090TG22.3% 101.1 / 2.5 59T224P (ACC→CCC) rarDputative chloramphenical resistance permease
Reads supporting (aligned to +/- strand):  ref base T (16/26);  new base G (14/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 9.76e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAA  >  NC_000913/4003276‑4003739
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         
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‑AGGTAAATTGCCGCCACGGGCAGCAGCcaca                                                                                                                   <  1:554688/239‑1 (MQ=255)
                                                                                                                         tAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGCAGGGACAATATCGGCGGCGGTTAGCAGTAAATTCAGCGACATCGGGCTTTGCCCCaca                                                                                                                                                                                 >  1:875733/1‑169 (MQ=255)
                                                                                                                                  aaagcgagaAGGAACAACAGCGTGGGGCAAAGGAACTGGAAAAAAACTAACGTGGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAATACAACAGCGGTACGGGAGTGACAATTGCGTGGTGGGTCAGCAGTAAATTCAGCGACATCGGGTTTTGCCGCATAAGGCTGGTTGAGCTGTCGGGAATAGCAAAC‑AGTTAAATTTTCGGCACGGGCAGAAGCCACATGGTTTCGATTAAc                                                                                                     <  2:338581/235‑1 (MQ=255)
                                                                                                                                      ccagcAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGAAAAACCAAAAAACCGTACGGGAGAGGCAAATACGGGGGCGGTCAGCAGCAAATTCCGGGACATCGGGGTTTTCCCCATATTGGTGGTTGGGCTGTCCGCAATAGCAAAC‑AGGGAAATTGCCGCCACCGGCAGTCGCCACCGTGTTTCTGTTAACAAG‑ccc                                                                                              >  2:1103836/1‑238 (MQ=255)
                                                                                                                                              aaCATCAGCGTCGGGCCAAGGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTg                                                                                         <  1:164876/236‑1 (MQ=255)
                                                                                                                                                     gCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACCACAGCGGGACGGGAGGGACAATAACGGCGGCGAGCAGCAGGAAATTCAGCGGCATCGGGGTTTGGCCCACATGGCCGGGTGAGATGTCCGCCACAGCAAAC‑AGGGAAATTGCCGCCCCGGGCAGCAGGCACATGGGTTCGCATAAAATG‑CCCGGTTGCCCCTCACCg                                                                               >  2:35987/1‑239 (MQ=255)
                                                                                                                                                       ggcGGTCCAAGGGACTGAAAAAAGCCTAACTGTGAGAGACGCAATCGCGTGGCAGAGGCTGTAAAATACAAGAGAGGGACGTTAGTGACAAGACCGGCGGCGAGCAGCAGAAAATGCAGCGACATCTGGTGTTGCCCCATATGGCGGGTTGAGCGGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGTCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACgt                                                                              <  2:760540/236‑2 (MQ=255)
                                                                                                                                                         cGGGCAAATGTACTGGAAAAAGCCTAACGTTGGGAGGCGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCa                                                                            <  2:976304/238‑1 (MQ=255)
                                                                                                                                                              cAATGTACTGGAAAAAGCCTAACGTAGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAAtcttc                                                                      <  1:948045/239‑1 (MQ=255)
                                                                                                                                                                aTGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATcttctt                                                                    <  1:280155/239‑1 (MQ=255)
                                                                                                                                                                      tGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGcggtacggt                                                                                                                                                                                                                                           >  1:146748/1‑68 (MQ=255)
                                                                                                                                                                      tGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGcggtacggt                                                                                                                                                                                                                                           <  2:146748/68‑1 (MQ=255)
                                                                                                                                                                      tGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAAAACAACAGCGGGACGGGAGAGACAATACCGGCGGCGGGAAGCAGTAAATTCCGCGGCGGCGGGGTTTGCCCCATATGGCGGGTTGGGCGTTCCGCAATAGCAAAA‑CAGTAAATTGCCTCCACGGGGAGACGCACCATGGTTTCGGATAACAAC‑ACGGTTTGGGGGTGAAAGGCAATCTTCTTTCGtc                                                               >  2:1174769/1‑237 (MQ=255)
                                                                                                                                                                      tGGAAAAAGACTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAc                                                              <  2:1085995/239‑1 (MQ=255)
                                                                                                                                                                          aaaaGCCTAACGGTGGGAGGCGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATCCCGGCGGGGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCGTATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGa                                                         <  2:22920/239‑1 (MQ=255)
                                                                                                                                                                             aGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGGACGGGAGTGACAATACCGGCGGCGGTCAGCAGTAAATTCAGCGACATCGGGGTTTGCCCCATATGGCTGGGTGAGCTTTCGGCAATAGCAAAA‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGGTTCGCTTAACAAG‑CCGGTTTGGGCCTCAACGGCCATCTTCTTGGCTAC‑CCGACCg                                                      >  2:977540/1‑239 (MQ=255)
                                                                                                                                                                             aGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAAAAGCGGTACGGTAGTGACAATACCGGCGGCGGTCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCCCGGGGAGCAGCCCAATGGTTTCGATTAACATGCCCGGTTTG‑GCCTCAACGGCAATCCTCTTTTGTAC‑Caaaacg                                                      >  2:1008510/1‑234 (MQ=255)
                                                                                                                                                                                       tGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAATCACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCaa                                            <  1:209657/239‑1 (MQ=255)
                                                                                                                                                                                       tGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCa                                             <  1:106475/238‑1 (MQ=255)
                                                                                                                                                                                         agagaCGCAAGCGCGTGGCAGCGGCGGTAAAACACCACAGCGGGACGGTAGTGACAATACCGGCGGCGATCAGCCGTAAATTCAGCGACATCGGGTTTTTCCCCATATGGCTGGGTGAACTGTCGGCAATAACAAACAAGGGAAATT‑TCGCCCCCGGCCGCAGCCACATGGGTTCGCTTAAACTGCCCGGTTTGGGC‑TCAACGGGAAACTTCTTTTGTACACAAACCGAAGAGGGCaaa                                           >  2:656526/1‑239 (MQ=255)
                                                                                                                                                                                         agagaCGCAAGCGCGTGGCAGCGGCGGTAAAAACCAACAGCGGTAAGGGAGTGAAAATAACGGCGGCGATCAGAAATAAAATTAGGGACATAGGGTTTTTCCCCATATGGCTGGTTGTGATGTCGGAGAAAGCAAAA‑AGGTAAATTTTCCCCACGGGAAGAATCCCCAACGTTTCCAATTACAAG‑CCGGGTTTCCCGTTAAAGAGAACTTTCTTTCGTAC‑AAAGCCGTAGAAAGCaaa                                           >  2:1073912/1‑238 (MQ=255)
                                                                                                                                                                                          gagaCGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCaaaa                                          <  1:155679/238‑1 (MQ=255)
                                                                                                                                                                                              cGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTaa                                      <  1:722052/238‑1 (MQ=255)
                                                                                                                                                                                               gCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAAt                                    <  1:161928/239‑1 (MQ=255)
                                                                                                                                                                                                 aaGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACGTCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGc                                  <  1:1174769/239‑1 (MQ=255)
                                                                                                                                                                                                           cAGCGGCGGTAAAATACAACAGCGGTACGGGAGTGACAATACCGGCGGCGATCAGCAGTAAATTAAGCGACATCGGGTTTTGCCCCAGATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAg                         <  2:541337/238‑1 (MQ=255)
                                                                                                                                                                                                           cAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAg                         <  1:902853/238‑1 (MQ=255)
                                                                                                                                                                                                           cAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAg                         <  1:358979/238‑1 (MQ=255)
                                                                                                                                                                                                             gcggggGTAAAACACAACGGCGGTACGGTAGTGACAATACCGACGGGGAGCAGCAGGAAATTCAGCGACATCGGGTTTTGCCGCATATGGATGGTTGAGCTGGCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTGCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAgcg                       <  2:12275/238‑1 (MQ=255)
                                                                                                                                                                                                               ggcggTAAAACACAACAGCGGTACGGTAGGGACAATCCCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAgcg                       <  2:937122/236‑1 (MQ=255)
                                                                                                                                                                                                               ggcggTAAAACACAACAGCGGTACGGTAGGGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAgcg                       <  2:230837/236‑1 (MQ=255)
                                                                                                                                                                                                                              caGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGTAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGAt                                                                                                         <  2:204395/141‑1 (MQ=255)
                                                                                                                                                                                                                              caGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGAt                                                                                                         >  1:204395/1‑141 (MQ=255)
                                                                                                                                                                                                                              caGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGaa      >  2:1070579/1‑238 (MQ=255)
                                                                                                                                                                                                                              caGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATA‑CCGGTTTGCGCCTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGCAGGCAAAACTACATGCCAGTCCCAGCGCAATAATAGGTAGCGaa      >  1:559719/1‑238 (MQ=255)
                                                                                                                                                                                                                                 cggtacggtaGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCaa  >  2:290563/1‑239 (MQ=255)
                                                                                                                                                                                                                                 cggtacggtaGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCaa  >  1:1173179/1‑239 (MQ=255)
                                                                                                                                                                                                                                 cggtacggtaGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGCTTTGCCCCATATGGCTGGTTGAACTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAAGTAGCGAACCaa  >  1:674691/1‑239 (MQ=255)
                                                                                                                                                                                                                                  gttactgtaGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCaa  <  1:1020762/236‑1 (MQ=255)
                                                                                                                                                                                                                                  gggacggtcGTGCCAGTCCCGGTGGCGATCAGCAGCAAACTCAGCGACATCGGGTTTTGCACATTATGGCTGGGTGAGCTGTCGGCAATAGCAAAC‑AGGTAAACTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGAGCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTCGCGAAc     <  2:450422/232‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         
TATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAAC‑AGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATG‑CCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAC‑CAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAA  >  NC_000913/4003276‑4003739

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: ‑