Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,131,769 T→C 7.2% H101R (CAT→CGT)  wcaC ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,131,7690TC 7.2% 211.8 / 2.7 65H101R (CAT→CGT) wcaCputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (34/26);  new base C (2/3);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

ACAGCCCAGCGCCAGCATCTCACGGAATAACTGGCGTTTGCCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACC  >  NC_000913/2131495‑2132022
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                             
aCAGCCCAGCGCCAGCATCTCACGGAATAACTGGCGTTTGCCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACAt                                                                                                                                                                                                                                                               >  2:377469/1‑275 (MQ=255)
    cccAGCGCCAGCATCTCACGGAATAACTGGCGTTTGCCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGaaa                                                                                                                                                                                                                                                           >  2:524802/1‑275 (MQ=255)
         cgcCAGCATCTCACGGAATAACTGGCGTTTGCCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGGCGGTACAGGCGCAGCGTCCGTTACCACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGATGTCCGGTTTAGGGTTTTTCACTTTTTCGCAAAATCGCACCACTCTCTTAAGATTTCGCCAATAGCTTTCCAGCACATGAAAATGCAgaa                                                                                                                                                                                                                                                   >  1:21195/1‑276 (MQ=255)
                            accTGGCGTTTGCCCGCCATCAGTTGGGGTGCGCGATCAATCTTTACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCtgt                                                                                                                                                                                                                                   <  4:377468/273‑1 (MQ=255)
                                         cgcctcccaGTTGGGGTGCGCGATCAACTTGACCCGGCGGATATTTTTTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTg                                                                                                                                                                                                                      <  4:524801/269‑1 (MQ=255)
                                                         tgCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACACGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATCTACCTCATTGAAAtt                                                                                                                                                                                                  >  1:217582/1‑279 (MQ=255)
                                                           cgcgATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGATGTTCAGGTAATGGGCCGACAAAACCCATTGaaa                                                                                                                                                                                                    >  4:165732/1‑275 (MQ=255)
                                                                     cTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAAc                                                                                                                                                                                           <  2:165733/274‑1 (MQ=255)
                                                                        cACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACCCATTGAACTTGCCAAACAGa                                                                                                                                                                                        >  2:394049/1‑274 (MQ=255)
                                                                               cggATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTaaa                                                                                                                                                                             >  1:53234/1‑279 (MQ=255)
                                                                               cggATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCAAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTGTTCGCACAAGCGCACCACGCCCTTAAGATTTACCCATTACCTGTGCAGCTCATGAAAATGCACGACCACCGGACCAGCGTTGCGAGTACTGATGCGATAACCGTAATTAAAATTGCAAAACacaccccgtc                                                                                                                                                                                 >  4:4328/1‑266 (MQ=255)
                                                                                        aTTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAAt                                                                                                                                                                                                                          >  2:87150/1‑225 (MQ=255)
                                                                                        aTTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACCCATTGAAATTGCCAAACAGATCGCGATTaaacaaac                                                                                                                                                                        >  4:314913/1‑275 (MQ=255)
                                                                                           tggggTGGGGAATTTTTGGCCGCGTGTTTTCNANCCGGGGCAACCGGCGGGAAAGGCTCAGGCACCGGTAACACTCCAGNGGTCGTGCGGCGACCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCANGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACa                                                                                                                                                                  <  3:53234/275‑1 (MQ=255)
                                                                                            aaGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTg                                                                                                                                                                                                                      <  4:87150/225‑1 (MQ=255)
                                                                                                   ggCATTTCTGGCAGCCTGTTTTCCAGCCCTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGCCGAAAAAGAGCa                                                                                                                                                             >  2:85282/1‑275 (MQ=255)
                                                                                                                        tCCAGCCTTCGCAACCGTCGGTAAAGGCGCGGCTTCCGGGAACACTCCCGGGGTCGTGCAGCGTCCAGACCAGTGTGACGTCCGGTTTATGGTTTTTGACTGTTTCGCAAAAGCGCACCACCCTAGTAAGATTTAGCCAGTAGCTGTGCAGCACACGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATc                                                                                                                                        <  4:19108/275‑1 (MQ=255)
                                                                                                                           aGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCa                                                                                                                                                             <  1:97027/251‑1 (MQ=255)
                                                                                                                                 cGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCa                                                                                                                                          >  4:347474/1‑264 (MQ=255)
                                                                                                                                  gCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGa                                                                                                                                                                                                                               >  2:94979/1‑178 (MQ=255)
                                                                                                                                     aCCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCcg                                                                                                                                      <  2:347475/264‑1 (MQ=255)
                                                                                                                                      ccGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTaa                                                                                                                                                                                                                           <  4:94979/178‑1 (MQ=255)
                                                                                                                                                aGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATCTTCGCCATCCCGTTCATCCGCGGCGTGTGTTTTTTGACCTGc                                                                                                                >  2:343443/1‑275 (MQ=255)
                                                                                                                                                      aGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGa                                                                                                                                                                                                                                           >  2:466438/1‑146 (MQ=255)
                                                                                                                                                          tCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCg                                                                                                                                                                                                                                       <  4:466437/146‑1 (MQ=255)
                                                                                                                                                                       cAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGAAAGTTCTGATGGCGGACGcc                                                                                         >  4:60585/1‑274 (MQ=255)
                                                                                                                                                                         gTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGAcc                                                                                                                                                                                                                                         >  4:190133/1‑129 (MQ=255)
                                                                                                                                                                             tCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCt                                                                                                                                                                                                                                     <  2:190134/129‑1 (MQ=255)
                                                                                                                                                                               gTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAAt                                                                                                                                                                                                                          >  4:39058/1‑138 (MQ=255)
                                                                                                                                                                               gTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACg                                                                                                                                                                       <  1:168162/189‑1 (MQ=255)
                                                                                                                                                                               gTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCGCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACACATTGAAATTGCCAAACAGATCGCGATTAAACAGACg                                                                                                                                                                       >  3:168162/1‑189 (MQ=255)
                                                                                                                                                                                   ggctgccagaccaggGTGCCGGCCTGTTTTTGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATCTAGCCAGTAGCTGAGCAGGACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTgccgcc                                                                             <  4:404623/270‑1 (MQ=255)
                                                                                                                                                                                   aGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTg                                                                                                                                                                                                                      <  2:39058/138‑1 (MQ=255)
                                                                                                                                                                                       tCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGAAGAAACAAAGCAAACTTCCCCTTCGCGGTACTCCGCTGCTCATGTTTTCGTACCGTCGTGTAGTTCTGATGGCTGACCCCCCCTTTGCCCCCgttc                                                                         >  4:589491/1‑274 (MQ=255)
                                                                                                                                                                                       tCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTCTGCAGCACACGAAAATGCAGGACCACCGGa                                                                                                                                                                                                                                           >  4:213214/1‑113 (MQ=255)
                                                                                                                                                                                          agaccagaGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCg                                                                      >  2:260980/1‑275 (MQ=255)
                                                                                                                                                                                           gaccagaTTGATGTACGGTTTGTGGTTTTTCACTTTTGCGCAAAAGCGGACCACGCGCTTGAGGTTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCTGACACGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCg                                                                      <  4:394048/274‑1 (MQ=255)
                                                                                                                                                                                           gaccagaGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTCTGCAGCACACGAAAATGCAGGACCACCGGACCCg                                                                                                                                                                                                                                       <  2:213215/113‑1 (MQ=255)
                                                                                                                                                                                                 agTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCAGACGCCCCTTTTGCAGCCTTTGCCGTAACCGt                                                               >  2:662685/1‑275 (MQ=255)
                                                                                                                                                                                                 agTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGAGTAAACAGACGAAACAGAGCAATATTCGCCAACGCGGTCACCCGCGGAGTTTGTTTGGTGCCCTGGGGGACGTTTTAATGGCTTCAGCCCTCTCTCCCCCATTTGCCGCAAcgct                                                               >  4:273709/1‑272 (MQ=255)
                                                                                                                                                                                                 agTGAAGGCCGGGTTATGGTTTTTAACTTTTGCGCATAAGAGAACCAGGCGCTTAAGATTTAGCCAGTAGCTGTGCAGCATATGAAAAGGCAGTACCACCGGACTATGTGTGGGAGGAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCg                                                                <  4:343442/274‑1 (MQ=255)
                                                                                                                                                                                                      cGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACTACCGGACCCGCCTTGCGAGTAATGGTTCGATATAACTCATTGACATTGCCAAACAGATCGCGAATAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCCCGGTGTATGTTTGATCACCTGCGGCAAGTTCTGATGGCTGACTATCCCTTTGCCCCCTTTTCCGTACCCGtcaac                                                           >  4:557185/1‑270 (MQ=255)
                                                                                                                                                                                                      cGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCCCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGACAGTTCTGATGGCTGGCGCTCTCTTTGCCGCCTTTGCCGGAACCGTAAAc                                                           >  4:189221/1‑274 (MQ=255)
                                                                                                                                                                                                        tCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACaa                                                         >  2:352675/1‑274 (MQ=255)
                                                                                                                                                                                                           ggTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAAtg                                                     <  2:189222/275‑1 (MQ=255)
                                                                                                                                                                                                           ggTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGAATGTTTGATGACCTGCGGCTAGTTCTGATGGATGACGCTATCTTTGCCGCCTTTGCCGTAACCGTAGAAAAAAtt                                                     >  2:148123/1‑274 (MQ=255)
                                                                                                                                                                                                            gTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAAtgt                                                    <  4:352674/275‑1 (MQ=255)
                                                                                                                                                                                                                 tGGTTTTTCACTTTTTCGCCACAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACACGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCGAACAGATCGCGATTAAACAGACGAATCAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACg                                                <  2:273710/274‑1 (MQ=255)
                                                                                                                                                                                                                  ggTTTTTCACTGTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGcc                                              <  2:626174/275‑1 (MQ=255)
                                                                                                                                                                                                                       ttCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCg                                                                                                                                                                                                                                       >  2:546487/1‑85 (MQ=255)
                                                                                                                                                                                                                       ttCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTaaa                                                                                                                                                                             >  4:161816/1‑143 (MQ=255)
                                                                                                                                                                                                                       ttCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCa                                                                                                                                          >  2:710253/1‑178 (MQ=255)
                                                                                                                                                                                                                           cTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCtgt                                                                                                                                                                                                                                   <  4:546486/85‑1 (MQ=255)
                                                                                                                                                                                                                           cTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGa                                                                                                                                                                         <  2:161817/143‑1 (MQ=255)
                                                                                                                                                                                                                           cTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCcg                                                                                                                                      <  4:710252/178‑1 (MQ=255)
                                                                                                                                                                                                                              tgtaccAAAAGCGAACCAGCGTCTAAAGGTTTGNNNAGTAGCTGTGCAGAACACGAAAATGAAGGACACCCGGACCCGCGGGGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTACGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAgcg                              <  3:217582/274‑1 (MQ=255)
                                                                                                                                                                                                                                            accacGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCtggtg                    <  2:60585/275‑1 (MQ=255)
                                                                                                                                                                                                                                                   tctTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAACGCCCCCTGGTGGAAATc              >  4:137337/1‑274 (MQ=255)
                                                                                                                                                                                                                                                           tttAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGGACACCGGAGCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAGCCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCcac     <  2:4328/275‑1 (MQ=255)
                                                                                                                                                                                                                                                               gCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTtgatg                                                                                                                      >  4:456765/1‑158 (MQ=255)
                                                                                                                                                                                                                                                               gCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCCCGCTGGTGGAGATCTAACGCCACAcc  >  2:297515/1‑274 (MQ=255)
                                                                                                                                                                                                                                                               gCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCCTTGCCGTAACCGTCGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACAcc  >  2:378745/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                   gTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCt                                                                                                                  <  2:456766/158‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gtgCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGc                                               >  2:206498/1‑219 (MQ=255)
                                                                                                                                                                                                                                                                            cagcaCATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACg                          >  1:198905/1‑237 (MQ=255)
                                                                                                                                                                                                                                                                             agcaCATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGc                                           <  4:206497/219‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   tGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGc                         <  2:132927/231‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                             
ACAGCCCAGCGCCAGCATCTCACGGAATAACTGGCGTTTGCCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGAGTAATGGTGCGATATAACTCATTGAAATTGCCAAACAGATCGCGATTAAACAGACGAAACAGAGCAATATTCGCCATCGCGGTCATCCGCGGCGTATGTTTGATGACCTGCGGATAGTTCTGATGGCTGACGCTCTCTTTGCCGCCTTTGCCGTAACCGTAGACAAAATGTGACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACC  >  NC_000913/2131495‑2132022

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