Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,511,862 T→G 18.1% V70G (GTG→GGG)  ydcS → putative ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,511,8620TG18.1% 150.4 / 3.2 80V70G (GTG→GGG) ydcSputative ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (48/14);  new base G (8/9);  total (56/24)
Fisher's exact test for biased strand distribution p-value = 3.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGCACACCTTACCAATGGGGG  >  NC_000913/1511639‑1512088
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      
cAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAAAAATACGACCGGGTAACGCAGTTCGAAAAAGAGAAGGGCTAGCGGGTGAATGGGAAAAccgccaggaa                                                                                                                                                                                                                       >  2:987413/1‑234 (MQ=255)
cAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCCGGATACATCGAACCCGGACAAACAGATAAACAATACGACTGGGTAAACCCGTTCGAAAAAAAGACAGGCAGCGCGAGGAATTGGAAAAcccacccgcg                                                                                                                                                                                                                        >  2:194717/1‑231 (MQ=255)
   gAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTCCCACTTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGCTAAACAATACGAATGGGTAACGCGGTTCTAAAAAAAGACAGGGCGCGCGGGGAGTTGGAAAACCCCcgcgcatc                                                                                                                                                                                                                     >  2:558610/1‑234 (MQ=255)
      tAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGa                                                                                                                                                                                                                 >  1:624230/1‑239 (MQ=255)
          actttAGGAGCAAAACATTGCCCCGCAGCAGCCTGTGTGCGCTCAGCATAACAATAAGGACCGCCACCGCCGCCGACCCGCCTCCCAATTAAGATAAACCGGAAGGGCGACTGGATATTATCGCCGGGCCGGGATACATCGAACGCGGAAAAACTGAAAAACAATCCGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGaa                                                                                                                                                                                                             <  1:878955/236‑1 (MQ=255)
             ttATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAAAAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCAGAGGaaagg                                                                                                                                                                                                          >  2:698974/1‑237 (MQ=255)
              tattagcaacggaTTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGCCCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCAACGGGCAATTATCGCCTGGCCGGGCTACATCGAACACGGCCAAGCTGTGAAACAATACGCCTGGGACACGCAGATGGAAAAGAGACCCGGCCACACGGTGAATTGAAAAACCCCCGCACCTTcccaaaaaagga                                                                                                                                                                                                         >  2:1040107/13‑229 (MQ=255)
              tATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGCAAAACTGATAACCCATACGCCCGGGCAAAGCAGTTCGAAAAAACAACAGGCCCCCCGGGGAAAGGGGAAAccccccccgaacccgaagaaacgg                                                                                                                                                                                                         >  2:959842/1‑220 (MQ=255)
              tATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGCAAAAGAGACAGGCTTCCCGGGGAATGGGAAAAACGCCCCGACCTCCCAAGAAAAgg                                                                                                                                                                                                         >  2:994882/1‑239 (MQ=255)
              tATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAAAACGACTGGGTACCGCAGTTCGAAAAAGAGACCGGCTGCGCGGAGAATGTGAAAACAGACGCGACCTCCAATGAAATgg                                                                                                                                                                                                         >  2:435008/1‑239 (MQ=255)
                    cAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAAAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGCTGAAATGGTCAg                                                                                                                                                                                                     >  1:966116/1‑237 (MQ=255)
                       gACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGCAAACCGCCGCGACTTCCGATGAAATGGTCAGTAtga                                                                                                                                                                                                >  1:518244/1‑239 (MQ=255)
                            ttGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGCCAACTGGATATTATCGCCTGCCCGGGATACATCCAACCCGGCCAAAATGTTAACCAATCGAATGGGGAAAGCCAATTCGAAAAGGAGACAGGCGCCGCGGCTAATTGGAAAACACCCGCGACATCCGAGAAAAGGGtaacccttgtgaca                                                                                                                                                                                           >  2:403155/1‑226 (MQ=255)
                                               tgCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCGGATGACCAAAGGGGGTTAAGATCTgt                                                                                                                                                                        >  1:473333/1‑238 (MQ=255)
                                                       gCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGCATGGGTAACGCAGTTCGAAAAAAAGACAGGCTTCGCGTTGAAGGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGAATGCCAAAGGGGCTAACAATCGGGGTACGGc                                                                                                                                                                 >  2:207327/1‑238 (MQ=255)
                                                       gCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACAGCCGCGACTTCCGATGAAATGGTCAGCCTGATGACCAAAGGGGGTTACGACCTGGTTACGGc                                                                                                                                                                 >  1:791457/1‑238 (MQ=255)
                                                           gACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACCAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCg                                                                                                                                                           >  1:1045534/1‑239 (MQ=255)
                                                              ataaatgcccccgccccccgccgAACCGCCTACAAATTTAAAAAAACGGGAAGGGGCAGTGGATCTTACCGCCGGGCAGGAAAAAGTCGAACGAGGAAAAGCTGAAAACAAATAGGACTGGTAAACGAAGTACGAAAAAGAGACAGGCTGCGGGGGGAATGGAAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCg                                                                                                                                                        <  2:518244/222‑1 (MQ=255)
                                                                taatGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGAt                                                                                                                                                      >  1:175470/1‑239 (MQ=255)
                                                                taatGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCAGCGCGGCGAATGTGAAAACCGCCGCGACTGCCGATGAAATGGTCAGTCTGGTGACCAAAGGGGGTTACGATCTGGTAACGGAATCCgggacg                                                                                                                                                      >  2:740465/1‑235 (MQ=255)
                                                                  atGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGTTg                                                                                                                                                     >  1:409389/1‑238 (MQ=255)
                                                                       cGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGTATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGc                                                                                                                                               >  1:168379/1‑239 (MQ=255)
                                                                         cTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCCCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGTGGGTTACGATCTGGTTACGGCATCCGCCGATGCCTCGCt                                                                                                                                              >  1:1074387/1‑238 (MQ=255)
                                                                          tCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACGCCCGCGACTTCCGATGAAATGGTCATGCTGAGAGCCAAAGGGGGTTACGTTCTGGTTACGGAGACCGGCGagacccccggcc                                                                                                                                            >  2:735742/1‑228 (MQ=255)
                                                                          tCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCAGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGCCGATGCCTCGCAGc                                                                                                                                            >  1:999438/1‑239 (MQ=255)
                                                                             ggcccccGACCCGCCAACGAAGTTAGAAAAGCGGGGAGGGGGAATGGTTATGTTGGGCGGGCCGGGAAGAATGGAAGGCGGAAAAACGGATACCGAAGACGTAGGGAGAAAGGAGTGGAAAAAAGGGACAGGGGGCGGGGTGAAGGAAAAAACCGCCGCGTCTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGtt                                                                                                                                         <  2:791457/238‑1 (MQ=255)
                                                                                 gccgAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGCCTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAAACGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGAACTGGTTACGGCATCCGGCGAGGCCCCGCCGCGGTTGAt                                                                                                                                     >  2:1063683/1‑239 (MQ=255)
                                                                                                    aaaaaaaaCAGGAAGGGCGACGGAATATTACCGCCGGGCGGGAAAACATCGAACGCGGAAAAACTGATAAACAATACGACTGGGAAACGCAGTTCAAAAAGGGGACAGGCTGCGCGGGGAATGTAAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAg                                                                                                                  <  1:1063683/234‑1 (MQ=255)
                                                                                                       aTAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAACCTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATTGGAAAACCGCCGCGACTTCCGAGGAAAGGGTCAGTCTGATGACCAACGGGGGTTCCGATCTGGTTACGGACTCCAGGCGGGCCGCGATGGGGTTGATTATGGGTaaaaaggggaaaccg                                                                                                               >  2:521406/1‑227 (MQ=255)
                                                                                                       aTAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCg                                                                                                               >  1:1049326/1‑239 (MQ=255)
                                                                                                         agcccGGAGGGGAGACGGAATTTTGCCGCCTGGCCGGAAACCAGCAAAGGGGGAAAAACTGAAAAAAAAGCCGACGGGGAACCGCAGTCCAAAAAAGAGACAGGCTGCGCGGTGAATGGGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGAt                                                                                                             <  1:46825/236‑1 (MQ=255)
                                                                                                         aacccGGAAGGGCGATTGTATATTCACGCCGGGCGGGAATACAACGAACGGGGACAAACTGAAAAACAATCCGACGGGGTAACGCGGTTCGAAAAAGGGACGGGGTGCGGGGTGAATGTGAAACCCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGAt                                                                                                             <  1:435008/236‑1 (MQ=255)
                                                                                                         aaaCCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCCAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTCCGTTCTGGTTACGGCATCCGGCGATGCCTCGCTGGGTTTGGTTATGGTGAAACGCGGGccccccct                                                                                                             >  2:1013378/1‑232 (MQ=255)
                                                                                                                aaGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGAACTGGTTCAGGCATCCGGCGATGCCTCCCTTCGTTTGATTATGGGGAAAAGGGGGCAGCCGCTAAAAAcc                                                                                                      >  2:470684/1‑239 (MQ=255)
                                                                                                                       aCTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGACTCCGGCGATGCCTCCCTGCGTTTGATTATGGGTAAACCCGGGCACCCGATAAAAACCCCa                                                                                                   >  2:46825/1‑235 (MQ=255)
                                                                                                                                tatCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGGTTCCCa                                                                                        >  1:429005/1‑236 (MQ=255)
                                                                                                                                tatCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAACCGATTAATACCGCCTTGATTCCCa                                                                                        >  1:1081955/1‑236 (MQ=255)
                                                                                                                                  tCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTgg                                                                                   >  1:803806/1‑239 (MQ=255)
                                                                                                                                  tCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGAGCGGTGCATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGCTCTGGTTACGGCATCCGGCGATGGCCACCTGCGTTTGATTTTGGGTAAACGCGTGCAGCCGATTAATACCCACTTGTTTCCCAAATgg                                                                                   >  2:341701/1‑239 (MQ=255)
                                                                                                                                     ccTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATCACGCATTGTTTCCCACCTGGaaa                                                                                >  1:433174/1‑239 (MQ=255)
                                                                                                                                        ggCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGGGGGAAGCCGCTAAAAACCGCATTGATTCCAAACGGGaaaaa                                                                              >  2:528904/1‑237 (MQ=255)
                                                                                                                                           cGGGATCCATCCAACGGGGAAAAAGTGAAAAACAATGCGACGGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGGGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCt                                                                           <  1:757799/238‑1 (MQ=255)
                                                                                                                                                   acagaacgggaaaaaacgaaaaaacaaacCGACGGGGAAACCAAGTTAAAAAAAGGGACGGGCTGGGGGGGGAATGGAAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCcgcgc                                                                 <  1:994882/218‑1 (MQ=255)
                                                                                                                                                    tggAACGGGGAAAAACTGAAAAAAAATACGACTGGATAACGAGTTTCAAAAAAGAGACAGGCTGCGGGGGGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCcgcgc                                                                 <  1:528904/237‑1 (MQ=255)
                                                                                                                                                    tagacccgggaaaactggaaaaaaaaTCCGATGGGGTAACGCAGTCCAAAAAGGGGACAGTCTGCGGGGGAAAGGGAAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAAGACCGCATTGATTCCCAACTGGAAAACGCTCGATCcgcgc                                                                 <  1:987413/219‑1 (MQ=255)
                                                                                                                                                    tCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAACGGGGGTTACGATCTGGTTACGGCATCCGGCGATTCCACCCTGCGTTTGGTTTTGGGTAAACGCGTGCCGCCGGTAAAAAACCACTTGTTTCCCACCTGGAAAACGCCCGACCcccgc                                                                 >  2:757799/1‑239 (MQ=255)
                                                                                                                                                    acGAACGCGGACAAACTGAAAAAAAAGACGACGGGATAACGCAGGTAGAAAAAGAGACAGGCTGGGGGGTAAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCcgcgc                                                                 <  1:959842/238‑1 (MQ=255)
                                                                                                                                                        aCGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGc                                                                                                                                                    >  1:295735/1‑153 (MQ=255)
                                                                                                                                                                   tgtaaaaaCAACCCGAGGGGGTACCCGAGTCCAAAAAAGGGGCGGGGTGCGGGGTAATGGGAAAAGCCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGAc                                                  <  1:698974/234‑1 (MQ=255)
                                                                                                                                                                            aaTACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAAt                                         >  1:91043/1‑239 (MQ=255)
                                                                                                                                                                              tgcGACTGGGTAACGAATGGCAAAAAAGAGCCAGGCTGCGAGGTGAATGGGAAACCCGGCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGt                                       <  1:194717/237‑1 (MQ=255)
                                                                                                                                                                              tACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGCTCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTTTGGGTAAACGCGGGCAGCCGGTTAATACCGTATTGCTTCCCACCTGGAAAATCCTCGTTCCGCGCGTGGGTAAACACCACTGGTTTAAGGt                                       >  2:250259/1‑239 (MQ=255)
                                                                                                                                                                              tACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGt                                       >  1:298195/1‑239 (MQ=255)
                                                                                                                                                                              tACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGGGGTTAAAGGCGACTGGTTTAATGt                                       >  1:204331/1‑239 (MQ=255)
                                                                                                                                                                              tACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTCCGTTTGATTATGGGTAAACCCGTGCAGCCCATTAATCCCGCATTGATTTCCAACTGAAAAAAACCCGATCCGCGCGAGGTTAAAggggcaggggttaagtt                                       >  2:919040/1‑224 (MQ=255)
                                                                                                                                                                              tACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCAGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGGTGCATCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGTTTCCAAACTGGAAAACGGTCGATCCGCCCAGGATTAAACGCGACCGGTTAAAAGt                                       >  2:168927/1‑239 (MQ=255)
                                                                                                                                                                                cGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGCTGCCTCGCTTCGTTTGATTATGGGTAAACGCGGGCAGCCGCTTAATACCCCATTTATTCCCAACTGGAAAACGCTGGATCCCCACGGGGTTAAAGGGGACCGGTTTAttgttg                                     >  2:694514/1‑239 (MQ=255)
                                                                                                                                                                                cGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAg                                                                                                                  >  2:934259/1‑163 (MQ=255)
                                                                                                                                                                                cGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAg                                                                                                                  <  1:934259/163‑1 (MQ=255)
                                                                                                                                                                                cGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCCCGCGTGGTTAAAGGCGACTGGTTTAATGTTg                                     >  1:991064/1‑239 (MQ=255)
                                                                                                                                                                                      ggTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAACGGGGGGTCCGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAAAGCGTGCCGCCGATTAAAAACGCCTTTATTCCCAACTGGAAAAAGGCGGAACCGCCCGCGGGTAAAGGCAACGTGTTTAATGTTGCCTGCa                               >  2:511932/1‑239 (MQ=255)
                                                                                                                                                                                      ggTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATa                                                                                                        >  1:206569/1‑167 (MQ=255)
                                                                                                                                                                                      ggTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCACCTGGAAAAGGCTCGATCGGCAGGTGGTTAAAGGCGACTGGTTTAAGTTTGGCGGCa                               >  1:755014/1‑239 (MQ=255)
                                                                                                                                                                                      ggTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGGGCAGCCGATTAATACCGCATTTATTCCCAACTGGAAAACCCTCGATCCGCCCGTGGTTAAAGGGGACCGGTTTAATGTGGGGGGCa                               >  2:441542/1‑239 (MQ=255)
                                                                                                                                                                                        taaacgaattccGAAAAACAGAAGGGGTGCGCGTGAATTGGAAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCa                               <  1:694514/227‑1 (MQ=255)
                                                                                                                                                                                                 ttCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAAc                                                                                                                          >  1:153982/1‑139 (MQ=255)
                                                                                                                                                                                                 ttCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGAtt                                                                                             >  1:405381/1‑168 (MQ=255)
                                                                                                                                                                                                 ttCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCCCGTGGTTAAAGGCGACTGGTTTAATTTTTGCGGCAAAGTTTACg                      >  1:830499/1‑238 (MQ=255)
                                                                                                                                                                                                 ttCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAAATCCGCATTGATTCCCCAACTGAAAACCCCCGGCCCGCGACGGGGTAAAGGCGAACGGTTTAAAGATGGAGGGAAATGTTACgg                     >  2:533519/1‑239 (MQ=255)
                                                                                                                                                                                                  tcaaaaaaGTGACGGGGGGCCGGGTGAATTGAAAACCCCGCGCGAGTTACGGATAAAGTGTCAGGCTGATTGCCAAAGGGGGTAACGCTCTGGATACGGCATCCGGCAATGCCTCACTGCGTTTGATTATGGGGAAACCGGTGCAGCCGATTAATACCGCAGTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACgg                     <  2:999438/235‑1 (MQ=255)
                                                                                                                                                                                                    tcaaaaaaggagcacggttcggcgGTAAATGGGAAACCCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACgg                     <  1:376583/218‑1 (MQ=255)
                                                                                                                                                                                                    aaaaaaGAGACGGGTGGGGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGca                   <  1:740465/237‑1 (MQ=255)
                                                                                                                                                                                                     gaaaaaatgcacgggCTGCGCGGTAAATGTGAAACCCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGca                   <  1:470684/226‑1 (MQ=255)
                                                                                                                                                                                                     cgaaaaaggacttgggtGCGGGGGGAATGTGAAGACCGGCGCTGCTTCCGATGAAACGGGCAGTCGGATGACCAAAGGGGGTTAGGATCGGGTTCCGGCAGCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCAGTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGca                   <  2:1045534/223‑1 (MQ=255)
                                                                                                                                                                                                     aaaaaGAGACAGGCTGCGGGGTAAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGcaca                 <  1:735742/239‑1 (MQ=255)
                                                                                                                                                                                                      aaaaaagggaagaggggggggggTGAGTGGAAAACCCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGca                   <  2:803806/218‑1 (MQ=255)
                                                                                                                                                                                                       aaaaaagcgacgtggcgtggggagAATGTGAAACCCGGCTGGTCTTCCGGTTAAAGGGGCAGTCTGATGACCAAAGGGGGTTACGTTCTGGTACCGGAATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTTATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGTGACCGGTTTAATGTTGGCGGCAAAGTTTACGGca                   <  2:1081955/215‑1 (MQ=255)
                                                                                                                                                                                                        aaGGGACAGGGCTGCGGGGTAAATGTAAAACCCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGCACAcc               <  1:974536/239‑1 (MQ=255)
                                                                                                                                                                                                                        cgcgGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGCACACCTTCCCAATggggg  >  1:735850/1‑238 (MQ=255)
                                                                                                                                                                                                                        cgcgGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGCACACCTTACCAATggggg  >  1:125664/1‑238 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      
CAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGTCTGATGACCAAAGGGGGTTACGATCTGGTTACGGCATCCGGCGATGCCTCGCTGCGTTTGATTATGGGTAAACGCGTGCAGCCGATTAATACCGCATTGATTCCCAACTGGAAAACGCTCGATCCGCGCGTGGTTAAAGGCGACTGGTTTAATGTTGGCGGCAAAGTTTACGGCACACCTTACCAATGGGGG  >  NC_000913/1511639‑1512088

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