Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8160GA44.7% 70.6 / 91.7 94intergenic (+75/+144)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (49/3);  new base A (0/42);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 1.48e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCT  >  NC_000913/1622667‑1622960
                                                                                                                                                            |                                                                                                                                                
tttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCg                                                                                                                                                  <  2:364742/150‑1 (MQ=255)
  tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGaa                                                                                                                                                >  2:1412337/1‑150 (MQ=255)
   aTGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGaaa                                                                                                                                               >  2:648206/1‑150 (MQ=255)
    tGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGaaaa                                                                                                                                              >  1:993199/1‑150 (MQ=255)
     ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGaaaaaa                                                                                                                                            <  1:230142/151‑1 (MQ=255)
     ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGaaaaa                                                                                                                                             >  1:1303681/1‑150 (MQ=255)
      gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGaaaaaa                                                                                                                                            >  1:187097/1‑150 (MQ=255)
      gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGACCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGaaaaaaa                                                                                                                                           <  2:478262/151‑1 (MQ=255)
           aCTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGc                                                                                                                                      >  1:1151051/1‑151 (MQ=255)
            cTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGcc                                                                                                                                     >  2:55762/1‑151 (MQ=255)
                gAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGt                                                                                                                                  >  1:299448/1‑150 (MQ=255)
                    aatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAAcc                                                                                                                             >  2:775611/1‑151 (MQ=255)
                     attttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcg                                                                                                                            >  1:963079/1‑151 (MQ=255)
                            tttttgtttttcgggCCGCTTTTTTCAGCCAAAAAGGGAATAACTAATGAAGGAAAATAACCCCTCCCCAACCCGCCCCT‑TTGTT‑ACG‑G‑TTGTTTACAG‑GGGGCT‑TT‑ATTTTTAGGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                      <  1:565069/150‑1 (MQ=255)
                            tttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                      >  2:330364/1‑150 (MQ=255)
                               ctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATg                                                                                                                  >  2:1214418/1‑151 (MQ=255)
                                   ttttctggCCGCTTATTCCACCCAAAAAGGGGATGACAAATGAAGGGAGATAACCCCTCCCCAACCCGCCCCT‑TTTTT‑ACA‑G‑TTGTGAAAAG‑GGGGCT‑TT‑TTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                               <  2:1311888/150‑1 (MQ=255)
                                    tttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                             >  1:62549/1‑151 (MQ=255)
                                    tttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                             >  1:833746/1‑151 (MQ=255)
                                     ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                             >  1:284625/1‑150 (MQ=255)
                                       ctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                          >  1:345912/1‑151 (MQ=255)
                                       ctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                          >  1:374102/1‑151 (MQ=255)
                                       cgggCCGGGTTTTCCACCCAAAAAGGGGATGATAAATGAAGGGAGATAACCCCTCCCCTACCCGCCCCT‑TTTTT‑AAA‑G‑TTGTGTACAG‑GGGCCT‑TG‑TTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                                           <  2:244384/148‑1 (MQ=255)
                                        tggcgtctttttccAGCCAAAAAGGGGATGACTAATGAAGGGAAATAACCCCCCCCCAACCCGCCCCT‑TTGTT‑ACT‑T‑TTGTGTACAG‑GGGGCT‑TG‑TTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                          <  1:1406942/137‑1 (MQ=255)
                                        gggCCGCTTTTTCCACCCAAAAATGGGATGACTAATGAAGGGAAAT‑ATCCCCCCCCCAACCCGCCCCTTTGTT‑CCA‑T‑TTGTGTAAAG‑GGGGCC‑TT‑TTTTTTATGCGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                         <  2:302587/150‑1 (MQ=255)
                                            cgcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                       >  1:633615/1‑149 (MQ=255)
                                               gTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                                  >  2:1347154/1‑151 (MQ=255)
                                                 aTTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtg                                                                                                >  2:673280/1‑151 (MQ=255)
                                                    tCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTgtgt                                                                                             >  2:940583/1‑151 (MQ=255)
                                                    tCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGGGAAAGCTTGCATGGTTAGGTTTgtgt                                                                                             >  2:453797/1‑151 (MQ=255)
                                                           aaaaGGGGATGACTAATGAGCGGAAATAACCCCTCACCAACCCGGCCCT‑TTGTT‑ACG‑T‑TTGTGAAAAG‑GGGGCT‑TT‑TTTTTTATGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                      <  1:373461/150‑1 (MQ=255)
                                                           aaaaGGGGATGACAAATGAACGGAGATAACCCCTCCCCAACCCGCCCCT‑TTGTT‑ACT‑G‑TTGTGTACAG‑GGGGCT‑TG‑TTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                      <  1:391321/150‑1 (MQ=255)
                                                            aatggggATAAATAATGAAAGGAAAAAACCCCCCCCCAACCCGCCCCC‑TTGTT‑ACA‑‑‑TTGGGGGAAA‑GGGGCCATT‑TTTTTTAGGCCGGAAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt                                                                                     <  1:677939/147‑1 (MQ=255)
                                                              aTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                                    >  1:396107/1‑150 (MQ=255)
                                                                ggtataattAATGAAGGGAAATACCCCCTCCCCAACCCGCCCCT‑TTTTT‑ACG‑T‑TTGTGTACAG‑GGGGCT‑TT‑TTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg                                                                                 <  1:1199014/143‑1 (MQ=255)
                                                                gggATAACTAATGAAGGAAGATACCCCCTCCCCTACCCGCCCCC‑TTGTC‑ACG‑T‑TTGTGTAAAG‑GGGGCG‑TG‑TTTTTTATGGGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg                                                                                 <  1:962709/151‑1 (MQ=255)
                                                                  gATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt                                                                               >  2:47606/1‑151 (MQ=255)
                                                                  gATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt                                                                               >  2:754240/1‑151 (MQ=255)
                                                                     gaatAATGAAGGGAGATAACCCCTCCCCTACCCGCCCCT‑TTGTT‑CCA‑G‑TTGTGTAAAG‑GGGGCT‑TG‑TTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                            <  1:112512/148‑1 (MQ=255)
                                                                     aattAATGAAGGGAGATATCCCCCCCCCTACCCGCCCCT‑TTTTT‑ACA‑T‑TTGTGTAAAG‑GGGGCT‑TT‑TTTTTTATGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                            <  1:5499/148‑1 (MQ=255)
                                                                     aacaaaTGAGAGGAGAAATCCCCCCCCCTACCCGGCCCT‑TTTTT‑ATA‑G‑TTGTG‑ACAAGGGGGCT‑TG‑TTTTTAATGAGGGCAAAAAAACTCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                            <  2:1060348/147‑1 (MQ=255)
                                                                     aaatAATGAAGGGAAATACCCCCCCACCAACCCGGCCCT‑TTTTT‑ACA‑G‑TGGTGAACAG‑GGGGCT‑TG‑ATTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                            <  1:1012317/148‑1 (MQ=255)
                                                                        tAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAgg                                                                         >  1:1332326/1‑151 (MQ=255)
                                                                        tAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAgg                                                                         >  1:734884/1‑151 (MQ=255)
                                                                         aaTGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGc                                                                        >  2:196021/1‑151 (MQ=255)
                                                                           gaaccggagaatttCCCCCCCCTACCCGGCACT‑TTGTT‑AGG‑G‑TTGTGTACAG‑GGGGCT‑TG‑GTTTTTATGCGGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAt                                                                      <  1:1080767/138‑1 (MQ=255)
                                                                             aaCGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAttt                                                                    >  2:332371/1‑151 (MQ=255)
                                                                             aaCGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAttt                                                                    >  1:599462/1‑151 (MQ=255)
                                                                                 gagaTAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCt                                                                >  2:690674/1‑151 (MQ=255)
                                                                                    aataccccctccccaaCCCGCCCC‑TTTTTAACA‑T‑TGGTGAAAAG‑GGGGCT‑TT‑TTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                              <  2:1407400/137‑1 (MQ=255)
                                                                                     tAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                            >  2:1376664/1‑151 (MQ=255)
                                                                                     tAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                            >  1:605201/1‑151 (MQ=255)
                                                                                     tAACCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                            >  1:137534/1‑151 (MQ=255)
                                                                                      aaTCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                           >  2:1438640/1‑151 (MQ=255)
                                                                                      aaTCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                           >  2:227279/1‑151 (MQ=255)
                                                                                       cccccTCCCCTACCCGCCCCC‑TTTTT‑ACA‑G‑TTGTGAACAG‑GGGGCT‑TG‑TTTTTTATGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg                                                          <  1:609460/149‑1 (MQ=255)
                                                                                       cccccTCCCCAACCCGCCCCT‑TTTTT‑ACT‑T‑TTGTGAAAAG‑GGGGCT‑TT‑TTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg                                                          <  1:648212/149‑1 (MQ=255)
                                                                                        ccccccccccaaCCCGCCCC‑TTTTT‑TCC‑GTTTGTGAAAAG‑GGGGCT‑TT‑TTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg                                                          <  1:605742/141‑1 (MQ=255)
                                                                                         ccccccccTACCCGGCCCT‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGCCT‑TG‑ATTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                        <  2:1251357/145‑1 (MQ=255)
                                                                                         ccccccccTACCCGCCCCT‑TTGTT‑ACG‑G‑TTGTGAACAG‑GGGCCT‑TG‑TTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                        <  1:280527/145‑1 (MQ=255)
                                                                                          cccccccaacccGCCCCT‑TTTTT‑ACG‑G‑TTGTGAACAG‑GGGGCT‑TG‑TTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATa                                                       <  1:1430216/141‑1 (MQ=255)
                                                                                            cccccTAACCGCCCCC‑TTGTT‑AGT‑G‑TTTTATCCAG‑GGGGCC‑TG‑TTTTTTAGAACGCAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                     <  2:1288780/148‑1 (MQ=255)
                                                                                             ccacccaCCCGCCCC‑TTTTT‑AAATG‑TTGTGAAAAA‑GGGGCT‑TG‑TTTTTTAGGGGGGAAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                     <  1:290352/145‑1 (MQ=255)
                                                                                              ccctccccGCCCCT‑TTTTT‑ACA‑T‑TTGTATAAAA‑GGGGCT‑TT‑TTTTTTAGGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                                   <  1:775618/145‑1 (MQ=255)
                                                                                               ccTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                                   >  2:18714/1‑150 (MQ=255)
                                                                                                cTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAAt                                                 >  1:365441/1‑151 (MQ=255)
                                                                                                 tAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc                                                >  2:223668/1‑151 (MQ=255)
                                                                                                 tAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc                                                >  1:704796/1‑151 (MQ=255)
                                                                                                 aacccGCCCCT‑TTGTT‑ACG‑G‑TTGTGTAAAG‑GGGGCT‑TT‑TTTTTTATGAGGGAAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc                                                <  1:49365/148‑1 (MQ=255)
                                                                                                 aacccGCCCCT‑TTGTT‑ACA‑G‑TGTTGAACAG‑GGGGCT‑TG‑TTTTTTAGGCCGGGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc                                                <  1:1418341/148‑1 (MQ=255)
                                                                                                  acccGCCCCT‑TTTTT‑ACA‑T‑TTGTGTAAAA‑GGGGCC‑TT‑TTTTTTATGCCGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc                                                <  1:301584/148‑1 (MQ=255)
                                                                                                      ggCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                            >  2:1029995/1‑150 (MQ=255)
                                                                                                              tttt‑ACA‑G‑TTGTGTAAAG‑GGGGCC‑TT‑TTTTTTAGGGGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt                                     <  2:585902/148‑1 (MQ=255)
                                                                                                                 t‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc                                  >  2:744408/1‑150 (MQ=255)
                                                                                                                    cA‑G‑TTGTGTACAG‑GGGGCC‑TG‑TTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                               <  2:579043/151‑1 (MQ=255)
                                                                                                                       g‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                               >  1:858917/1‑149 (MQ=255)
                                                                                                                            tgtACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt                         >  1:1080949/1‑151 (MQ=255)
                                                                                                                            gtgaaaa‑GGGGCCTTG‑TTTTTTAAGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtt                           <  2:1343015/145‑1 (MQ=255)
                                                                                                                             gtACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtc                        >  2:743196/1‑151 (MQ=255)
                                                                                                                               aCAA‑GGGGCT‑TG‑TTTTTATTGAGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                      <  2:187090/151‑1 (MQ=255)
                                                                                                                                 cagggggCC‑TT‑TTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                     <  2:920566/148‑1 (MQ=255)
                                                                                                                                 aa‑GGGGCC‑TGTATTTTGAGGGGGGGAAAAAAAGACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                     <  2:262745/151‑1 (MQ=255)
                                                                                                                                    ggggCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatg                  >  2:657564/1‑151 (MQ=255)
                                                                                                                                            ttttttttATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaatg              <  1:860236/146‑1 (MQ=255)
                                                                                                                                            g‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt            >  1:980382/1‑150 (MQ=255)
                                                                                                                                             tttttttAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt            <  1:1214425/148‑1 (MQ=255)
                                                                                                                                             gttttttGAGACGACAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt            <  1:580527/148‑1 (MQ=255)
                                                                                                                                             gttttttATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt            <  2:702411/148‑1 (MQ=255)
                                                                                                                                              aTTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTa           >  1:1083197/1‑150 (MQ=255)
                                                                                                                                               tttttATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaa          <  1:459385/150‑1 (MQ=255)
                                                                                                                                               tttttAGGACGGAAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaa          <  1:828289/150‑1 (MQ=255)
                                                                                                                                               tttttAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaa          <  1:695352/150‑1 (MQ=255)
                                                                                                                                                ttttATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAc        >  2:107217/1‑151 (MQ=255)
                                                                                                                                                 tttATGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAAc        <  2:291409/150‑1 (MQ=255)
                                                                                                                                                      gACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCt  <  1:1347152/151‑1 (MQ=255)
                                                                                                                                                            |                                                                                                                                                
TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTT‑ACA‑G‑TTGTGTACAA‑GGGGCC‑TG‑ATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCT  >  NC_000913/1622667‑1622960

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