Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,259,696 T→C 5.9% intergenic (‑400/+23) psuK ← / ← fruA putative pseudouridine kinase/fructose‑specific PTS multiphosphoryl transfer protein FruA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,259,6960TC 5.9% 369.7 / 9.1 119intergenic (‑400/+23)psuK/fruAputative pseudouridine kinase/fructose‑specific PTS multiphosphoryl transfer protein FruA
Reads supporting (aligned to +/- strand):  ref base T (48/64);  new base C (4/3);  total (52/67)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 122 total aligned reads displayed.

GAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAGGGATC  >  NC_000913/2259547‑2259845
                                                                                                                                                     |                                                                                                                                                      
gAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAAttt                                                                                                                                                        >  2:625704/1‑150 (MQ=255)
    ataatCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGc                                                                                                                                                   <  1:23102/150‑1 (MQ=255)
     taatCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGCAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGcc                                                                                                                                                  <  2:1115100/150‑1 (MQ=255)
     taatCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGcc                                                                                                                                                  <  1:1536655/150‑1 (MQ=255)
            cccTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAac                                                                                                                                           <  2:1416585/150‑1 (MQ=255)
            cccTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAac                                                                                                                                           <  1:1435042/150‑1 (MQ=255)
             ccTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTa                                                                                                                                             <  1:190783/147‑1 (MQ=255)
             ccTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAaca                                                                                                                                          >  1:614008/1‑150 (MQ=255)
                tGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAacacac                                                                                                                                       >  2:1206048/1‑150 (MQ=255)
                 gAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATNNNATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACAcc                                                                                                                                      <  2:1448987/150‑1 (MQ=255)
                  aGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCt                                                                                                                                     >  2:1688522/1‑150 (MQ=255)
                   ggAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAacatttgat                                                                                                                                    <  1:602797/150‑7 (MQ=38)
                      aaTAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAacatttgat                                                                                                                                    <  1:1880586/147‑7 (MQ=38)
                      aaTAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGGCAGAAATTTCCTGTCCTGTAacatttgat                                                                                                                                    >  2:1880586/1‑141 (MQ=38)
                      aaTAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTtta                                                                                                                                 >  2:1941338/1‑150 (MQ=255)
                        tAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttatt                                                                                                                               <  2:381765/150‑1 (MQ=255)
                        tAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttatt                                                                                                                               <  1:117480/150‑1 (MQ=255)
                          aCCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTAc                                                                                                                             >  2:477565/1‑150 (MQ=255)
                            cTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACgc                                                                                                                           <  1:271069/150‑1 (MQ=255)
                               aGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACgctgc                                                                                                                        <  1:743544/150‑1 (MQ=255)
                                    ttGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCg                                                                                                                   <  2:1454734/150‑1 (MQ=255)
                                      gatgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCt                                                                                                                 <  1:520837/150‑1 (MQ=255)
                                      gatgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCt                                                                                                                 <  1:481597/150‑1 (MQ=255)
                                       atgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTa                                                                                                                >  1:1787635/1‑150 (MQ=255)
                                        tgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTAc                                                                                                               >  1:1468153/1‑150 (MQ=255)
                                        tgaACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTAc                                                                                                               >  2:130018/1‑150 (MQ=255)
                                          aaCGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTg                                                                                                             <  1:1031388/150‑1 (MQ=255)
                                             tttttAAGTTTATACTTAAATATAAAACAATTTAAATTAATAAAATTACATTCAAAAACCGGTCATTCTCTTTGCGATCATCAAATGTTACAGGACAGTAAATTTCTGCCCTGTAACACACATTTTATTACGCTGCTTTCGCTACTGCGt                                                                                                          <  1:807360/149‑1 (MQ=255)
                                             gTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGt                                                                                                          <  2:1371740/150‑1 (MQ=255)
                                                tGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTccc                                                                                                       >  2:1374016/1‑149 (MQ=255)
                                                tGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCa                                                                                                       >  2:1037760/1‑150 (MQ=255)
                                                   aGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACtt                                                                                                    >  2:241452/1‑150 (MQ=255)
                                                    gTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTc                                                                                                   >  1:139705/1‑150 (MQ=255)
                                                     tttCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTcc                                                                                                  <  2:295644/150‑1 (MQ=255)
                                                      ttCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCg                                                                                                 <  2:1610945/150‑1 (MQ=255)
                                                        cTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGa                                                                                               <  2:356808/150‑1 (MQ=255)
                                                        cTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCACCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGa                                                                                               >  2:1536238/1‑150 (MQ=255)
                                                         tACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGAc                                                                                              <  2:1904265/150‑1 (MQ=255)
                                                         tACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGAc                                                                                              >  1:1295837/1‑150 (MQ=255)
                                                         tACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGAc                                                                                              <  1:625704/150‑1 (MQ=255)
                                                           cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                                            <  1:428272/150‑1 (MQ=255)
                                                           cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                                            >  2:135392/1‑150 (MQ=255)
                                                           cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                                            <  1:293231/150‑1 (MQ=255)
                                                           cTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCAGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGt                                                                                            >  2:1083785/1‑150 (MQ=255)
                                                            tGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGtt                                                                                           <  1:1037760/150‑1 (MQ=255)
                                                              aaaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGGTTTCGCTACTGCGTCCACTTCCGGACGtttg                                                                                         <  1:1941338/150‑2 (MQ=255)
                                                               aaGATAAACCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCa                                                                                        <  1:961312/150‑1 (MQ=255)
                                                               aaGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCa                                                                                        >  1:1950339/1‑150 (MQ=255)
                                                                aGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAg                                                                                       <  1:1648049/150‑1 (MQ=255)
                                                                 gATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttatta                                                                                                                              >  1:148693/1‑110 (MQ=255)
                                                                 gATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTttatta                                                                                                                              <  2:148693/110‑1 (MQ=255)
                                                                     aaaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAgg                                                                                  >  2:1680615/1‑150 (MQ=255)
                                                                      aaCCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGc                                                                                 <  1:1206048/150‑1 (MQ=255)
                                                                              aaaTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACCCCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCa                                                                         >  2:1612900/1‑150 (MQ=255)
                                                                               aaTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaa                                                                        >  1:609046/1‑150 (MQ=255)
                                                                                aTTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCaaa                                                                       <  1:39272/150‑1 (MQ=255)
                                                                                 ttCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCCAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc                                                                      <  1:314982/150‑1 (MQ=255)
                                                                                 ttCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc                                                                      <  1:1871625/150‑1 (MQ=255)
                                                                                 ttCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAAc                                                                      <  1:63655/150‑1 (MQ=255)
                                                                                   cATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACgc                                                                                                                           <  1:1646856/95‑1 (MQ=255)
                                                                                   cATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCATGTAACACACCTTTTATTACgc                                                                                                                           >  2:1646856/1‑95 (MQ=255)
                                                                                    aTAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCg                                                                   >  1:1902865/1‑150 (MQ=255)
                                                                                     tAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGc                                                                  >  2:1941401/1‑150 (MQ=255)
                                                                                      aaaaTTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGTCCAAACCCGcc                                                                 <  1:1232696/150‑1 (MQ=255)
                                                                                      aaaaTTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGcc                                                                 <  1:1673454/150‑1 (MQ=255)
                                                                                       aaaTTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCAGGACGTTTCAGGAATGCATAGGCCAAACACGcca                                                                >  2:900348/1‑150 (MQ=255)
                                                                                               ttCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTa                                                        >  1:338990/1‑150 (MQ=255)
                                                                                                  acaAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCg                                                     >  2:194890/1‑150 (MQ=255)
                                                                                                   caAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCgg                                                    <  2:793821/150‑1 (MQ=255)
                                                                                                    aaaCCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGTATAGGCCAAACCCGGCACCAGCGTACCggg                                                   <  2:867201/150‑2 (MQ=255)
                                                                                                    aaaCCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGc                                                   <  1:900348/150‑1 (MQ=255)
                                                                                                    aaaCCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGc                                                   <  1:898053/150‑1 (MQ=255)
                                                                                                       ccGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCaat                                                <  1:494734/150‑1 (MQ=255)
                                                                                                             aTTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGc                                          <  1:1770693/150‑1 (MQ=255)
                                                                                                              ttCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCt                                         <  2:1950339/150‑1 (MQ=255)
                                                                                                                   tttGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTAGTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAg                                    >  1:476325/1‑150 (MQ=255)
                                                                                                                     tGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGt                                  >  1:1072056/1‑150 (MQ=255)
                                                                                                                     tGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGt                                  <  1:1792392/150‑1 (MQ=255)
                                                                                                                     tGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCCGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGt                                  <  1:1083785/150‑1 (MQ=255)
                                                                                                                       cgcTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTaa                                >  2:1303518/1‑150 (MQ=255)
                                                                                                                       cgcTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTaa                                >  2:429432/1‑150 (MQ=255)
                                                                                                                        gcTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAAc                               <  1:1941401/150‑1 (MQ=255)
                                                                                                                         cTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAAcc                              >  2:894061/1‑150 (MQ=255)
                                                                                                                         cTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAAcc                              <  2:284983/150‑1 (MQ=255)
                                                                                                                          tcatcaAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAAccc                             <  1:130018/150‑1 (MQ=255)
                                                                                                                          tcatcaAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGACAATAATTGCTACCAGGTAAccc                             <  2:1308368/150‑1 (MQ=255)
                                                                                                                               aaaTGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATAc                        >  2:337542/1‑150 (MQ=255)
                                                                                                                                  tGTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCCCCAGCGTACAGGCAATAATTGCTACCAGGTAGCCCAATACCgg                     >  1:1905822/1‑150 (MQ=255)
                                                                                                                                   gTTACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGc                    >  1:1990072/1‑150 (MQ=255)
                                                                                                                                    ttACAGGGCAGAAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCCGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCg                   >  1:869340/1‑150 (MQ=255)
                                                                                                                                    ttACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATATTTGCTACCAGGTAACCCAATACCGGCg                   <  1:904760/150‑1 (MQ=255)
                                                                                                                                    ttACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCg                   <  2:1787635/150‑1 (MQ=255)
                                                                                                                                         ggacaggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGCCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATa              >  2:1438898/1‑150 (MQ=255)
                                                                                                                                         ggacaggaAATTTCTGCCCTGTAAAACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATa              >  2:1982221/1‑150 (MQ=255)
                                                                                                                                            caggaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAgcg           >  1:472771/1‑150 (MQ=255)
                                                                                                                                              tgaAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGcc         >  1:1035719/2‑150 (MQ=255)
                                                                                                                                               gaAATTTCCTGTCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGcc         >  2:1729745/1‑150 (MQ=255)
                                                                                                                                                 aaTTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAgg      <  2:139705/150‑1 (MQ=255)
                                                                                                                                                   tttCTGCCCTGTAACACACCTTTTATTATGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAGGGa    <  1:1859657/150‑1 (MQ=255)
                                                                                                                                                     tCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAGGGATc  <  1:1594820/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCTGCCCTGTAACACACCTTTTATTACGCTGCTTTCGCTACTGCGTCCACTTCCGGACGTTTCAGGAAGGCATAGGCCAAACCCGCCACCAGCGTACCGGCAATAATTGCTACCAGGTAACCCAATACCGGCGTAATAGCGCCAGGGATC  >  NC_000913/2259547‑2259845

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