Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,628,6250AC39.1% 198.6 / 31.7 116V105G (GTC→GGC) yfgFcyclic‑di‑GMP phosphodiesterase, anaerobic
Reads supporting (aligned to +/- strand):  ref base A (54/16);  new base C (0/45);  total (54/62)
Fisher's exact test for biased strand distribution p-value = 1.86e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

ATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATC  >  NC_000913/2628477‑2628771
                                                                                                                                                       |                                                                                                                                                   
aTTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAcc                                                                                                                                                   >  2:332156/1‑151 (MQ=255)
   ttCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCt                                                                                                                                                 <  1:1456270/150‑1 (MQ=255)
   ttCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCt                                                                                                                                                 <  2:596162/150‑1 (MQ=255)
   ttCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTg                                                                                                                                                >  1:1695843/1‑151 (MQ=255)
         gCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAAc                                                                                                                                           <  2:1252274/149‑1 (MQ=255)
            cGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGccc                                                                                                                                      >  2:335474/1‑151 (MQ=255)
            cGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGccc                                                                                                                                      <  1:255669/151‑1 (MQ=255)
                  tAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAAc                                                                                                                                  >  1:1458865/1‑149 (MQ=255)
                    aTCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACag                                                                                                                                >  1:2079813/1‑149 (MQ=255)
                      ccTACAAATGCAGAAAACTGGAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACTTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACagag                                                                                                                              <  2:901123/149‑1 (MQ=255)
                       caaaaaaggcacaaaaCGGGAAAATTATCAGAAACAGCGTTGCAGAAAACACAATCTCCAAGAATATCCGCTGGAAAATAAAACGGGTACCCCCATGAAAAACGTTGTTCCGCCGGGAAGTAAAGACCCTGTAACCGCCCCAACaga                                                                                                                               <  1:12384/134‑1 (MQ=255)
                       cTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaat                                                                                                                           >  2:2068205/1‑151 (MQ=255)
                        tACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatg                                                                                                                          >  1:11675/1‑151 (MQ=255)
                         aaaaaGGCAGAAAACTGAAAAAGTATCAAAAACAGCGTTGCAGAAAACAAAATCTGCCAAAATATACGCGGGAAAATAAAACGGTTACCCCCATGTAAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGaa                                                                                                                            <  1:1090984/146‑1 (MQ=255)
                          cAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgat                                                                                                                        >  1:1121179/1‑151 (MQ=255)
                           aaaTGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgat                                                                                                                        >  2:532962/1‑150 (MQ=255)
                           aaaTGCAGCAAACTGAAAAAGTATCAAAAACAGCGTTGCAGGAAACACAACCTGCCAGAATATACGCTGGAAAATTAAACGGGTATCACCATGTGAAACGTTGTTCCGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGaatgat                                                                                                                        <  2:1763653/150‑1 (MQ=255)
                           aaaTGCAGAAAACTGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATCCGCGGGGAAATTAAACGGGTATCCCCATGTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTTTAACCGCCCCAACAGAGaatgat                                                                                                                        <  2:927933/150‑1 (MQ=255)
                           aaaTGCAGAAAACGAAAAAAGTATCAGAAACAGCGTGGAAGAAAACAAAATCTGCCAAAATAGCCGGGGGAAAATAAAAGGGGTACCCCCATGTAAACCGTTTTTACGCGGGGGAAAAAACCCCCTATAACCGCCCCAACAGAGaatgata                                                                                                                       <  1:1694804/151‑1 (MQ=255)
                           aaaTGAAGAAAACTGAAAAAGTATGAGAAACAGCGTTGCAGAAAACACAACCTGCAAGAATATCCGCGGGAAAATAAACGGGTTACCCCCATGTAAAACTTTTTTCCGCGGGGGAGAAAACCCCTTTTAACCGCCCCAACAGAGaatgat                                                                                                                        <  1:1069081/150‑1 (MQ=255)
                            aaagcaacaAACGGGAAAAGTAAAAAAAACAGGGTTGCAGGAAAACAAATCCACAAGAAAAACGGCGGGAAAATTAAACGGGTACCCCCAAGAAAAAGGTTTTTACGCGGGGGGGAAAAGCCCCTGTAACCGCCCCAACAGAGaatgataa                                                                                                                      <  1:1167053/148‑1 (MQ=255)
                             aTGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgataat                                                                                                                     >  2:859078/1‑151 (MQ=255)
                                cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgataatgat                                                                                                                  >  2:1247585/1‑151 (MQ=255)
                                 agcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgataatgat                                                                                                                  >  2:2046273/1‑150 (MQ=255)
                                   caAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATgataatgataat                                                                                                               >  1:572545/1‑151 (MQ=255)
                                   aaatcaaaaaaaTTAACAGAAACAGCGTTCCAACAAACAAAACCCGCCGGAAAAAACTCTGGGAAATAAAACGAATACCCCCTGGTTAAACGTTATTTCCCTGGGGAGTAGAACCCTTTAACCCGCCCCAACAGAGaatgataatgat                                                                                                                  <  1:1708001/140‑1 (MQ=255)
                                    aaaCGGAAAAAGTATAAGAAACAGCGTTGCAGAAAACAAAATCTCCCAAAATATCCGCTGGAAAATAAACCGGGTATCCCCATGTAAAACTTTTTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATgataatgataat                                                                                                               <  2:597028/150‑1 (MQ=255)
                                     aaCTGAAAAATTATAAAAAACAGGGTTGCAGAAAACACAACCTGCAAAAAAAGACGCGGGAAAATAAACCGGGTATCACCTTGTAAAAGGTTGTTACGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGAtaataa                                                                                                             <  1:38154/151‑1 (MQ=255)
                                      aCTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGAtaataa                                                                                                             >  2:244488/1‑150 (MQ=255)
                                       cggaaaaaTTATCAGAAACAGCGTTGCAGAAAACACACTCTGCCAGAATATACGTGGGAAATTAAAACGGAGATCACTATGTGAAACTTTTTTCCGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATaaa                                                                                                            <  2:1076280/146‑1 (MQ=255)
                                        tGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAAt                                                                                                          >  1:1932938/1‑151 (MQ=255)
                                           aaaaGTATCAGAAAAAGCGTTGCAGAAAACACAATCGGCCAGAATATACGCGGGAAAATAAACCGGGTATCCCCATGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTa                                                                                                        <  1:388536/150‑1 (MQ=255)
                                             aaGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGcc                                                                                                     >  2:1270567/1‑151 (MQ=255)
                                               gTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                                   >  2:1526469/1‑151 (MQ=255)
                                                tataaaaaaCAGCGTGCGAGAAAACCCAATCCCCCAAAATATGCGGGGGGAAATAAAACGGTTACCCCTATTTAAAAGGTTGTAACGCGGGGGAGTAAAGACCCTGTAGCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCg                                                                                                    <  2:343364/143‑1 (MQ=255)
                                                tataAGAAACAGCGTTGCAGGAAACAAAACCTGCAAAAATATACGCTGGAAAATAAAACGGGTCTCACCATGAGAAACTTTTTAAGGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCg                                                                                                    <  1:1052984/149‑1 (MQ=255)
                                                    aaaaccAGCGTTCCAGAAACCACAATCTCCAAAAAAATCCGCGGGAAAATAAACGGGTTACCACCATGAGAAGTTTTTTACCGCGGGGGAAAAAACCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATg                                                                                              <  1:1766352/146‑1 (MQ=255)
                                                    aaaaaCAGCGTTGAAGAAAACACAATCGCCAAGAAAATCCGCTGGAAAATAAAAGGGTTACCACCATGTAAAACGTTTTTACGCGGGGGAATAAAACCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATg                                                                                              <  1:1790456/149‑1 (MQ=255)
                                                    aGAAACAGCTTTGCAGAAAACACAACCTGCCAAAATATCCGCTGGAAAATAAAACGGGTACCACAATGTAAACCTTTGTTACGCGGGGGAGTAAAGCCCTTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTaa                                                                                                <  2:1881083/149‑1 (MQ=255)
                                                    aGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAAt                                                                                               >  2:403250/1‑150 (MQ=255)
                                                    aGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGAAAATAAACCGGGTATCACCATGTAAAAGGTTGTTACGCGGGGGATAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATg                                                                                              <  1:1683332/151‑1 (MQ=255)
                                                       aaaaccGTTGAAGAAACCACATTCGGCCAAAATATCCGCGGGAAAATTAAGCGTTTATCACCGTCTGAAACTTTGTTACCCGGGGAAAAAAAGCCCCTATAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATa                                                                                           <  1:1639138/146‑1 (MQ=255)
                                                       aaCAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGAt                                                                                            >  1:285629/1‑150 (MQ=255)
                                                          aacGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCg                                                                                        >  2:1727874/3‑151 (MQ=255)
                                                          aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCg                                                                                        >  2:1925047/1‑151 (MQ=255)
                                                          aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCg                                                                                        >  2:194159/1‑151 (MQ=255)
                                                           gCGTTGCAGGAAACACAATCTGCCAGAATATCCGCGGGAAAATAAACCGGGTATCACCATGTGAAACGTTTTTACGCCGGGGAAAAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc                                                                                         <  1:1472185/149‑1 (MQ=255)
                                                             gTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                                     >  1:109667/1‑151 (MQ=255)
                                                              ttGCAGAAAACACAATCCGCAAAAATATCCGCTGGAAAATAAAAGGGTTATCACCATGTAAACCTTTGTAGCGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGt                                                                                       <  2:956682/148‑1 (MQ=255)
                                                               tGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCg                                                                                   >  2:1833984/1‑151 (MQ=255)
                                                                gCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGt                                                                                  >  2:655360/1‑151 (MQ=255)
                                                                 cAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGt                                                                                  >  2:565634/1‑150 (MQ=255)
                                                                     aaaCACAATCTGCAAGAATACACTTGGGAAAATTAAACGGGTACCCCCCTGGAAAGCTTTTTCATGCGGGGGAGTAAAGCCCCTTAATCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                                     <  2:1101882/143‑1 (MQ=255)
                                                                         acaATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc                                                                          >  1:1794244/1‑150 (MQ=255)
                                                                                 ccAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGACACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAAAAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGc                                                                  >  1:920770/1‑150 (MQ=255)
                                                                                 caaaaataaccgcgGGAAAATAAAACGGTTACCCCCATGTAAAACGGTTTTCCGCGGGGGAAAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTg                                                                   <  1:1297213/139‑1 (MQ=255)
                                                                                 caaaaaTATACGCGGGAAAATAAAACGGGTATCCCCATGTAAAAGGTTGTTACGCGGGGAAAAAAAGACCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTg                                                                   <  1:811724/145‑1 (MQ=255)
                                                                                 caaGAATATACGCGGGAAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCTTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGc                                                                  <  2:1137081/148‑1 (MQ=255)
                                                                                  cAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCat                                                                >  1:954635/1‑151 (MQ=255)
                                                                                    gAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                              >  2:428484/1‑151 (MQ=255)
                                                                                     aaTATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                             >  2:1597011/1‑151 (MQ=255)
                                                                                      atatACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAttt                                                            >  1:213889/1‑151 (MQ=255)
                                                                                      atatACGCTGGAAAATAAAACGGGTATCCCCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAttt                                                            <  2:1384572/151‑1 (MQ=255)
                                                                                         tgcgcgGGAAAATAAAACGGGTATCCCCATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCCTTAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                              <  1:1669303/144‑1 (MQ=255)
                                                                                         taCGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCg                                                         >  1:660736/1‑151 (MQ=255)
                                                                                         taCGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCg                                                         >  2:487609/1‑151 (MQ=255)
                                                                                           cgcgGGAAAATAAAACGGGTACCCCCATGTAAACCGTTGTTGCGCGGGGGAGAAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGa                                                       <  2:1826934/151‑1 (MQ=255)
                                                                                            gcgggaaaaTAAAACGGGTATCACCATGGGAAAGTTTGTTACGCGGGGAAGAAAAACCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGc                                                          <  1:485708/141‑1 (MQ=255)
                                                                                             gctggAAATTAAAGCGGGTATCACCATGGAAAACGTTGTCACGCCGGGAAGAAAAGACCCTAAAGCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGa                                                       <  2:1121204/147‑1 (MQ=255)
                                                                                             cTGGGAAATTAACCGGTTTTATCGATGTAAAGCTTTATCCCGCCGGGGAATAAGATCCCGGTAACCGCCCCAACAGAGAGTGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTGTGCATATTTg                                                           <  2:353217/145‑1 (MQ=255)
                                                                                             cTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGAc                                                      >  2:1236242/1‑150 (MQ=255)
                                                                                                ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                                   >  1:2025267/1‑150 (MQ=255)
                                                                                                   aaatAAACGGGTTATCACCATGTGAAACGTGGTTACGCGGGGGAGAAAGGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaaa                                                  <  2:1380798/146‑1 (MQ=255)
                                                                                                   aaTTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATc                                               >  1:869695/1‑151 (MQ=255)
                                                                                                   aaTAAAAGGGGTATCACAATGTAAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATc                                               <  2:1635725/151‑1 (MQ=255)
                                                                                                    aTTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATcc                                              >  2:94657/1‑151 (MQ=255)
                                                                                                        aaaCGGGTATCACCATGAGAAACGTTGTTACGCGGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTa                                           <  1:413516/151‑1 (MQ=255)
                                                                                                        aaaCGGGTACCACCAGGTAAACGGTTGTTACGCGGGGGAGAAAAGACCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTa                                           <  1:1311629/151‑1 (MQ=255)
                                                                                                         aaCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTa                                           >  1:1264116/1‑150 (MQ=255)
                                                                                                          aCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTa                                           >  2:1567647/1‑149 (MQ=255)
                                                                                                               gtatgccgattgtgAAGAGTTGTGACCCCGGGTAGAAAAACCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTaa                                          <  1:181782/137‑1 (MQ=255)
                                                                                                                    cacaaggtaaaaCTTTTTTACCCGGGGGAGAAAAACCCTTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATc                                   <  1:1436654/140‑1 (MQ=255)
                                                                                                                      ccATGTGAAACTTTGTTACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCt                                <  2:1799973/150‑1 (MQ=255)
                                                                                                                          gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCa                            >  2:1782733/1‑150 (MQ=255)
                                                                                                                            aaaaCGTTTTTACGCGGGGGAAAAAACCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATg                         <  1:1225936/150‑1 (MQ=255)
                                                                                                                             aactttttttACCCGGGGTAGTAAACCCCCTTAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGc                             <  1:1260511/138‑1 (MQ=255)
                                                                                                                                  tttttccgcgGGGGAAAAAACCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgcc                    <  2:1946779/144‑1 (MQ=255)
                                                                                                                                    gttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc                   >  2:1111953/1‑149 (MQ=255)
                                                                                                                                     ttACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAg                 <  1:771430/150‑1 (MQ=255)
                                                                                                                                        cGCCGGGGAGTAAAGCCCTTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCt               <  2:1070754/149‑1 (MQ=255)
                                                                                                                                         gCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGa             >  1:1808870/1‑150 (MQ=255)
                                                                                                                                           gggggAGAAAGCCCCTGGATCCGGCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCt               <  2:93509/144‑1 (MQ=255)
                                                                                                                                            ggggAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa         <  1:1683973/151‑1 (MQ=255)
                                                                                                                                             gggAGTAAACCCCTTGTATCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAc        <  2:1617846/151‑1 (MQ=255)
                                                                                                                                             ggaaGTAAAGACCCTTAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAc        <  2:1436703/148‑1 (MQ=255)
                                                                                                                                              ggAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAca       <  2:1931752/151‑1 (MQ=255)
                                                                                                                                                  tAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcat   >  1:43370/1‑151 (MQ=255)
                                                                                                                                                   aaaGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatc  <  1:1000388/151‑1 (MQ=255)
                                                                                                                                                   aaaGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatca    <  1:91904/149‑1 (MQ=255)
                                                                                                                                                    aaGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatc  >  1:1091512/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
ATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATC  >  NC_000913/2628477‑2628771

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