Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,091,193 T→A 19.3% I392F (ATT→TTT)  frvB ← putative PTS enzyme, IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,091,1930TA19.3% 163.6 / 8.7 76I392F (ATT→TTT) frvBputative PTS enzyme, IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base A (12/3);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 7.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTTTTATTGATGGGGCCGCCCATATCGAACGCCAGCATCGCA  >  NC_000913/4090958‑4091431
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
cggcgAAACCACAATCCTTCCAGGACGATAAGGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGAAAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATCCCCAAACCAGGTGCCAGACGTTTGATCCCGAAGGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAg                                                                                                                                                                                                                                                <  2:603093/238‑1 (MQ=255)
 ggcgTAACCCCAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGCCGTTTGATTCCGAAGGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGc                                                                                                                                                                                                                                               <  1:1034692/238‑1 (MQ=255)
    gTACCCCCAATCCTTCCAGGAGAATAATGATTGAAATGTTCCCCCCGCGGCCAACCGCCAGCCCCAAAAAAAACGAACCCCCCGCGGACACCAGCCCGATCGGCGGTTCAAAAATCCCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGACCACCAGCACGCCGGAAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCcggt                                                                                                                                                                                                                                           <  2:993940/239‑1 (MQ=255)
     tAACCACAATCCTTTCAGGACGAAAAGGAAGGAAATGTTAAGCGCGAGGCAAACCGCCAGCACCAGATAAAACGACCCCACCGGCGACATCAGCCCGATCAGCGGGTCAAAGATCCCCAAACCAGGGGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCcggt                                                                                                                                                                                                                                           <  2:360662/238‑1 (MQ=255)
            aaTCCTTCCAGGACGTAAATGAAGAAAATGTTCAGCGGGGGGCCAATCCCCAGCACCAAATAAAACGAACCCACCGGCGAAAACCGGCCGATCCGCGGGTAAAAATCACCCAAACCAGGTCCCGGCCGTTGGATTCCGAAAGCGATCCCCGGCACGCCGTTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGc                                                                                                                                                                                                                                   <  2:1127668/239‑1 (MQ=255)
             aTCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCCACGGAGCTGCCCGGGGGGAAGGAAAAACcgggtcgggttgg                                                                                                                                                                                                                                  >  2:86710/1‑234 (MQ=255)
                 tttCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGCAATCATCGGCAACAGGACTGCCAGCGCGTAAGGAAAAGCCGATTCTGTTGCGGCtg                                                                                                                                                                                                                              >  2:761677/1‑239 (MQ=255)
                   tCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGCGCCGCCCGAGCGGAAGGAATAAACAGTTCGGGTgaggccgtg                                                                                                                                                                                                                            >  2:507642/1‑231 (MQ=255)
                         cGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCAcc                                                                                                                                                                                                                       <  2:130620/238‑1 (MQ=255)
                          gATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAAACGAGCTGCCAGGGGGCAAGGAAGAGCCGGGTCGGCTGCGAGTGGGGCACCg                                                                                                                                                                                                                      >  1:1043625/1‑238 (MQ=255)
                          gATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATACTCGGCAACGGAGGTGCCCGAGCGGAAGGAAAAGCCGGTTCGGTTGCGGCCTTGGCCCCg                                                                                                                                                                                                                      >  1:1062046/1‑238 (MQ=255)
                            tAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGGAAGGAATAGc                                                                                                                                                                                                                                               >  1:153950/1‑211 (MQ=255)
                            tAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGGCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATAAACCGCAACGGAGACGCCGGAGCGTAAGGAATAGCCGGTTCGGGTGCGGCAGTgggcccagac                                                                                                                                                                                                                   >  2:406278/1‑231 (MQ=255)
                              atgaatgaAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACcac                                                                                                                                                                                                                 <  1:1058207/239‑1 (MQ=255)
                                  atgaaATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATg                                                                                                                                                                                                             <  1:897649/239‑1 (MQ=255)
                                     aaaTGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTAAAAGATCCCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATgg                                                                                                                                                                                                            <  1:1155817/237‑1 (MQ=255)
                                      aaTGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAAAGGAGCTGCCCGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCa                                                                                                                                                                                                          >  1:431799/1‑238 (MQ=255)
                                      aaTGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGAAATCATCGGCAACGGAGCTGCCCGAGCGGAAGGAAAAGCCGGGTTGGGTGGGACTGTGGGGCCCACcacacatgga                                                                                                                                                                                                          >  1:512837/1‑232 (MQ=255)
                                      aaTGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCACGCCAGCGTATTGGCGGCACGCATCCGCGACCGAGCTGACCGAGCGGAAGGGAAAACCCGGTTGGGTGCGACTGGGGCGCCGGCCACACGGGCa                                                                                                                                                                                                          >  2:269630/1‑238 (MQ=255)
                                         gTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGCAATCATCGGCAAAGGGGCCGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAAAGGCacta                                                                                                                                                                                                       >  1:856099/1‑238 (MQ=255)
                                           tCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAAACATCGGCAAAGGAGCTGCCCGGGCGTAGGGAAAAGCCGGGTGGGTTGGGGCCGTGGCACCGGACAAAATgggccaaaag                                                                                                                                                                                                    >  1:1097077/1‑231 (MQ=255)
                                            cAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCCGCAACGGAGCAGCCAGAGCGGAAGGAATAGCCGGTTCGGTTTCGACTGGGGCCCCGAACACAATGGGACCACTg                                                                                                                                                                                                    >  1:1323707/1‑238 (MQ=255)
                                             aGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAAACATCGGCAACGGAGCTGCCCGCGCGGAAGGAATAGCCGGTTCGGGTGGGGCTGTGGCACCGACCACAAAGGCACTACCgg                                                                                                                                                                                                   >  2:278101/1‑238 (MQ=255)
                                               cgcgAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCCGCCACGGAGCGGCCAGAGCGTAAGGAAGAACCGGGTCGGGTGGGACCGTGGGACCGGCCCCCAAGGCCCTCCTGGCa                                                                                                                                                                                                 >  2:754034/1‑238 (MQ=255)
                                               cgcgAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGGAATCATCGGCCACGGAGCTGCCAGAGCGGAAGGAAAAACCGGGTCGGGTGCGGCCGGGGGCCCGCACACAAAGGCACCACTGGCa                                                                                                                                                                                                 >  2:1155817/1‑238 (MQ=255)
                                                gcaagccaaaTCGCCAGCCCCAGAAAAAACGCACCCCCCGGCGACAACAGGCCGATCAGCGGGTCAAAGATCCCCAAACCAGGTGCCAGCCGTTTGATTCCGAAAGCGACCACCAGCCCGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAg                                                                                                                                                                                                <  2:1300141/230‑1 (MQ=255)
                                                    ggCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAAAGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGGGACTGTGGCACCGACCACACAGGGACTACTgggcgccact                                                                                                                                                                                           >  1:1168749/1‑231 (MQ=255)
                                                           cGCCAGCACCAGATAAAACGAACCCGCCGGCGACATCAGGCCGACCAGCGGGTCAAAGACACCCAAACCAGGTGCCAGACTTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCt                                                                                                                                                                                    <  2:1054994/239‑1 (MQ=255)
                                                             ccacccccAGATAAAACGACCCCACCGGGGACATCGGGCCGACCGGGGGGTCAAAGATCCCCAACCCAGTTGCCAGACGTTTGATCCCAAAAGCGACCCCCACCACGCCGTTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGt                                                                                                                                                                                  <  2:512837/233‑1 (MQ=255)
                                                              cAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGGTCGGTTGCGACTGTGGCACCGACCACAAAGGCACCACTGGCAGCCTCTTTTTCCTGTc                                                                                                                                                                                 >  1:280667/1‑239 (MQ=255)
                                                               aGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCCGAGCGGAAGGAATAGCCGGTTTGGGTGCGACTGGGGCACCGACCCCAATGGGACTACACGCAGCCTCTTTTTcccttcg                                                                                                                                                                                >  2:175343/1‑234 (MQ=255)
                                                                gCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATTCTCCGCCACGGAGCTGCCAGAGCGGAAAGGATAACCGGGTCGGTTGCGACTGGGGCAACGGCCACAAAGGCACTACTGGCAACCTCTTTTTCCTGTCgg                                                                                                                                                                               >  2:594096/1‑239 (MQ=255)
                                                                  acaaGAAAAAACAAACCCACCGGCGACATCAGGCCAATCGGCGGGTAAAAGATCCCAAACCAGGGTGCCAGACGTTTGATCCCGAAAGCGATCCCCAGCACGCGGGTAAGGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCg                                                                                                                                                                                <  2:569369/233‑1 (MQ=255)
                                                                   ccAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTa                                                                                                                                                                            <  1:603688/239‑1 (MQ=255)
                                                                      gATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCCACGGAGCTGCCAGAGCGGAAGGAAAAACCGGTTCGGTTGCGACTGTGGCACCGCCCACCATGGCACTAATGGCAGCCTCCTTTTTCTGTCGGGTAAAc                                                                                                                                                                         >  1:670224/1‑239 (MQ=255)
                                                                      caaaaaaCGAACCCCCCGGAGACAGCAGGCCGCACAGGGGGTAAAAGATCCCCAACCCGGGTGCCAGCCGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAAc                                                                                                                                                                         <  2:1062046/236‑1 (MQ=255)
                                                                        tAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAAc                                                                                                                                                                         <  1:618769/237‑1 (MQ=255)
                                                                               aaCCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGGAATCATCGGCCACGGGGCTGCCAGGGCGGAAGGAATAGCCGGTTCGGTTGCGACTGTGGCGCCGACCACCAAGGCACCACTCGCCGCCCCCTTTTTCTGTCTGGTAAAAAGTTTCGgc                                                                                                                                                                >  2:314843/1‑239 (MQ=255)
                                                                                aCCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGgc                                                                                                                                                                >  1:368111/1‑238 (MQ=255)
                                                                                     ccGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCCAAGGAGCTGCCCGAGCGTAAGGAAAAACCGGGTCGGTTGCGACTGTGGGACCGACCACAATGGGAACACTGGCGGCCTCTTTTTTCTGTCGGGGAAACAGATTTGGCGCAATg                                                                                                                                                          >  2:529782/1‑239 (MQ=255)
                                                                                       ggCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGaa                                                                                                                                                        <  1:86710/239‑1 (MQ=255)
                                                                                          gACATCAGGCCGACCAGCGGGCCAAAGATACCCAACCCGGGTGCCAGAGGTTTGATCCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACg                                                                                                                                                     <  2:724574/238‑1 (MQ=255)
                                                                                            cATCAGGCCGATCAGCGGGTAAAAGATCCCCAACCCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGAc                                                                                                                                                                                                                    <  2:397133/174‑1 (MQ=255)
                                                                                                     caccggggggTAAAAGATCCCAAAGCCGGTCTCAGGTCTTTTATTCCGAAAAGCGCCCCCAGCCCCGCCGTTAATCCCGCCCGCCGGCGTATTGGCGGTAATCATGGGCAACGGAGCTCCCAGAGCGTAAGAAATAGCCGGTTCGTTTGCGACTGTGGCCCCGACCACAATGGCACTCCTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTGGGGCGCAATGAACGTCGCAAGACccg                                                                                                                                         <  2:1291617/232‑1 (MQ=255)
                                                                                                      aTCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCCAACGAGCTGCCAGAGCGGAAGGAATAGCCGGGTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGGCCCTTTTTTCCGTCGGGGAAACAGGTTTGGCGCAATGAACATCCCAAGAaccc                                                                                                                                         >  2:742344/1‑237 (MQ=255)
                                                                                                      aTCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCCGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGGTCGGTTGCGACTGTGGCACCGACCACAATGGCa                                                                                                                                                                                                          >  1:397133/1‑174 (MQ=255)
                                                                                                           cGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTTGGCGCAATGAACGTCGCAAAACCCGCcgcg                                                                                                                                    >  1:938951/1‑238 (MQ=255)
                                                                                                             ggTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATTCTCCGCAACGGAGCCGCCCAAGCGTAAGGGAAAGACGGTTCCGCTGCGACCGCGGCACCGACCCCAAATGCACCACCGGCAGCCTCTTTCTTCTGTCCGGTAAAACGATCGGGCGCCATGACCGTCGCACTACCcgccgcccac                                                                                                                                 >  2:1060536/1‑235 (MQ=255)
                                                                                                                        cccAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGGAATCATCGGCCACGGAGCTGCCCGAGCGGAAGGAAAAGCCGGTTCGGGTGCGACTGTGGCACCGACCACAATGGCACTACCTGCAGCCCCTTTTTTCTGTTCGGTGAACAGGTTTGGCGCAAAGAACGGCGCCAGACCCCCCCCGTCgggggggctga                                                                                                                       >  2:1039856/1‑233 (MQ=255)
                                                                                                                             aCCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGGAAACAGGTTTGGCGCCATGAACGTCGCCAGACCCGCCGCGACGGGCGGCATCAGcgcg                                                                                                                  >  1:23870/1‑238 (MQ=255)
                                                                                                                             aCCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGGAAGGAATAGGCGGTTCGGGTGCGACTGTGGGACCGACCACAATGGCACTACTGGCAGCCTCCTTTTTCTGTCGGGGAAACAAGTTTGGCGCAAATAACCGCGCCAAGCCCGGCGCGACGGGCGGCagacggcgg                                                                                                                  >  2:65039/1‑230 (MQ=255)
                                                                                                                                               tGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCCGCAACGGAGCTGCCAGAGCGTAAGGAAAAACCGGGTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCCCTTTTTTCTGTCGGGGAAACAGTTTGGGCGCAACGAACGGCGCCAAACCCGCCCCGGCGGGCGGCCTCAGGCCGACAAAGcccgcgcagggg                                                                                                >  2:1145576/1‑229 (MQ=255)
                                                                                                                                                      gAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCCGCAACGGAGCTGCCAGAGCGGAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGGAAACAGTTTGGGCGCAATGAACGTCGCCAGACCCCCCCCGACGGGCGGCCTCCGCGCGACCACGCCGCCCATGGCGTACCac                                                                                          >  2:68690/1‑237 (MQ=255)
                                                                                                                                                        aagggcacccccAGCCCGGCGGAAGGGCCGCCCGCCAGCGTATTGGGGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTc                                                                                        <  1:68690/229‑1 (MQ=255)
                                                                                                                                                        aaaGCGACCCCCGGCCCGCGGGTAGGCCCGCCCGCCAGCGTATGGGCGGTAATCACCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCGGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTGGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTc                                                                                        <  1:622097/239‑1 (MQ=255)
                                                                                                                                                             gATCCCGAGCACGCGGGTAATGCCCCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCCCCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATaaa                                                                                   <  1:742344/239‑1 (MQ=255)
                                                                                                                                                                      cACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCCGCCAGAGCGGAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGGTTGGGCGCAATGAAAGGCGCAAAACCCGCCCCGGCGGGCGGCAACACGGCGCACACCCCGCCCACGGCGGACCAGCCATAAAAGAGttttt                                                                          >  2:1257764/1‑239 (MQ=255)
                                                                                                                                                                        acgcggtaaaggccgccCGCCAGCGTTTTGGGGGTAAGCAGCGGCAACGGGGCTGCCAGGGGTTAAGAAGTAGCGGGTGCGGTTGCGAGTGGGGCACCGCCCAAAATGGCCATACTGGCAGCCTCTTTTCCCTGTGGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGCCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTc                                                                         <  2:1142804/228‑1 (MQ=255)
                                                                                                                                                                            cGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTTCTGTCGGGGAAAAAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCcagc                                                                     >  1:342983/1‑239 (MQ=255)
                                                                                                                                                                                      gcccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGGAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGGAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGGACCAGTCATAAATGTGTTTTTTCAGCAGTGAGAAGc                                                           >  1:349611/1‑239 (MQ=255)
                                                                                                                                                                                       cccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCCCAATGGCACTACTGGCAGCCTCTTTTTCCTTTCGGGGAAACCGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCCTAAATGTGTTTTTCCAGCAGTGCGAAGCa                                                          >  1:354151/1‑239 (MQ=255)
                                                                                                                                                                                        ccgcgggggTATGGGGGGAAACCGTCGGCAACGGGCCTGCCAGGGCGTAAGGAATGGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCCCTACTGGCAGCCTTTTTTTCCTGTCGGGAAAACAGTTTGGGCGAAATGAACGCCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAGGCag                                                         <  2:670224/231‑1 (MQ=255)
                                                                                                                                                                                                 ttttGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGAAAGAAATAGCCGGTTCGTTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCa                                                 <  1:1232954/236‑1 (MQ=255)
                                                                                                                                                                                                  agtGGCGGAAATCACCGGCAAGGGGGCTGCCAGGGCGGAAGAAAAGGCGGGTTCGGTGGCGACTGTGCCCCCACCCAAAAGGGCCCTACGGGCGCCTTCTTTTCCCTGTGGGAAAAACAGTTGGGCCGAAATGAACGTCGCAAGCCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAgg                                               <  2:288620/237‑1 (MQ=255)
                                                                                                                                                                                                  aTTGGCGGTAATCATCGGCAACGGGGCTGCCAGAGCGTAAGGAATGGCGGGTTCGGTTGCGACTGTGGCCCCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAgg                                               <  1:754034/239‑1 (MQ=255)
                                                                                                                                                                                                  aTTGGCGGTAATCACCGCCAACGGAGCGGCCAGAGCGTAAGGAATAGCCGGTTCGGTGGCGACGGTGGCACCGACCACAATGGCACTACTGGCGCCCTCTTTTTCCTGTCGGATAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAgg                                               <  1:314843/239‑1 (MQ=255)
                                                                                                                                                                                                          gaaTCACCGGCAACGGAGCGGCAAGAGCGGAGAGAAAAGCCGGTTCGGTGGCGACTGTGGCACCGACCACAAGGCCACTTGTGGCAGCCTCTTTTCCCGGTCGGGTAAACAGTTTGGGCGAAAGAACGGCCGCAAGCCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGtttt                                        <  2:703682/237‑1 (MQ=255)
                                                                                                                                                                                                              cctcGGCACGGGAGCTGCAAGAGGGTAAGGAAGAGCGGGTGCGGTGGCGACGGCGGCCCCGACCACAAGGGCACTACTGGCAGCCTCTTTTTCCTGCCGGTAAAACAGTTTGGGCGCAAGGACCGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTTTTATTGa                                   <  2:1304864/237‑1 (MQ=255)
                                                                                                                                                                                                                tcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTTGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATAGCGTCCCCGTCATCAATGTGTTTTTCCAGCAGTGTGAAGCCGAGGCACCAGGCGGATTTCTTGCTg                                >  2:603339/1‑239 (MQ=255)
                                                                                                                                                                                                                tcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTTTTATTGATg                                >  1:1194718/1‑239 (MQ=255)
                                                                                                                                                                                                                          aGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGCTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTTTTCTTGATTGGGgccgcc                        >  2:528106/1‑238 (MQ=255)
                                                                                                                                                                                                                              gCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGCCGGGCGGCATCAGCGCGACAACGCCGAAGATGGCGCACCAGTCATAAATGTGTTTTTCAAGCAGTGAGAAGCAGAAGAACCAGGCGGGTTTATTGTTGGGGCCCGCCATATc                  >  2:1043736/1‑239 (MQ=255)
                                                                                                                                                                                                                                     cgaaaGGAAAAGCCGGTTCGGTTGCGACGGGGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGAAAGCCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTTTTATTGATGGGGCCGCCCATATCGAACGc            <  1:1012662/235‑1 (MQ=255)
                                                                                                                                                                                                                                         aGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAgaagaa                                                    >  1:607227/1‑195 (MQ=255)
                                                                                                                                                                                                                                           ggAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTTTTATTGATGGGGCCGCCCATATCGAACGCCAGCa       >  2:1094124/1‑238 (MQ=255)
                                                                                                                                                                                                                                              aTAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTTTTATTGATGGGGCCGGCCCTAACGGACCCCCGccatcc   >  2:1007969/1‑234 (MQ=255)
                                                                                                                                                                                                                                               tAGCCGGTTCGGTTGGGATTGTGGCACCGACCAAAATGGCACTACTGGCCGCCTCTTTTTCCTGGCGGGAAAACAGTTTGGGCGCAATGAACGCCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTTTTATTGATGGGGCCGCCCATATCGAACGCCAGCATCGCa  <  1:295941/239‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
CGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTTTTATTGATGGGGCCGCCCATATCGAACGCCAGCATCGCA  >  NC_000913/4090958‑4091431

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