Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,091,186 T→G 24.2% Y394S (TAC→TCC)  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,1860TG24.2% 128.8 / 4.4 73Y394S (TAC→TCC) frvBputative PTS enzyme, IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base T (21/28);  new base G (19/2);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 1.90e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTTTTATTGATGGGGCCGCCCATATC  >  NC_000913/4090949‑4091415
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                        
ttCGCTTTGGGGAGTAACCACAGTCCGTTCAGGACGAGAATGAATGAAAGGTTCGGCGCGAGGCCAATGCCCATGCCCAGAAAAACCGAACCCACCGGCGCCATCAGGGCGATCAGCGGGTCAAAGATCCCCAAACCAGGTCCCGGCCGTTTGATTCCGAAAGCGTTCACCAGGCCGCCTGTAAGGCCGCCTGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGt                                                                                                                                                                                                                                          <  2:158451/238‑1 (MQ=255)
        gcggcgTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGGAAAGGGGCCGCCAGAGCGAAAGGaaaa                                                                                                                                                                                                                                 >  1:1075013/1‑236 (MQ=255)
         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)
                           ttCGGGACGAAAAGGAATGAAATTTCCAGCGCGAGGCCAATCGCCAGCCCCAGAAAAAACGAACCCCCCGGCGACATCAGGCCGATCAGCGGGTCAAAGATCCCCAACCCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGt                                                                                                                                                                                                                                          <  2:153950/211‑1 (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)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                        
TTCGCTTTGCGGCGTAACCACAATCCTTTCAGGACGATAATGAATGAAATGTTCAGCGCGAGGCCAATCGCCAGCACCAGATAAAACGAACCCACCGGCGACATCAGGCCGATCAGCGGGTCAAAGATACCCAAACCAGGTGCCAGACGTTTGATTCCGAAAGCGATCACCAGCACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTTTTATTGATGGGGCCGCCCATATC  >  NC_000913/4090949‑4091415

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: