Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,056,1150AC42.5% 159.2 / 25.0 94G75G (GGT→GGG) glnLsensory histidine kinase in two‑component regulatory system with GlnG
Reads supporting (aligned to +/- strand):  ref base A (45/9);  new base C (0/40);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 4.25e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGT  >  NC_000913/4055965‑4056263
                                                                                                                                                      |                                                                                                                                                    
ggCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTaaaa                                                                                                                                                      <  1:282224/151‑1 (MQ=255)
  cGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAAcc                                                                                                                                                    >  1:612293/1‑151 (MQ=255)
     gCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTg                                                                                                                                                 >  2:962173/1‑151 (MQ=255)
     gCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTg                                                                                                                                                 >  2:867899/1‑151 (MQ=255)
     gCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTg                                                                                                                                                 <  2:841168/151‑1 (MQ=255)
      cTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGc                                                                                                                                                >  2:832634/1‑151 (MQ=255)
            cTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGcc                                                                                                                                          >  1:570704/1‑151 (MQ=255)
              gTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTc                                                                                                                                        >  2:638732/1‑151 (MQ=255)
                 ccTGATTTAGGCGGCGCTGGTTATCCACCGGAGCCATCCCCAGAAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGAATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAg                                                                                                                                     <  1:957226/151‑1 (MQ=255)
                 ccTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGCCCGGCATACGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCTTTATCGGTAAACCCTTGCCCCGCCTCCAg                                                                                                                                     <  2:160628/151‑1 (MQ=255)
                   tGACTGAAGCGGCGCTGTTTACCCATCGGAGCCACCCCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAAAAAGGATATGCGGGCGCCCGTCGATGCCCAGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCa                                                                                                                                      <  2:796205/148‑1 (MQ=255)
                       ttAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTc                                                                                                                               >  1:342149/1‑151 (MQ=255)
                         aaGCGGCGCTGTTTATCCAGCGGAGCCATCTCCAGCAGGATCATGCCGCCCGGCATACGCTGGGCCGCCACAGAAAGGATATGGGAGCGCCCGTCGATGCCCAGCGTCACTCCTTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCtt                                                                                                                             <  1:622184/151‑1 (MQ=255)
                          aGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTg                                                                                                                            >  2:503398/1‑151 (MQ=255)
                           gcggcgCTGGTTACCCACCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCTTTATTGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTg                                                                                                                            <  2:179501/150‑1 (MQ=255)
                            cggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGc                                                                                                                           >  1:618504/1‑150 (MQ=255)
                             ggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCAt                                                                                                                         >  2:333596/1‑151 (MQ=255)
                             ggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGCCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGCCCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCAt                                                                                                                         <  1:1180996/151‑1 (MQ=255)
                                gcTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATc                                                                                                                        >  2:784043/1‑149 (MQ=255)
                                gcTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATc                                                                                                                        >  2:1268949/1‑149 (MQ=255)
                                   ggTTATCCATCGGAGCCTTCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCt                                                                                                                    <  1:873240/150‑1 (MQ=255)
                                   ggTTATCCATCGGACCCATCTCCAGCAGGATCATGCCGCCCGGCATACGCTGGGCCGTCACAAAAAGGATATGCGGGCCCCCGTCGATGCCCAGCGTCTCTTCTTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTc                                                                                                                   <  2:865146/151‑1 (MQ=255)
                                    gTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTc                                                                                                                   >  2:667921/1‑150 (MQ=255)
                                     ttACCCATCGGAGCCATCCCCAGCAGGATCATGCCGTCCGGCATCGGCTGGGCCGTCACAAAAAGGATATGCCAGCCCCCGCCGATGACCAGCGCCCTTTCGTTATCGGAAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCt                                                                                                                    <  2:567259/148‑1 (MQ=255)
                                      tATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGat                                                                                                                >  2:447983/1‑151 (MQ=255)
                                      tATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGat                                                                                                                >  2:1273244/1‑151 (MQ=255)
                                        tCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGatat                                                                                                              >  2:297542/1‑151 (MQ=255)
                                         cgaccgggcccATCCCCAGCAGGATCATCCCGCCCGGCATACACTGGGCCGTCAAAGAAAGGAAATGCGAGCCCCCGGTGATGCCCAGCCTCCCTCCGTTACCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATAtt                                                                                                             <  1:370110/142‑1 (MQ=255)
                                          cATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCTGCCTCCAGACTTTCTTGCATCAGCTCGATAttt                                                                                                            >  1:1099431/1‑151 (MQ=255)
                                          cATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATAttt                                                                                                            >  2:286404/1‑151 (MQ=255)
                                          cATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATAttt                                                                                                            >  1:706595/1‑151 (MQ=255)
                                            tCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTa                                                                                                           >  2:560348/1‑150 (MQ=255)
                                             cGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAAt                                                                                                         >  1:829282/1‑151 (MQ=255)
                                              ggAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATg                                                                                                        >  1:338805/1‑151 (MQ=255)
                                                aGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATg                                                                                                        >  2:1133951/1‑149 (MQ=255)
                                                aGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATg                                                                                                        >  1:1046177/1‑149 (MQ=255)
                                                  ccATCCCCAGCAGGATCATGCCGCCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATgaga                                                                                                     <  1:376938/150‑1 (MQ=255)
                                                          agcagGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTc                                                                                             >  1:420880/1‑150 (MQ=255)
                                                             agGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAAc                                                                                          >  2:138089/1‑150 (MQ=255)
                                                             agGATCATGCCGTCCGGCATACGCGGGGCCGTCAAAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCTTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCa                                                                                            <  1:525376/148‑1 (MQ=255)
                                                              ggATCATGCCGTCCGGCATACGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACa                                                                                         <  2:376870/150‑1 (MQ=255)
                                                              ggATAATGCCGCCCGGCATACGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGCCCGGCGTCATTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAg                                                                                                <  1:506844/143‑1 (MQ=255)
                                                                aTCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGtt                                                                                      >  1:871062/1‑151 (MQ=255)
                                                                 tCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTc                                                                                     >  2:269836/1‑151 (MQ=255)
                                                                  cATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGTGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTcc                                                                                    >  2:1189362/1‑151 (MQ=255)
                                                                  cATGCCGTCCGGCATACGCGGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGTTAAACCCTTCCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTc                                                                                     <  1:254963/150‑1 (MQ=255)
                                                                   aggccgccCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCCCTTCGTTATCGGTAAAACCTGGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTcc                                                                                    <  2:661450/148‑1 (MQ=255)
                                                                      ccGTCCGGCATACGCGGGGCCGTCACAGAAAGGATAGGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCgg                                                                                  <  1:546647/149‑1 (MQ=255)
                                                                      ccGCCCGGCATACGCGGGGCCGTCACAGAAAGGTTATGCGAGCGCCCGTCGTTGACCAGCGTCACTTCGTTATCGTTAAACCCTTCCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGt                                                                                 <  1:83198/150‑1 (MQ=255)
                                                                          ccGGCATACGCGGGGCCGTCACAGAAGGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCTTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAAc                                                                              <  1:1149940/149‑1 (MQ=255)
                                                                                tccgcgGGGCCGTCAAAAAAAGGATATGCGAGCGCCCGTCGATGCCCAGCGTCATTTCGTTATCGGTAAACCCTTCCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTAc                                                                       <  1:43531/148‑1 (MQ=255)
                                                                                tACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTAcc                                                                      >  1:193482/1‑151 (MQ=255)
                                                                                tACGCTGGGCCGTCAAAGAAAGGATATGCGAGCGCCCGGCGATGCCCAGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTa                                                                        <  1:1270584/149‑1 (MQ=255)
                                                                                    cggggCCGCAACAAAAAGGATATGCGACCGCCCTTCGATGCCCAGCGTCCCTTCGTTATCGGTAAACCCTCCCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGgtgt                                                                         <  1:1156756/142‑1 (MQ=255)
                                                                                       ggCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAAt                                                               <  1:333599/151‑1 (MQ=255)
                                                                                             cacaGAAAGGATATGCGCGCGCCGGTCGATGACCGGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGc                                                           <  1:628904/149‑1 (MQ=255)
                                                                                             cacaGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCgg                                                         >  2:929703/1‑151 (MQ=255)
                                                                                               caGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCggg                                                        >  1:649561/1‑150 (MQ=255)
                                                                                                aGAAGGGATATGCGGGCTCCCGTCGATGACCAGCATCCGTTCTTTATGGATAAACCCTGGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGa                                                       <  2:726325/150‑1 (MQ=255)
                                                                                                aGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTg                                                            <  1:1035965/145‑1 (MQ=255)
                                                                                                   aaGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCtt                                                   >  2:70604/1‑151 (MQ=255)
                                                                                                     ggATATGCGAGCGCCCGTCGATGACCAGCTTCCCTTCTTTATCGTTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCttt                                                  <  1:1036469/150‑1 (MQ=255)
                                                                                                      gATATGCGAGCGCCCGCCGATCGCCAGCGTCTCTTCTTTATCGGTAAACCTTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTg                                                 <  2:1191050/150‑1 (MQ=255)
                                                                                                           gCGAGCGCCCTTCGATGCCCAGCGTCACTTCGTTATGGGAAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCg                                             <  2:1107910/149‑1 (MQ=255)
                                                                                                             gAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCa                                         <  1:770556/151‑1 (MQ=255)
                                                                                                                cgcCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCa                                         <  2:391695/148‑1 (MQ=255)
                                                                                                                 gcCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttg                                     >  1:842731/1‑151 (MQ=255)
                                                                                                                 gcCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttg                                     >  2:329762/1‑151 (MQ=255)
                                                                                                                  cccGTGGATGCCCAGCGTCACTTGTTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgtt                                      <  1:541830/149‑1 (MQ=255)
                                                                                                                  cccGTCGATGACCAGCGTCACTTCGTTATCGGAAAACCCTGGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGc                                          <  2:1096115/145‑1 (MQ=255)
                                                                                                                   ccGTCGATGACCAGCGTCACTTCTTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAgttgtt                                   <  1:1260957/151‑1 (MQ=255)
                                                                                                                    cGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGttgttg                                  >  2:1005446/1‑151 (MQ=255)
                                                                                                                     gTCGATGACCAGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTgc                                 <  2:555205/151‑1 (MQ=255)
                                                                                                                      tCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTgcg                                >  2:974309/1‑151 (MQ=255)
                                                                                                                        gATGACCAGCGTCACTTCGTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTgcgc                               <  1:1199320/150‑1 (MQ=255)
                                                                                                                        gATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGcgcg                              >  2:1147210/1‑151 (MQ=255)
                                                                                                                             ccAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAg                          <  1:798716/150‑1 (MQ=255)
                                                                                                                             ccAGCGTCACTTCGTTATCGATAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAgg                         <  2:149403/151‑1 (MQ=255)
                                                                                                                              cggCGTCACTTCTTTATCGGTAAACCCTTCCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCgg                             <  2:58042/144‑1 (MQ=255)
                                                                                                                              cAGCGTCACTTCTTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAggg                        <  2:1160170/151‑1 (MQ=255)
                                                                                                                                gcgcCACTTCGTTACCGGAAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAgg                         <  2:1042679/148‑1 (MQ=255)
                                                                                                                                 cGTCACTTCTTTATCGGAAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCa                           <  2:10411/145‑1 (MQ=255)
                                                                                                                                 cGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTg                     >  1:829608/1‑151 (MQ=255)
                                                                                                                                  gTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGGTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGtt                      >  2:634108/1‑149 (MQ=255)
                                                                                                                                   tCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGc                   >  2:1156247/1‑151 (MQ=255)
                                                                                                                                    cccTTCGTTATCGGAAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCg                  <  2:387692/149‑1 (MQ=255)
                                                                                                                                    cACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCg                  >  2:826007/1‑151 (MQ=255)
                                                                                                                                    cACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCg                  >  1:269402/1‑151 (MQ=255)
                                                                                                                                       ttCTTTATCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCg                  <  2:666538/148‑1 (MQ=255)
                                                                                                                                       ttCGTTATCGGAAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGt                 <  1:1087667/149‑1 (MQ=255)
                                                                                                                                          gTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATgg            >  1:788791/1‑151 (MQ=255)
                                                                                                                                             aTCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGAt          <  2:1053342/150‑1 (MQ=255)
                                                                                                                                             aTCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGAt          <  1:218388/150‑1 (MQ=255)
                                                                                                                                              tCGGTAAACCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGAt          <  2:658390/149‑1 (MQ=255)
                                                                                                                                                ggTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGTGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGcc      >  1:1154517/1‑151 (MQ=255)
                                                                                                                                                   aaacCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCg        <  1:1249999/146‑1 (MQ=255)
                                                                                                                                                   aaacCCTGGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGa           <  2:1261061/143‑1 (MQ=255)
                                                                                                                                                    aaaCCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGt  >  1:1050615/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAACGGTGTACCAAACAATTTGCGGGAGCTTTGGGCGAGCAGTTGTTGCGCGGCAGGGTTGGCGTAATGGATCGCCAGGT  >  NC_000913/4055965‑4056263

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