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

TATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAA  >  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)
                                                                                     aaaGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGGAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCTGCCATaacaaacac                                                                                                                                                >  2:948045/1‑238 (MQ=255)
                                                                                            ccATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACCAAACAATGGTACGGGCAGGACAAAAACCGCGGGGGTCAGTAACAAATTCAACGCAACGGGGGTTTGCCCCCACTTGCTGGGTGAGCTGTCGGGAATCACTAACAgcgaacg                                                                                                                                          >  2:240202/1‑232 (MQ=255)
                                                                                                   aTCGACGCCCGGTTTTTCACCATAAAGCGGCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCGGTCGGCAATAGCAAACAGGTAAATTGCCGCCa                                                                                                                                  <  2:797611/239‑1 (MQ=255)
                                                                                                     cGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAAAATCGGGAAGGGAGTGACAATACCCTCGGCGATCAGCCGGAAATTCCGCTGCAACGGGGTTTGCCCCATATTGGTGGTTCAGCTGTTGGCAATCACAAAAATGGAAATTGCAGCCACg                                                                                                                                >  2:134455/1‑239 (MQ=255)
                                                                                                          cccGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGGAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAACAAAAAAGGAAATTGCCGCCACGGgca                                                                                                                            >  1:976304/1‑238 (MQ=255)
                                                                                                           ccGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGGAAATTCAGAGACATCCGGTTTTGCCCCATATGGCTGGTTGAACTGTCGGCAATAGCAAACAGGTAAAATGCCGCCACGGgcagc                                                                                                                          >  1:1098458/1‑239 (MQ=255)
                                                                                                                  ttCACCAAAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTCACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCcaca                                                                                                                   <  1:554688/239‑1 (MQ=255)
                                                                                                                         tAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGCAGGGACAATATCGGCGGCGGTTAGCAGTAAATTCAGCGACATCGGGCTTTGCCCCaca                                                                                                                                                                                 >  1:875733/1‑169 (MQ=255)
                                                                                                                                  aaagcgagaAGGAACAACAGCGTGGGGCAAAGGAACTGGAAAAAAACTAACGTGGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAATACAACAGCGGTACGGGAGTGACAATTGCGTGGTGGGTCAGCAGTAAATTCAGCGACATCGGGTTTTGCCGCATAAGGCTGGTTGAGCTGTCGGGAATAGCAAACAGTTAAATTTTCGGCACGGGCAGAAGCCACATGGTTTCGATTAAc                                                                                                     <  2:338581/235‑1 (MQ=255)
                                                                                                                                      ccagcAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGAAAAACCAAAAAACCGTACGGGAGAGGCAAATACGGGGGCGGTCAGCAGCAAATTCCGGGACATCGGGGTTTTCCCCATATTGGTGGTTGGGCTGTCCGCAATAGCAAACAGGGAAATTGCCGCCACCGGCAGTCGCCACCGTGTTTCTGTTAACAAGccc                                                                                              >  2:1103836/1‑238 (MQ=255)
                                                                                                                                              aaCATCAGCGTCGGGCCAAGGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTg                                                                                         <  1:164876/236‑1 (MQ=255)
                                                                                                                                                     gCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACCACAGCGGGACGGGAGGGACAATAACGGCGGCGAGCAGCAGGAAATTCAGCGGCATCGGGGTTTGGCCCACATGGCCGGGTGAGATGTCCGCCACAGCAAACAGGGAAATTGCCGCCCCGGGCAGCAGGCACATGGGTTCGCATAAAATGCCCGGTTGCCCCTCACCg                                                                               >  2:35987/1‑239 (MQ=255)
                                                                                                                                                       ggcGGTCCAAGGGACTGAAAAAAGCCTAACTGTGAGAGACGCAATCGCGTGGCAGAGGCTGTAAAATACAAGAGAGGGACGTTAGTGACAAGACCGGCGGCGAGCAGCAGAAAATGCAGCGACATCTGGTGTTGCCCCATATGGCGGGTTGAGCGGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGTCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACgt                                                                              <  2:760540/236‑2 (MQ=255)
                                                                                                                                                         cGGGCAAATGTACTGGAAAAAGCCTAACGTTGGGAGGCGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCa                                                                            <  2:976304/238‑1 (MQ=255)
                                                                                                                                                              cAATGTACTGGAAAAAGCCTAACGTAGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAAtcttc                                                                      <  1:948045/239‑1 (MQ=255)
                                                                                                                                                                aTGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATcttctt                                                                    <  1:280155/239‑1 (MQ=255)
                                                                                                                                                                      tGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGcggtacggt                                                                                                                                                                                                                                           >  1:146748/1‑68 (MQ=255)
                                                                                                                                                                      tGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGcggtacggt                                                                                                                                                                                                                                           <  2:146748/68‑1 (MQ=255)
                                                                                                                                                                      tGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAAAACAACAGCGGGACGGGAGAGACAATACCGGCGGCGGGAAGCAGTAAATTCCGCGGCGGCGGGGTTTGCCCCATATGGCGGGTTGGGCGTTCCGCAATAGCAAAACAGTAAATTGCCTCCACGGGGAGACGCACCATGGTTTCGGATAACAACACGGTTTGGGGGTGAAAGGCAATCTTCTTTCGtc                                                               >  2:1174769/1‑237 (MQ=255)
                                                                                                                                                                      tGGAAAAAGACTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTAc                                                              <  2:1085995/239‑1 (MQ=255)
                                                                                                                                                                          aaaaGCCTAACGGTGGGAGGCGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATCCCGGCGGGGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCGTATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGa                                                         <  2:22920/239‑1 (MQ=255)
                                                                                                                                                                             aGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGGACGGGAGTGACAATACCGGCGGCGGTCAGCAGTAAATTCAGCGACATCGGGGTTTGCCCCATATGGCTGGGTGAGCTTTCGGCAATAGCAAAAAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGGTTCGCTTAACAAGCCGGTTTGGGCCTCAACGGCCATCTTCTTGGCTACCCGACCg                                                      >  2:977540/1‑239 (MQ=255)
                                                                                                                                                                             aGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAAAAGCGGTACGGTAGTGACAATACCGGCGGCGGTCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCCCGGGGAGCAGCCCAATGGTTTCGATTAACATGCCCGGTTTGGCCTCAACGGCAATCCTCTTTTGTACCaaaacg                                                      >  2:1008510/1‑234 (MQ=255)
                                                                                                                                                                                       tGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAATCACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCaa                                            <  1:209657/239‑1 (MQ=255)
                                                                                                                                                                                       tGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCa                                             <  1:106475/238‑1 (MQ=255)
                                                                                                                                                                                         agagaCGCAAGCGCGTGGCAGCGGCGGTAAAACACCACAGCGGGACGGTAGTGACAATACCGGCGGCGATCAGCCGTAAATTCAGCGACATCGGGTTTTTCCCCATATGGCTGGGTGAACTGTCGGCAATAACAAACAAGGGAAATTTCGCCCCCGGCCGCAGCCACATGGGTTCGCTTAAACTGCCCGGTTTGGGCTCAACGGGAAACTTCTTTTGTACACAAACCGAAGAGGGCaaa                                           >  2:656526/1‑239 (MQ=255)
                                                                                                                                                                                         agagaCGCAAGCGCGTGGCAGCGGCGGTAAAAACCAACAGCGGTAAGGGAGTGAAAATAACGGCGGCGATCAGAAATAAAATTAGGGACATAGGGTTTTTCCCCATATGGCTGGTTGTGATGTCGGAGAAAGCAAAAAGGTAAATTTTCCCCACGGGAAGAATCCCCAACGTTTCCAATTACAAGCCGGGTTTCCCGTTAAAGAGAACTTTCTTTCGTACAAAGCCGTAGAAAGCaaa                                           >  2:1073912/1‑238 (MQ=255)
                                                                                                                                                                                          gagaCGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCaaaa                                          <  1:155679/238‑1 (MQ=255)
                                                                                                                                                                                              cGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTaa                                      <  1:722052/238‑1 (MQ=255)
                                                                                                                                                                                               gCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAAt                                    <  1:161928/239‑1 (MQ=255)
                                                                                                                                                                                                 aaGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACGTCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGc                                  <  1:1174769/239‑1 (MQ=255)
                                                                                                                                                                                                           cAGCGGCGGTAAAATACAACAGCGGTACGGGAGTGACAATACCGGCGGCGATCAGCAGTAAATTAAGCGACATCGGGTTTTGCCCCAGATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAg                         <  2:541337/238‑1 (MQ=255)
                                                                                                                                                                                                           cAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAg                         <  1:902853/238‑1 (MQ=255)
                                                                                                                                                                                                           cAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAg                         <  1:358979/238‑1 (MQ=255)
                                                                                                                                                                                                             gcggggGTAAAACACAACGGCGGTACGGTAGTGACAATACCGACGGGGAGCAGCAGGAAATTCAGCGACATCGGGTTTTGCCGCATATGGATGGTTGAGCTGGCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTGCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAgcg                       <  2:12275/238‑1 (MQ=255)
                                                                                                                                                                                                               ggcggTAAAACACAACAGCGGTACGGTAGGGACAATCCCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAgcg                       <  2:937122/236‑1 (MQ=255)
                                                                                                                                                                                                               ggcggTAAAACACAACAGCGGTACGGTAGGGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAgcg                       <  2:230837/236‑1 (MQ=255)
                                                                                                                                                                                                                              caGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGTAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGAt                                                                                                         <  2:204395/141‑1 (MQ=255)
                                                                                                                                                                                                                              caGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGAt                                                                                                         >  1:204395/1‑141 (MQ=255)
                                                                                                                                                                                                                              caGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGaa      >  2:1070579/1‑238 (MQ=255)
                                                                                                                                                                                                                              caGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATACCGGTTTGCGCCTCAACGGCAATCTTCTTGCGTACCAGACCGTAGCAGGCAAAACTACATGCCAGTCCCAGCGCAATAATAGGTAGCGaa      >  1:559719/1‑238 (MQ=255)
                                                                                                                                                                                                                                 cggtacggtaGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCaa  >  2:290563/1‑239 (MQ=255)
                                                                                                                                                                                                                                 cggtacggtaGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCaa  >  1:1173179/1‑239 (MQ=255)
                                                                                                                                                                                                                                 cggtacggtaGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGCTTTGCCCCATATGGCTGGTTGAACTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAAGTAGCGAACCaa  >  1:674691/1‑239 (MQ=255)
                                                                                                                                                                                                                                  gttactgtaGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCaa  <  1:1020762/236‑1 (MQ=255)
                                                                                                                                                                                                                                  gggacggtcGTGCCAGTCCCGGTGGCGATCAGCAGCAAACTCAGCGACATCGGGTTTTGCACATTATGGCTGGGTGAGCTGTCGGCAATAGCAAACAGGTAAACTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGAGCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTCGCGAAc     <  2:450422/232‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         
TATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATCGGGTTTTGCCCCATATGGCTGGTTGAGCTGTCGGCAATAGCAAACAGGTAAATTGCCGCCACGGGCAGCAGCCACATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAA  >  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: