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

TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCT  >  NC_000913/1622667‑1622960
                                                                                                                                                            |                                                                                                                                                
tttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCg                                                                                                                                                  <  2:364742/150‑1 (MQ=255)
  tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaa                                                                                                                                                >  2:1412337/1‑150 (MQ=255)
   aTGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaa                                                                                                                                               >  2:648206/1‑150 (MQ=255)
    tGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaa                                                                                                                                              >  1:993199/1‑150 (MQ=255)
     ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaaa                                                                                                                                            <  1:230142/151‑1 (MQ=255)
     ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaa                                                                                                                                             >  1:1303681/1‑150 (MQ=255)
      gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaaa                                                                                                                                            >  1:187097/1‑150 (MQ=255)
      gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGACCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaaaa                                                                                                                                           <  2:478262/151‑1 (MQ=255)
           aCTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGc                                                                                                                                      >  1:1151051/1‑151 (MQ=255)
            cTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGcc                                                                                                                                     >  2:55762/1‑151 (MQ=255)
                gAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGt                                                                                                                                  >  1:299448/1‑150 (MQ=255)
                    aatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAAcc                                                                                                                             >  2:775611/1‑151 (MQ=255)
                     attttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcg                                                                                                                            >  1:963079/1‑151 (MQ=255)
                            tttttgtttttcgggCCGCTTTTTTCAGCCAAAAAGGGAATAACTAATGAAGGAAAATAACCCCTCCCCAACCCGCCCCTTTGTTACGGTTGTTTACAGGGGGCTTTATTTTTAGGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                      <  1:565069/150‑1 (MQ=255)
                            tttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                                      >  2:330364/1‑150 (MQ=255)
                               ctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATg                                                                                                                  >  2:1214418/1‑151 (MQ=255)
                                   ttttctggCCGCTTATTCCACCCAAAAAGGGGATGACAAATGAAGGGAGATAACCCCTCCCCAACCCGCCCCTTTTTTACAGTTGTGAAAAGGGGGCTTTTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                               <  2:1311888/150‑1 (MQ=255)
                                    tttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                             >  1:62549/1‑151 (MQ=255)
                                    tttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                             >  1:833746/1‑151 (MQ=255)
                                     ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                             >  1:284625/1‑150 (MQ=255)
                                       ctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                          >  1:345912/1‑151 (MQ=255)
                                       ctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                          >  1:374102/1‑151 (MQ=255)
                                       cgggCCGGGTTTTCCACCCAAAAAGGGGATGATAAATGAAGGGAGATAACCCCTCCCCTACCCGCCCCTTTTTTAAAGTTGTGTACAGGGGCCTTGTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                                           <  2:244384/148‑1 (MQ=255)
                                        tggcgtctttttccAGCCAAAAAGGGGATGACTAATGAAGGGAAATAACCCCCCCCCAACCCGCCCCTTTGTTACTTTTGTGTACAGGGGGCTTGTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                          <  1:1406942/137‑1 (MQ=255)
                                        gggCCGCTTTTTCCACCCAAAAATGGGATGACTAATGAAGGGAAATATCCCCCCCCCAACCCGCCCCTTTGTTCCATTTGTGTAAAGGGGGCCTTTTTTTTATGCGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                                         <  2:302587/150‑1 (MQ=255)
                                            cgcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                       >  1:633615/1‑149 (MQ=255)
                                               gTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                                  >  2:1347154/1‑151 (MQ=255)
                                                 aTTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtg                                                                                                >  2:673280/1‑151 (MQ=255)
                                                    tCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTgtgt                                                                                             >  2:940583/1‑151 (MQ=255)
                                                    tCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGGGAAAGCTTGCATGGTTAGGTTTgtgt                                                                                             >  2:453797/1‑151 (MQ=255)
                                                           aaaaGGGGATGACTAATGAGCGGAAATAACCCCTCACCAACCCGGCCCTTTGTTACGTTTGTGAAAAGGGGGCTTTTTTTTTATGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                      <  1:373461/150‑1 (MQ=255)
                                                           aaaaGGGGATGACAAATGAACGGAGATAACCCCTCCCCAACCCGCCCCTTTGTTACTGTTGTGTACAGGGGGCTTGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                      <  1:391321/150‑1 (MQ=255)
                                                            aatggggATAAATAATGAAAGGAAAAAACCCCCCCCCAACCCGCCCCCTTGTTACATTGGGGGAAAGGGGCCATTTTTTTTAGGCCGGAAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt                                                                                     <  1:677939/147‑1 (MQ=255)
                                                              aTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                                    >  1:396107/1‑150 (MQ=255)
                                                                ggtataattAATGAAGGGAAATACCCCCTCCCCAACCCGCCCCTTTTTTACGTTTGTGTACAGGGGGCTTTTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg                                                                                 <  1:1199014/143‑1 (MQ=255)
                                                                gggATAACTAATGAAGGAAGATACCCCCTCCCCTACCCGCCCCCTTGTCACGTTTGTGTAAAGGGGGCGTGTTTTTTATGGGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg                                                                                 <  1:962709/151‑1 (MQ=255)
                                                                  gATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt                                                                               >  2:47606/1‑151 (MQ=255)
                                                                  gATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt                                                                               >  2:754240/1‑151 (MQ=255)
                                                                     gaatAATGAAGGGAGATAACCCCTCCCCTACCCGCCCCTTTGTTCCAGTTGTGTAAAGGGGGCTTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                            <  1:112512/148‑1 (MQ=255)
                                                                     aattAATGAAGGGAGATATCCCCCCCCCTACCCGCCCCTTTTTTACATTTGTGTAAAGGGGGCTTTTTTTTTATGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                            <  1:5499/148‑1 (MQ=255)
                                                                     aacaaaTGAGAGGAGAAATCCCCCCCCCTACCCGGCCCTTTTTTATAGTTGTGACAAGGGGGCTTGTTTTTAATGAGGGCAAAAAAACTCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                            <  2:1060348/147‑1 (MQ=255)
                                                                     aaatAATGAAGGGAAATACCCCCCCACCAACCCGGCCCTTTTTTACAGTGGTGAACAGGGGGCTTGATTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                            <  1:1012317/148‑1 (MQ=255)
                                                                        tAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAgg                                                                         >  1:1332326/1‑151 (MQ=255)
                                                                        tAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAgg                                                                         >  1:734884/1‑151 (MQ=255)
                                                                         aaTGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGc                                                                        >  2:196021/1‑151 (MQ=255)
                                                                           gaaccggagaatttCCCCCCCCTACCCGGCACTTTGTTAGGGTTGTGTACAGGGGGCTTGGTTTTTATGCGGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAt                                                                      <  1:1080767/138‑1 (MQ=255)
                                                                             aaCGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAttt                                                                    >  2:332371/1‑151 (MQ=255)
                                                                             aaCGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAttt                                                                    >  1:599462/1‑151 (MQ=255)
                                                                                 gagaTAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCt                                                                >  2:690674/1‑151 (MQ=255)
                                                                                    aataccccctccccaaCCCGCCCCTTTTTAACATTGGTGAAAAGGGGGCTTTTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                              <  2:1407400/137‑1 (MQ=255)
                                                                                     tAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                            >  2:1376664/1‑151 (MQ=255)
                                                                                     tAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                            >  1:605201/1‑151 (MQ=255)
                                                                                     tAACCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                            >  1:137534/1‑151 (MQ=255)
                                                                                      aaTCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                           >  2:1438640/1‑151 (MQ=255)
                                                                                      aaTCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                           >  2:227279/1‑151 (MQ=255)
                                                                                       cccccTCCCCTACCCGCCCCCTTTTTACAGTTGTGAACAGGGGGCTTGTTTTTTATGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg                                                          <  1:609460/149‑1 (MQ=255)
                                                                                       cccccTCCCCAACCCGCCCCTTTTTTACTTTTGTGAAAAGGGGGCTTTTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg                                                          <  1:648212/149‑1 (MQ=255)
                                                                                        ccccccccccaaCCCGCCCCTTTTTTCCGTTTGTGAAAAGGGGGCTTTTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg                                                          <  1:605742/141‑1 (MQ=255)
                                                                                         ccccccccTACCCGGCCCTTTGTTACAGTTGTGTACAAGGGCCTTGATTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                        <  2:1251357/145‑1 (MQ=255)
                                                                                         ccccccccTACCCGCCCCTTTGTTACGGTTGTGAACAGGGGCCTTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGAt                                                        <  1:280527/145‑1 (MQ=255)
                                                                                          cccccccaacccGCCCCTTTTTTACGGTTGTGAACAGGGGGCTTGTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATa                                                       <  1:1430216/141‑1 (MQ=255)
                                                                                            cccccTAACCGCCCCCTTGTTAGTGTTTTATCCAGGGGGCCTGTTTTTTAGAACGCAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                     <  2:1288780/148‑1 (MQ=255)
                                                                                             ccacccaCCCGCCCCTTTTTAAATGTTGTGAAAAAGGGGCTTGTTTTTTAGGGGGGAAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                     <  1:290352/145‑1 (MQ=255)
                                                                                              ccctccccGCCCCTTTTTTACATTTGTATAAAAGGGGCTTTTTTTTTAGGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                                   <  1:775618/145‑1 (MQ=255)
                                                                                               ccTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                                   >  2:18714/1‑150 (MQ=255)
                                                                                                cTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAAt                                                 >  1:365441/1‑151 (MQ=255)
                                                                                                 tAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc                                                >  2:223668/1‑151 (MQ=255)
                                                                                                 tAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc                                                >  1:704796/1‑151 (MQ=255)
                                                                                                 aacccGCCCCTTTGTTACGGTTGTGTAAAGGGGGCTTTTTTTTTATGAGGGAAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc                                                <  1:49365/148‑1 (MQ=255)
                                                                                                 aacccGCCCCTTTGTTACAGTGTTGAACAGGGGGCTTGTTTTTTAGGCCGGGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc                                                <  1:1418341/148‑1 (MQ=255)
                                                                                                  acccGCCCCTTTTTTACATTTGTGTAAAAGGGGCCTTTTTTTTATGCCGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATc                                                <  1:301584/148‑1 (MQ=255)
                                                                                                      ggCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                            >  2:1029995/1‑150 (MQ=255)
                                                                                                              ttttACAGTTGTGTAAAGGGGGCCTTTTTTTTAGGGGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt                                     <  2:585902/148‑1 (MQ=255)
                                                                                                                 tACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc                                  >  2:744408/1‑150 (MQ=255)
                                                                                                                    cAGTTGTGTACAGGGGGCCTGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                               <  2:579043/151‑1 (MQ=255)
                                                                                                                       gTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                               >  1:858917/1‑149 (MQ=255)
                                                                                                                            tgtACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt                         >  1:1080949/1‑151 (MQ=255)
                                                                                                                            gtgaaaaGGGGCCTTGTTTTTTAAGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtt                           <  2:1343015/145‑1 (MQ=255)
                                                                                                                             gtACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtc                        >  2:743196/1‑151 (MQ=255)
                                                                                                                               aCAAGGGGCTTGTTTTTATTGAGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                      <  2:187090/151‑1 (MQ=255)
                                                                                                                                 cagggggCCTTTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                     <  2:920566/148‑1 (MQ=255)
                                                                                                                                 aaGGGGCCTGTATTTTGAGGGGGGGAAAAAAAGACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa                     <  2:262745/151‑1 (MQ=255)
                                                                                                                                    ggggCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatg                  >  2:657564/1‑151 (MQ=255)
                                                                                                                                            ttttttttATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaatg              <  1:860236/146‑1 (MQ=255)
                                                                                                                                            gATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGtt            >  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)
                                                                                                                                                            |                                                                                                                                                
TTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCT  >  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: