Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,287,7190TA42.4% 66.9 / 32.8 59intergenic (‑348/‑45)acs/nrfAacetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base T (25/9);  new base A (0/25);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 2.19e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTG  >  NC_000913/4287570‑4287869
                                                                                                                                                       |                                                                                                                                                      
aaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACGTAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAAt                                                                                                                                                        >  1:273938/1‑150 (MQ=255)
aaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATa                                                                                                                                                       <  1:689422/151‑1 (MQ=255)
 aaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATaa                                                                                                                                                      <  1:917386/151‑1 (MQ=255)
   agagGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATaaa                                                                                                                                                     >  1:368297/1‑150 (MQ=255)
               tGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGc                                                                                                                                         >  2:821225/1‑150 (MQ=255)
                gACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGca                                                                                                                                        <  2:181098/150‑1 (MQ=255)
                   tAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGcaaca                                                                                                                                     >  1:485647/1‑150 (MQ=255)
                    aaCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGcaa                                                                                                                                       <  2:823776/147‑1 (MQ=255)
                          aaaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAAAAAAAACCGCCATTGCAACAATGgc                                                                                                                                <  2:853476/148‑1 (MQ=255)
                           aggTGGTATTTAACATGCATTTACAATTGATTAAAGAAAACATCCAAAGTGGGGTTATTTGTTACAAATGGGGGCGAGCAATGTATTGCC‑AGTGGAGGTGCGGTTATTGT‑CGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCaa                                                                                                                            <  2:212864/149‑1 (MQ=255)
                            aGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAAt                                                                                                                           >  2:966200/1‑151 (MQ=255)
                                 tATTTTCCATGCCCTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACAAAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTATGGT‑CGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCg                                                                                                                        <  2:47273/149‑1 (MQ=255)
                                   ttttACATTCACTTACAATTGATTAAAGACAACATTCACGGTTTGGTTTTTTGTTCCACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGAt                                                                                                                     <  1:715810/150‑1 (MQ=255)
                                   ttttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATg                                                                                                                    >  1:398358/1‑151 (MQ=255)
                                     ttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGa                                                                                                                   >  2:138130/1‑150 (MQ=255)
                                       aCATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGc                                                                                                                >  1:364223/1‑151 (MQ=255)
                                            cACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCt                                                                                                            >  2:344244/1‑150 (MQ=255)
                                            cACTTACAATTGATTAAAGACAACATTCACAGGGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTa                                                                                                           <  1:237107/151‑1 (MQ=255)
                                             aCTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt                                                                                                          >  2:315734/1‑151 (MQ=255)
                                             aCTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGGGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt                                                                                                          <  2:1105778/151‑1 (MQ=255)
                                               ttaaatttGATTAAAGACAACATTCAAAGGGGGGTTATTTGTTACACATAGGGGCGAGCAATGTCAGGAC‑AGTGTAGGTGGGGTTACGGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt                                                                                                          <  1:1086225/143‑1 (MQ=255)
                                               ttACAATTGATTAAAGACAACATTCCCAGTGTGGTTATTTGTTACACATAGGGGCGGGAAATGTCATGAA‑AGTGTTGGTGGGTTAACTGT‑GGTAAAAAAAAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt                                                                                                          <  2:1024040/149‑1 (MQ=255)
                                               ttACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt                                                                                                          <  1:343187/149‑1 (MQ=255)
                                                     ttGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAgg                                                                                                   >  2:191969/1‑150 (MQ=255)
                                                       gATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGAt                                                                                                 <  1:672773/150‑1 (MQ=255)
                                                         ttAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaa                                                                                               >  2:841059/1‑150 (MQ=255)
                                                          tAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                              >  2:630572/1‑150 (MQ=255)
                                                          tAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                              >  1:879890/1‑150 (MQ=255)
                                                              aaaaacaTTCACAGTGTGGTTATTTGTTACACATCCGGGCGAGCAATGTCATGTC‑AGTGTAGGTCCGGTTACTGT‑CGTGAAAAACACTAAAAACCGCCATTGAAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaaaa                                                                                            >  1:377043/4‑148 (MQ=255)
                                                                   caTTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgc                                                                                     >  2:317403/1‑150 (MQ=255)
                                                                           gtgtgtTTATTTTTTACCCATGGGGGCGAGAAATTTCATGCT‑ATTGTAGGTCTGGTTCCTGT‑CGTAAAAAACATTAAAAACCGCCATTGCAACAATGGCCCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt                                                                              <  2:911605/149‑1 (MQ=255)
                                                                              tgtttATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTGGGTGCGTTTACTGT‑CGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCt                                                                          <  2:896401/147‑1 (MQ=255)
                                                                                 tttttttttACAAATGGGGGGGGGAAATGTCATGAA‑AGTTTAGGTGCGTTTCCTGC‑CGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCa                                                                       <  2:405780/147‑1 (MQ=255)
                                                                                 ttATTTGTTACACATGGGGGCGAGAAATGTCATGAC‑AGTGTGGGGGGGTTTACTGT‑CGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCa                                                                       <  2:342169/150‑1 (MQ=255)
                                                                                     ttttaaaaaagaGGGGCAAGAAATGTAAGAAA‑TGTGTAGGGGCGTTTACTGT‑CGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                                                   <  1:1041569/139‑1 (MQ=255)
                                                                                     tgttgacactgaGGGGCGAGCAAAGTCA‑GACAAGTGGAGGGGGGGTTACTGC‑GTTGAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCt                                                                    <  1:600743/138‑1 (MQ=255)
                                                                                        ttACACATGGGGGCGAGAAATGTCATGAC‑AGTGTGGGTGCGGTTACTGT‑CGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAtt                                                                <  2:193208/150‑1 (MQ=255)
                                                                                        taacaaATGGGGGCGAGAATTTTAATGAC‑GGTGTGGGGGGGTTTCCTGT‑CGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAtt                                                                <  2:585016/148‑1 (MQ=255)
                                                                                                aGGGCCGAGCACTGTCACGAC‑AGCGTAGCTGCGGTTACTGT‑CGTGAAAACCAATAAAAACCGCCATTGCAACAATGGCGCAATTAGGATGAAGCTCCTATGACAAGGATAAAAATAAACGAACGCCGTATCTTCAGCTTATTGATTCCttt                                                       >  1:1059182/1‑151 (MQ=255)
                                                                                                 ggggCGAGAAATGTCATGAC‑AGTGTAGGTGGGTTTACTGT‑GGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtt                                                        <  2:885004/149‑1 (MQ=255)
                                                                                                 ggggCGACAAATGTAATGAC‑AGTGTGGGGGCGGTTACTGT‑CGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtt                                                        <  1:175392/149‑1 (MQ=255)
                                                                                                   ggCGAGAAATGTTTTGAC‑GGT‑TGGGGGGGGTTTCCGTCCGTAAAAAACAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttttt                                                    <  2:104051/151‑1 (MQ=255)
                                                                                                       agaaaTGTCATGAC‑AGTGTAGGTGGGTTTCCTGT‑CGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctt                                                 <  2:1056854/147‑1 (MQ=255)
                                                                                                          aaTGTAATGAA‑AGTTTGGGTGGGTTTACTGG‑CGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                                                  <  1:233754/146‑1 (MQ=255)
                                                                                                            tGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc                                           >  1:1050122/1‑151 (MQ=255)
                                                                                                                aTGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt                                       >  1:14664/1‑151 (MQ=255)
                                                                                                                 ttaC‑AGTGTAGGTGCGTTTACTTT‑CGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACt                                          <  2:981568/145‑1 (MQ=255)
                                                                                                                 tGAC‑GGTGTGGGGGCGGTTACTGT‑GGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt                                       <  2:1009277/150‑1 (MQ=255)
                                                                                                                 tGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTg                                      >  2:581149/1‑151 (MQ=255)
                                                                                                                   gaaggtgtAGGCGGGGTTCTGGT‑CGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttttc                                             <  1:443269/139‑1 (MQ=255)
                                                                                                                            ggTGCGTTTACTGT‑CGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTg                             <  1:483076/150‑1 (MQ=255)
                                                                                                                             gggCGGTTACTGT‑CGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGaa                           <  1:506545/149‑1 (MQ=255)
                                                                                                                              agCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAAc                          >  1:239015/2‑151 (MQ=255)
                                                                                                                               gCGGTTCCTGT‑CGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGaa                           <  2:38834/149‑1 (MQ=255)
                                                                                                                                 ggTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa                        >  1:331874/1‑150 (MQ=255)
                                                                                                                                 ggTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa                        >  1:296010/1‑150 (MQ=255)
                                                                                                                                       tGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctg                 >  1:513833/1‑151 (MQ=255)
                                                                                                                                            gtaaaaaaCAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACgg                    <  2:491604/141‑1 (MQ=255)
                                                                                                                                              gAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCg           >  2:139045/1‑151 (MQ=255)
                                                                                                                                              gAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCg           >  2:181869/1‑151 (MQ=255)
                                                                                                                                              aaaaaaaaaaaaaACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCg           <  1:138100/150‑1 (MQ=255)
                                                                                                                                                       tAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTg  >  2:480137/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
AAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGAC‑AGTGTAGGTGCGGTTACTGT‑CGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTG  >  NC_000913/4287570‑4287869

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