Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8160GA46.9% 78.9 / 114.8 113intergenic (+75/+144)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (54/6);  new base A (0/53);  total (54/59)
Fisher's exact test for biased strand distribution p-value = 1.40e-25
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

TTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTT  >  NC_000913/1622666‑1622961
                                                                                                                                                            |                                                                                                                                                  
ttttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCg                                                                                                                                                    >  2:1712865/1‑151 (MQ=255)
 tttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCg                                                                                                                                                    >  2:137748/1‑150 (MQ=255)
  ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑a                                                                                                                                                  >  1:1130169/1‑150 (MQ=255)
    aTGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑aaa                                                                                                                                                <  1:1758582/150‑1 (MQ=255)
     tGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑aaaa                                                                                                                                               >  1:723411/1‑150 (MQ=255)
     tGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑aaaa                                                                                                                                               >  1:1856673/1‑150 (MQ=255)
      ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑aaaaaa                                                                                                                                             <  2:324290/151‑1 (MQ=255)
      ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑aaaaa                                                                                                                                              >  2:1390200/1‑150 (MQ=255)
         cGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAAc                                                                                                                                          >  2:1225085/1‑151 (MQ=255)
          gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAAcc                                                                                                                                         <  2:332208/151‑1 (MQ=255)
           tACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCg                                                                                                                                        >  1:630737/1‑151 (MQ=255)
            aCTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGc                                                                                                                                       <  1:1327045/151‑1 (MQ=255)
               gggAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAg                                                                                                                                    <  1:1124002/151‑1 (MQ=255)
               gggAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAg                                                                                                                                    <  1:1162289/151‑1 (MQ=255)
                   ataatTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGt                                                                                                                                   >  2:697640/1‑148 (MQ=255)
                     aatTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAAcc                                                                                                                              >  1:1063246/1‑151 (MQ=255)
                      attttattttCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACcg                                                                                                                             >  2:1575822/1‑151 (MQ=255)
                         ttattttCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACcggcg                                                                                                                          >  1:1002399/1‑151 (MQ=255)
                            ttttctggttttctggCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGt                                                                                                                        >  2:1260685/1‑150 (MQ=255)
                            ttttctggttttcgggCCG‑CTTATTTCAGCCAAAAAG‑GGGATGACTAATGAACGGAAATAACCCCCCCCCTAACCGGCCCT‑TTTTTACGGTTGTGTAAAGGGGGCT‑TTTTTTTT‑AGGAG‑GGCA‑AAAAAAAACCGCCAGTAAACcggcgg                                                                                                                         <  2:1928804/149‑1 (MQ=255)
                               tctggttttctggCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGt                                                                                                                        >  2:823260/1‑147 (MQ=255)
                                ctggttttctggCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATg                                                                                                                   >  2:914231/1‑151 (MQ=255)
                                  gtttttcgggCCG‑TGTATTTCAGCAAAAAAG‑GGGATGATAAATGAAGGGAAATACCCCCTCCCCTACCCGCCCCC‑TTGTTACGGTTGTGTACAGGGGGCT‑TGTTTTTT‑ATGCC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATg                                                                                                                   <  1:960025/147‑1 (MQ=255)
                                    ttttctggCCG‑TTTATTCCACCCACAATG‑GGGATGATAAATGACCGAAGATACCCCCTCACCAACCCGCCCCT‑TTTTTACGGTTGTGTAAAGGGGGCT‑TGTTTTTT‑ATCAG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                                <  2:590167/150‑1 (MQ=255)
                                    ttttctggCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCCGTAAACCGGCGGGGGATGCtt                                                                                                                >  1:403619/1‑150 (MQ=255)
                                    ttttctgcCCG‑GTTTTTCCACCCAAAAAG‑GGGATGACAAATGAGGGGAAATAACCCCCCCCCTACCCGCCCCT‑TTTTTACATTGGTGTAAAGGGGGCT‑TTTTTTTT‑GGGAG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                                <  2:1221473/150‑1 (MQ=255)
                                     tttctggCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                              >  2:533473/1‑151 (MQ=255)
                                     tttctggCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                              >  2:1135988/1‑151 (MQ=255)
                                      tttcggCCGCCGTATTACA‑CCCCAAATGGGGATGGGTAATAAGCGGAAATAACCCCCCCCCAACCCGGCCCT‑TTCTTACGGTTTTCTA‑AGGGGGTC‑TTGTTTTT‑TTGGCGGGCA‑AAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTg                                                                                                               <  2:399201/146‑1 (MQ=255)
                                      ttcgggCCG‑CTTTTTTCAGCCAAAAAG‑GGGATGACTAATGAACGGAGATACCCCCTCCCCTACCCGCCCCC‑TTTTTACGGTTGTGTACAGGGGGCT‑TTTTTTTT‑ATGGG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                              <  2:773667/150‑1 (MQ=255)
                                          ggCCG‑CTTTTTCCAGCCAAAAAG‑GGGATGATTAATGAAGGGAAAAAACCCCCCCCCAACCCGCCCCC‑TTTTTACAGTTGTGAACAGGGGGCT‑TGTTTTTT‑ATGAG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa                                                                                                         <  2:1327039/151‑1 (MQ=255)
                                            ccG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                        >  2:1509576/1‑150 (MQ=255)
                                            ccG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                        >  2:161423/1‑150 (MQ=255)
                                             cggggtttttCACCCAAAAGG‑GGGATGATTAATAAACGAAAAAACCCCCTCCCCTACCCGCCCCT‑TTGTAACAGTTGTGTAAAGGGGGCT‑TGTTTTTT‑AGGAG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATg                                                                                                           <  2:1238231/140‑1 (MQ=255)
                                              g‑cgTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                        >  2:160719/1‑148 (MQ=255)
                                                 gTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                                   >  2:1349001/1‑151 (MQ=255)
                                                  tATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt                                                                                                  >  1:966282/1‑151 (MQ=255)
                                                   aTTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt                                                                                                  >  1:1044347/1‑150 (MQ=255)
                                                    tttCACCAAAAAAG‑GGGATGAATAAGGAACGGAAATAATCCCCCCCCTACCCGCCCCT‑TTGTTACATTTGTTTAAAGGGGGCC‑GTATTTTT‑ATGCG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAt                                                                                                    <  2:1633732/147‑1 (MQ=255)
                                                     tccAGCCACAAGG‑GGGATGACAAATGAGGGAAGATACCCCCCCCCCAACCCGCCCCT‑TTTTTACTGTTGTGTAAAGGGGGCT‑TTTTTTTT‑ATGGG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtg                                                                                                 <  1:792325/147‑1 (MQ=255)
                                                     tccAGCCAAAAAG‑GGGATGACAAATGAACGGAAAAACCCCCTCCCCTACCCGCCCCT‑TTTTTACAGTTGTGTAAAGGGGGCC‑TTTTTTTT‑ATGAG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtgt                                                                                                <  2:1097394/148‑1 (MQ=255)
                                                         ccaaaaagt‑GGGATTATTAAAGAAGGAAAAT‑ATCCCCCCCCAAACCCGCCCCTTTGTTACGGTTTTTTACAGGGGGCT‑TGTTTTTA‑AGGCG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtttt                                                                                           <  2:276834/143‑1 (MQ=255)
                                                          ccACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGGTAGATTTGTGTTTTg                                                                                          >  1:1345128/1‑151 (MQ=255)
                                                          caaaaaaG‑GGGATGAATAAAGAAGGGAAATAACCCCTCCCCTACCCGCCCCT‑TTGTTAGATTTGTTTAAAGGGGGCT‑TGTTTTTT‑ATGAG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTg                                                                                          <  2:1505948/147‑1 (MQ=255)
                                                            acaAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGCTAGATTTGTGTTTTGCt                                                                                        >  2:1112808/1‑151 (MQ=255)
                                                            acaAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCt                                                                                        >  1:1624742/1‑151 (MQ=255)
                                                            acaAAG‑GGGATAAATAAGGAAGGGAAATACCCCCCCCCCAACCCGCCCCT‑TTTTTACTTTTGT‑TAAAAGGGGGC‑TTTTTTTTGATGCC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGc                                                                                         <  1:1796857/150‑1 (MQ=255)
                                                             aaaaT‑GGGATGATTAAGGAAGGGAAAAAACCCCCCCCCTACCCGCCCCT‑TTGTTACGGTTGTGTAAAGGGGTCT‑TGTTTTTT‑ATGCG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                       <  1:427259/150‑1 (MQ=255)
                                                                aT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                                     >  2:580911/1‑150 (MQ=255)
                                                                   gggATGAAAAATGAAGGAAAAAAACCCCTCCCCAACCCGGCCTT‑TTTTTACAGTTGTGAAAAGGGGTCT‑TTTTTTTT‑ATGCC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                                     <  1:1867173/148‑1 (MQ=255)
                                                                   ggaatgaaaaaggAACGAAAATATCCCCTCCCCTAACCGCCCCT‑TTTTTACAGTTGTGTACAGGGGGCT‑TTTTTTTT‑ATGAG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTa                                                                                    <  1:57123/137‑1 (MQ=255)
                                                                     gATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt                                                                                >  2:735606/1‑151 (MQ=255)
                                                                        gACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTaa                                                                              >  2:1252450/1‑150 (MQ=255)
                                                                        gACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTaa                                                                              >  2:1957537/1‑150 (MQ=255)
                                                                         aCTAATGAAGGGAAATACCCCCCCCCCAACCCGCCCCT‑TTGTAACTGTTGTATACAGGGGCCT‑TTTTTTTT‑ATGAG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTaa                                                                              <  1:722701/149‑1 (MQ=255)
                                                                          cTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                             >  1:825924/1‑149 (MQ=255)
                                                                           tAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAgg                                                                          >  1:1049825/1‑151 (MQ=255)
                                                                            aaTGAACGGAGATAATCCCTCCCCAACCCGCCCCT‑TTTTTAGATTTGTAAACAGGGGGTG‑TGTTTTTT‑AGGAC‑GACA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAg                                                                           <  1:1941854/149‑1 (MQ=255)
                                                                            aaTGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAgg                                                                          >  1:973175/1‑150 (MQ=255)
                                                                              tGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCGAAAAAAAAACCCCCAGTAAACCGCCGGGGAATGCTTGAATGGTTAGATTTGTGTTTTGTTTTTAcccaaaaag                                                                           >  1:1030718/1‑140 (MQ=255)
                                                                              tGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAt                                                                       >  2:916758/1‑151 (MQ=255)
                                                                                aggggAAATACTCCCCCACCTACCCGCCCCT‑TTGTTACGGTTGTGAAAAGGGGGCT‑TTTTTTTT‑AGGGG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtt                                                                      <  1:953631/147‑1 (MQ=255)
                                                                                  cGGAAATAACCCCACACCTACCCGCCCCT‑TTGTTACTGTTGTGTACAGGGGGCT‑TGTTTTTT‑ATGGC‑GGCA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt                                                                    <  2:171693/150‑1 (MQ=255)
                                                                                   ggaaatacccccTCACCAACCCGGCCCT‑TTTTTAAGTTTGTGTAAAGGGGGCC‑TTTTTTTT‑ATGAG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCa                                                                        <  2:1642858/136‑1 (MQ=255)
                                                                                     agaTAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCt                                                                 >  2:28495/1‑150 (MQ=255)
                                                                                               cccccTACCCGCCCCT‑TTGTTACGGTTGTGAACAGGGGCCC‑TTATTTTT‑ATGCC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAAc                                                      <  2:1584821/148‑1 (MQ=255)
                                                                                                  ccaaccccccccc‑TTGTGACGGTTGAGCAAAAGGGCCT‑TTTTTTTT‑ACGGC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGaa                                                   <  1:1915555/142‑1 (MQ=255)
                                                                                                  ccaacccGCCCCT‑TTTTTACATTTGTGTACAGGGGGCT‑TTTTTTTT‑ATGAG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                                    <  2:1066613/145‑1 (MQ=255)
                                                                                                  ccTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                                    >  2:1106558/1‑150 (MQ=255)
                                                                                                       ccGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                              >  1:1292311/1‑151 (MQ=255)
                                                                                                       ccGCCCCT‑TTTTTACAGTTGTGTAAAGGGGGCT‑TTTTTTTT‑AGGAC‑GGCA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGt                                               <  2:580745/150‑1 (MQ=255)
                                                                                                       ccGCCCCT‑TTTTTACAGTTGTAAAAAGGGGGCT‑TTTTTTTT‑AGGCG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                                    <  1:1892821/145‑1 (MQ=255)
                                                                                                         ggCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACCCTAAAAGGCATTTTCCTGCACTGATAACGAATCGTTg                                             >  1:348719/1‑150 (MQ=255)
                                                                                                           ccct‑TTTTTACATTTGTGTAAAGGGGGCT‑TGTTTTTT‑ATGAG‑GGCA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                             <  2:1269950/148‑1 (MQ=255)
                                                                                                           cccg‑TTGTTTTATTTGTGGAGAGGGGGCC‑TTTTTTTA‑AGGAC‑GCGA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGaca                                          <  2:1141824/151‑1 (MQ=255)
                                                                                                                tttttAGAGTTGTGTAAAGGGGGCT‑TTTTTTTT‑ATGCG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAg                                       <  1:1260662/147‑1 (MQ=255)
                                                                                                                   ttacttttGTGTA‑AAGGGGCCTTTTTTTTT‑ATGGC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTa                                     <  1:1340853/143‑1 (MQ=255)
                                                                                                                   ttACGGTTGTAGAAAGGGGGCT‑TTTTTTTT‑AGGAG‑GGAA‑AAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt                                      <  1:90307/148‑1 (MQ=255)
                                                                                                                     actgTTGTGGACAGGGGGCT‑TTTTTTTT‑ATGAC‑GAAA‑AAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc                                   <  1:709936/146‑1 (MQ=255)
                                                                                                                     aCAGTTGTGTA‑AAGGGGCCTTGTTTTTT‑AGGAC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                                 <  2:534096/151‑1 (MQ=255)
                                                                                                                      cAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                                >  2:102547/1‑151 (MQ=255)
                                                                                                                         ttgtaaaaagggggCC‑TTTTTTTT‑AGGAG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCa                               <  1:671645/139‑1 (MQ=255)
                                                                                                                             gtACAGGGGTCC‑TTTTTTTT‑GGGAC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCa                               <  2:1602868/145‑1 (MQ=255)
                                                                                                                               aaaaGGGGCCTTGTTTTTT‑AGGAC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtt                            <  2:1505460/145‑1 (MQ=255)
                                                                                                                               aCAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc                       >  2:1494153/1‑151 (MQ=255)
                                                                                                                                 agggggCG‑TGTTTTTT‑ATGAC‑TGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaa                     <  1:1047328/149‑1 (MQ=255)
                                                                                                                                  aGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaat                    >  2:632974/1‑151 (MQ=255)
                                                                                                                                   ggtttc‑TGTTTTTT‑AAGAC‑GAAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaa                     <  1:1575791/144‑1 (MQ=255)
                                                                                                                                   ggggCT‑TGTTTTTT‑ATGCG‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct                        <  2:926582/146‑1 (MQ=255)
                                                                                                                                    gggCT‑TGTTTTTT‑AAGCC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatga                  <  2:1252265/151‑1 (MQ=255)
                                                                                                                                        cttttttttt‑ATGAC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat                <  2:1492618/145‑1 (MQ=255)
                                                                                                                                        c‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaatg               >  1:1901477/1‑150 (MQ=255)
                                                                                                                                          tttttttt‑tAGAC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaat                <  2:1778800/145‑1 (MQ=255)
                                                                                                                                          tgtttttt‑AGGAC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaatg               <  1:1111098/146‑1 (MQ=255)
                                                                                                                                             ttttt‑ATGAC‑GGAA‑AAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGACATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTa            <  1:1957318/149‑1 (MQ=255)
                                                                                                                                             ttttt‑ATGAC‑GAAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaatg               <  1:1593580/147‑1 (MQ=255)
                                                                                                                                                tt‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg       >  1:1509973/1‑151 (MQ=255)
                                                                                                                                                 t‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg       >  1:1434992/1‑150 (MQ=255)
                                                                                                                                                     gAC‑GGAA‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCGTTCCTGGATAGATTTGTGTTTTGTTTTTACGCTAACAGGCATTTTCCCGCACTGATAACGAAACGGTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGa      <  1:262112/148‑1 (MQ=255)
                                                                                                                                                       c‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCtt  >  1:428658/1‑150 (MQ=255)
                                                                                                                                                            |                                                                                                                                                  
TTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCG‑CGTATTTCAGCCACAAAT‑GGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCC‑TTGTTACAGTTGTGTACAAGGGGCC‑TGATTTTT‑ATGAC‑GGCG‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTT  >  NC_000913/1622666‑1622961

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: ‑