Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 114,507 A→T 5.1% intergenic (+20/+15) guaC → / ← hofC GMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913114,5070AT 5.1% 397.0 / 7.6 144intergenic (+20/+15)guaC/hofCGMP reductase/assembly protein in type IV pilin biogenesis, transmembrane protein
Reads supporting (aligned to +/- strand):  ref base A (65/68);  new base T (5/3);  total (73/71)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

Only 100 of 146 total aligned reads displayed.

GAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAGCGC  >  NC_000913/114360‑114656
                                                                                                                                                   |                                                                                                                                                     
gagaTATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGc                                                                                                                                                     >  2:1730006/1‑150 (MQ=255)
gagaTATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGc                                                                                                                                                     <  1:539511/150‑1 (MQ=255)
 agaTATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGTGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCa                                                                                                                                                    <  1:337657/150‑1 (MQ=255)
 agaTATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGTGCa                                                                                                                                                    >  1:1284963/1‑150 (MQ=255)
   atatTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGCGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAac                                                                                                                                                  <  2:167084/150‑1 (MQ=255)
     atTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAacac                                                                                                                                                >  2:1858183/1‑150 (MQ=255)
       tttGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACATGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGc                                                                                                                                              >  2:1702218/1‑150 (MQ=255)
         tGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCCACAAACTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCa                                                                                                                                            >  1:453710/1‑150 (MQ=255)
         tGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATATTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCa                                                                                                                                            <  2:865353/150‑1 (MQ=255)
          gggCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGTTGCTCCACGCCAg                                                                                                                                           <  2:1095014/150‑1 (MQ=255)
                ccTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTat                                                                                                                                     >  2:1951415/1‑150 (MQ=255)
                ccTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTat                                                                                                                                     >  2:1918089/1‑150 (MQ=255)
                 cTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatc                                                                                                                                    <  2:1767899/150‑1 (MQ=255)
                  tGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcc                                                                                                                                   >  1:1848930/1‑150 (MQ=255)
                   gCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatccc                                                                                                                                  >  2:1797492/1‑150 (MQ=255)
                     gTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGAGGAGCAACACGCCAGCGTTatccaat                                                                                                                                >  2:1274432/1‑150 (MQ=255)
                      ttCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccatc                                                                                                                               <  1:1971532/150‑1 (MQ=255)
                       tCAGCTTGTACATACGTTGGGGCCTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTatcccatcc                                                                                                                              <  1:1129810/150‑1 (MQ=255)
                           cTTGTACATACGTTGGGGCTTCACGCCTGAAAGCGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGATATCCCATCCCACt                                                                                                                          >  1:1812495/1‑150 (MQ=255)
                            ttGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTc                                                                                                                         <  2:1686845/150‑1 (MQ=255)
                               tACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCGCGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatc                                                                                                                      <  2:2085779/150‑1 (MQ=255)
                                     cGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGTGTGGAGCAACACGCCAGCGTTATCCCATCCCACTcatcgcatcg                                                                                                                >  2:990784/1‑150 (MQ=255)
                                     cGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCCTAATCCCATCCCACTcctcgcatcc                                                                                                                >  1:99008/1‑149 (MQ=255)
                                       ttGGGGCTTCACGGATGAAAGATCTGACCAAGCGCACCACGTTTATTCGAGTGCAGGAACAAGAAACCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGGGTAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGcc                                                                                                              <  2:824381/150‑1 (MQ=255)
                                       ttGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGcc                                                                                                              <  2:941447/150‑1 (MQ=255)
                                         ggggCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTa                                                                                                            >  2:1045900/1‑150 (MQ=255)
                                          gggCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAACACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGAGTGGAGCAACACGCCAGCGTTAACCCATCCAACTCATCGCATCGCCTaa                                                                                                           >  1:1357518/1‑150 (MQ=255)
                                           ggCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACCACCTGTAATCTCCCAACGCTGGCGGGGCGCAACACGCCAGCGTTATCCCCTCCCCCTCATCGCATCGCCGaaa                                                                                                          >  2:310248/1‑150 (MQ=255)
                                            gCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAAt                                                                                                         <  2:1974992/150‑1 (MQ=255)
                                             gttCACGCCTGAAAGAGCTGACCAAGCGCACCTCGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATg                                                                                                        <  2:214848/149‑1 (MQ=255)
                                             cTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATg                                                                                                        <  2:335144/150‑1 (MQ=255)
                                                  cGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGaaaa                                                                                                   >  2:953673/1‑150 (MQ=255)
                                                     cTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGGGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTg                                                                                                <  2:1226049/150‑1 (MQ=255)
                                                       gAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCCGCGTTATCCCAGCCCACTCATCGCATCGCCTAAATGGAAAATTGGc                                                                                              >  1:655424/1‑150 (MQ=255)
                                                       gAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGc                                                                                              <  2:1296954/150‑1 (MQ=255)
                                                         aaGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAg                                                                                            <  2:2102753/150‑1 (MQ=255)
                                                         aaGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAg                                                                                            <  2:494779/150‑1 (MQ=255)
                                                         aaGAGCTGACCAAGCGCACCACGGTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAg                                                                                            <  2:725624/150‑1 (MQ=255)
                                                          agagCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGa                                                                                           <  2:1848930/150‑1 (MQ=255)
                                                             gCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCCACAACCTGTAATCTCCCAACGCTGGCGTGGAGCACCACGCCCGCGTTATCCCATCCCACGCATCGCAACGCCGAAATGGAAAATTGGCAGAGAc                                                                                        >  1:171968/1‑150 (MQ=255)
                                                                gACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACAtt                                                                                     >  1:1630673/1‑150 (MQ=255)
                                                                    aaGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGccc                                                                                 >  1:1721553/1‑149 (MQ=255)
                                                                    aaGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGcca                                                                                 >  2:770742/1‑150 (MQ=255)
                                                                         caccacGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGGCCGCGTTATCCGCTACCACTCATCGCATCGCCTAAATGGATAATTGGCAGCTGCATTGccgccacc                                                                            >  2:1509927/1‑150 (MQ=255)
                                                                          accacGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCATCAACAACCTGAAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCAACCCACTCATCGCATCTCCTAAATTTAATATTGGCAGATACATTGccaccacca                                                                           <  1:1562588/150‑1 (MQ=255)
                                                                           ccacGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAg                                                                          <  1:1403264/150‑1 (MQ=255)
                                                                           ccacGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCACCCAACGCTGGCGTGGAGCAACACCCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAg                                                                          <  1:79028/150‑1 (MQ=255)
                                                                                tttATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCCCTCATCGCATCGCCTAAATGGAAAATTGGCAGCTACATTGCCCCCACCAGCGtgc                                                                     >  2:1631610/1‑148 (MQ=255)
                                                                                   aTTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaa                                                                  >  2:540926/1‑150 (MQ=255)
                                                                                    ttCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTATATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaat                                                                 >  2:676219/1‑150 (MQ=255)
                                                                                     tCGTGTGCAGGAACAAGAAAACCGCATCATCAACAACCAGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCaata                                                                <  2:203090/150‑1 (MQ=255)
                                                                                        tgtgCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAAtt                                                             <  1:1961146/150‑1 (MQ=255)
                                                                                        tgtgCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAAtt                                                             <  1:602564/150‑1 (MQ=255)
                                                                                          tgCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATtcc                                                           >  2:2040313/1‑150 (MQ=255)
                                                                                           gCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCCCCACCAGCGTACCAATAATtcct                                                          >  1:2074009/1‑150 (MQ=255)
                                                                                              gggaCAAGAACACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCCTCCCCCTCAGCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTGCCCATAATTCCGccc                                                       >  1:1456493/4‑150 (MQ=255)
                                                                                              ggAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGCGCAACACGCCAGCGTTATCCCATCCCACTCATCGCCTCGCCTAAATGGAAAATTGGCAGATACATTGCCACCGCCCGCGTCCCAATAATTCCGccc                                                       >  1:1265283/1‑150 (MQ=255)
                                                                                              ggAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGGACCAATAATTCCTccc                                                       >  1:2095354/1‑150 (MQ=255)
                                                                                                  cAAGAAAACCGCAACTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACCCGCCAGCGTTATCCCAACCCACTCATCGCATCGCCTAAATGCAAAATAGGCACATACATTCCCACCCCCCGCGTACCAAGTAATCCTCCCGTTa                                                   >  1:225566/1‑150 (MQ=255)
                                                                                                   aaGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTAt                                                  <  2:651875/150‑1 (MQ=255)
                                                                                                    aGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGACACCACCAGCGTACCAATAATTCCTCCCGTTATg                                                 >  1:577173/1‑150 (MQ=255)
                                                                                                       aaaCTGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATc                                              <  1:836555/150‑1 (MQ=255)
                                                                                                        acccGCATCTTCAACAACCTGAAATCTCCCAACGCGGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCa                                             <  2:301236/148‑1 (MQ=255)
                                                                                                        aaCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCa                                             <  1:1045900/150‑1 (MQ=255)
                                                                                                        aaCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCa                                             <  1:1684671/150‑1 (MQ=255)
                                                                                                        aaCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCa                                             <  1:226253/150‑1 (MQ=255)
                                                                                                        aaCCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCa                                             <  1:948911/150‑1 (MQ=255)
                                                                                                          ccGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGc                                           >  2:695478/1‑150 (MQ=255)
                                                                                                          ccGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGc                                           <  2:1407301/150‑1 (MQ=255)
                                                                                                          ccGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGc                                           >  2:980340/1‑150 (MQ=255)
                                                                                                               tCTTCAACACCCTGTAATCTCCCAACCCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCTCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACAAGAGTACCAATAATTCCTCCAGTTATTATCAGCAACGc                                      >  1:1290643/1‑150 (MQ=255)
                                                                                                                cTTCAACAACCTGTAATCTCCCAACTCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAAATGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGcc                                     <  2:38787/150‑1 (MQ=255)
                                                                                                                cTTCAACAACCTGTAATCTCCCAACGCGGGCGTGGCGCAACACGCCAGCGTTATCCCATCCCACTCATGGCCTCGCCTAAATGGAAAAGTGGCAGATACATTGCCACCCCCCGCGGCCCAATAATTCCTCCCGGGAGGATCAGCACCGcc                                     >  1:1489557/1‑150 (MQ=255)
                                                                                                                  tCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCATAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCgg                                   <  2:1258721/150‑1 (MQ=255)
                                                                                                                   caacaaCCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCCCCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGt                                  >  1:922422/1‑150 (MQ=255)
                                                                                                                   caacaaCCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATACTTCCTCCCGTTATGATCAGCAACGCCGGt                                  >  2:436753/1‑150 (MQ=255)
                                                                                                                   caacaaCCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGt                                  <  2:747448/150‑1 (MQ=255)
                                                                                                                    aacaacCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGAAACATTGCCACAACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGtt                                 >  1:1274761/1‑150 (MQ=255)
                                                                                                                      caacCTGTAATCTCCCAACGCTGGCGTGTTGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTcc                               <  2:1520055/150‑1 (MQ=255)
                                                                                                                      caacCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTcc                               <  2:1039152/150‑1 (MQ=255)
                                                                                                                      caacCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTcc                               <  2:1124479/150‑1 (MQ=255)
                                                                                                                      caacCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTcc                               <  2:1478910/150‑1 (MQ=255)
                                                                                                                        acattaaaTCTCCCAACGCTGGCGTGTAGCAACACTCCAGCTTTATCCCATCCCACTCATCGCATCGCCTAAATTGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAg                             <  1:333390/144‑1 (MQ=255)
                                                                                                                              aaTCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGc                       >  1:1228817/1‑150 (MQ=255)
                                                                                                                                tctcCCAACGCTGGCGTGGAGCAACACGCAAGCGTTATCCCATCCCACTCATCGCATCGCCTAACGGGAAAATTGGCAGGTACAATGCCCCAACCCGCGTACAAATAATGCCTCCGGTTATGATCAGCAACGCGGGTACAAGAAACGCTg                     >  1:1026343/1‑150 (MQ=255)
                                                                                                                                  tcCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCCCCAGCGTACCAAGAATTCCACCCGGTATGATCAGCAACGACGGTGCCAGTAAGGCTGcc                   >  2:925265/1‑149 (MQ=255)
                                                                                                                                  tcCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCg                   <  2:675485/150‑1 (MQ=255)
                                                                                                                                  tcCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACACATCGCAACGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCg                   >  1:1333836/1‑150 (MQ=255)
                                                                                                                                  tcCCAACGCTGGCGGGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCCTCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCg                   <  2:1218925/150‑1 (MQ=255)
                                                                                                                                      aaCGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGt               <  2:99008/150‑1 (MQ=255)
                                                                                                                                       aCGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAACTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGtt              >  1:515448/1‑150 (MQ=255)
                                                                                                                                        cGCTGGCGTGGATCAACACGCCAGCTTTATCCCATCACACTCATCGCATCGCCTAAATGGAAACTTGGCAGATACATTGCCACCACCAGCGTACAAATAATTCCTCCCTTTATGATCAGCAACGCCGGTTCCAGTAAGTCTGCGAGGTTa             <  1:485242/150‑1 (MQ=255)
                                                                                                                                         gCTGGCGTGTTGCTCCACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTAt            >  1:2060139/1‑150 (MQ=255)
                                                                                                                                          cTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATc           <  1:1256190/150‑1 (MQ=255)
                                                                                                                                            ggtgTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCg         <  1:990784/147‑1 (MQ=255)
                                                                                                                                            ggCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCg         <  2:1609910/150‑1 (MQ=255)
                                                                                                                                              cGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATTCATTGCCCCCACCAGCGTACCAATAATTCCTCCCGTGATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGCGGTTATCCGcc       >  1:2084582/1‑150 (MQ=255)
                                                                                                                                                  gAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAgcg   >  1:1525195/1‑150 (MQ=255)
                                                                                                                                                   aGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCCACGCCGCTTCCAGTAAGGCTGCGAGGTTATCCGCCAgcgc  >  2:852349/1‑150 (MQ=255)
                                                                                                                                                   aGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAgcgc  <  1:2071707/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GAGATATTTTGGGCGGCCTGCGTTCAGCTTGTACATACGTTGGGGCTTCACGCCTGAAAGAGCTGACCAAGCGCACCACGTTTATTCGTGTGCAGGAACAAGAAAACCGCATCTTCAACAACCTGTAATCTCCCAACGCTGGCGTGGAGCAACACGCCAGCGTTATCCCATCCCACTCATCGCATCGCCTAAATGGAAAATTGGCAGATACATTGCCACCACCAGCGTACCAATAATTCCTCCCGTTATGATCAGCAACGCCGGTTCCAGTAAGGCTGCGAGGTTATCCGCCAGCGC  >  NC_000913/114360‑114656

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