Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,287,7100GA43.4% 77.6 / 48.9 77intergenic (‑339/‑54)acs/nrfAacetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base G (36/7);  new base A (0/33);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 5.84e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGA  >  NC_000913/4287562‑4287843
                                                                                                                                                    |                                                                                                                                     
gcTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGa                                                                                                                                      >  2:1141254/1‑150 (MQ=255)
     tGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAAc                                                                                                                                 <  2:110580/150‑1 (MQ=255)
        aaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                                                                                                              >  1:952136/1‑150 (MQ=255)
        aaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                                                                                                              >  1:1292142/1‑150 (MQ=255)
           agagGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaaa                                                                                                                           >  1:497615/1‑150 (MQ=255)
             agGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTGGGTGCGGTTACTGTCGTGAAAAACAATAAAAAc                                                                                                                        >  2:1441423/1‑151 (MQ=255)
             agGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAAc                                                                                                                        >  2:963341/1‑151 (MQ=255)
                 aGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGcc                                                                                                                    <  1:1178553/151‑1 (MQ=255)
                   aTACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAAACAACATTCACAGTGTGGTTTTTTGTTACACATGGGGGCGAGCAATGTCATGACAGTGTAGGTGCGTTTACTGTCGTAAAAAACAAAAAAAACCGCCAt                                                                                                                  <  1:415427/151‑1 (MQ=255)
                    tACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCAtt                                                                                                                 >  1:792375/1‑151 (MQ=255)
                      cTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGc                                                                                                               >  1:876483/1‑151 (MQ=255)
                         aCTAACTCAAAAGTGGTATTTTACATGCACTTAAAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTTACACAAAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGca                                                                                                              <  2:1348447/149‑1 (MQ=255)
                            aaCTCTAAAGTGGTTTTTTACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGTTTTTTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGGGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGca                                                                                                              <  2:873814/146‑1 (MQ=255)
                              ctctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATgg                                                                                                       >  1:927381/1‑151 (MQ=255)
                                 tAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGgcgc                                                                                                    >  1:359541/1‑151 (MQ=255)
                                  aaaGTGGTATTTTACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGGGCGGTTACTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCa                                                                                                   <  2:125808/151‑1 (MQ=255)
                                   aaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCa                                                                                                   >  2:839060/1‑150 (MQ=255)
                                     gTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCa                                                                                                   >  2:784719/1‑148 (MQ=255)
                                      tGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTc                                                                                               >  1:1094481/1‑151 (MQ=255)
                                       gttATTTTACATGCACTTACAATTTATAAAAGACAACATTCACGGTTTGGTTTTTTGTTACACATAGGGGCGAGCATTTTCAGGACAGTGTAGGTGGGTTTACTGCGGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCg                                                                                              <  1:5841/149‑1 (MQ=255)
                                         tATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCgg                                                                                             <  1:821349/150‑1 (MQ=255)
                                           ttttACTTGCATTTACAATTAATAAAAGAAACCATTCATGGTTTGGTTTTTTGTAACACAGAGGGGCGAGCAATGTCATAACAGGTTAGTGGCGGTTACGATCATAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTc                                                                                               <  2:694052/146‑1 (MQ=255)
                                           ttttACATGCCTTAAAAATTAATTAAAGAAAACATCCACGGTTTGGTTTTTTTTTACAAATAGGGGCGAGCAATCTAAAAACAGTGTGGGTGCTGTTTCTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGAt                                                                                           <  2:1327288/150‑1 (MQ=255)
                                           ttttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATg                                                                                          >  1:246649/1‑151 (MQ=255)
                                             ttACATGCACTTACAATTGATTAAAGACAACATTCACGGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGa                                                                                         <  1:713383/150‑1 (MQ=255)
                                             ttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGa                                                                                         >  1:1073354/1‑150 (MQ=255)
                                             ttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGa                                                                                         >  2:492626/1‑150 (MQ=255)
                                               acagccatttACATTTGATTAAAGAAATTATGCTAAGTGTGGTTATTTTTTAAACAAAGGGGCTAAAAATGTAATGACAGTGTAGGTGCGGTTGTTGACGTTAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGaa                                                                                        <  1:90074/141‑1 (MQ=255)
                                               aCATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGc                                                                                      >  1:363513/1‑151 (MQ=255)
                                               aCATGCACTTAAAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTTACCCATAGGGGGGAGAAATTTCATGCCAGTGTAGGTGGGTTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGc                                                                                      <  2:1187967/151‑1 (MQ=255)
                                                 aTGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGccc                                                                                    >  1:167112/1‑151 (MQ=255)
                                                    cACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCt                                                                                  >  2:749213/1‑150 (MQ=255)
                                                         aaaaTTGATAAAAGACAACATTCACAGTGTGGTTATTTGTTCCACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAc                                                                             <  1:1156077/148‑1 (MQ=255)
                                                         aCAATTGATTAAAGACAACATTCACAGTGTGGTTTTTTGTTACAAATGGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAc                                                                             <  1:6899/150‑1 (MQ=255)
                                                           aaTTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACa                                                                            >  2:888239/1‑149 (MQ=255)
                                                           aaTTGATTAAAAACAACATTCACAGTGTGTTTTTTTGTTACACATAGGGGCGAGAAATGTCAGGACAGTGTAGGTGCGGTTACTGTCGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACaa                                                                           <  1:1141284/150‑1 (MQ=255)
                                                               gATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGAt                                                                       >  2:1225486/1‑150 (MQ=255)
                                                                 ttAAAGACAACATTCACAGTGTGGTTATTTGTTACACATGGGGGCGAGAAATGTCATGACGGTGTAGGTGCGGTTCCTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATa                                                                      <  1:1409974/149‑1 (MQ=255)
                                                                   aaaGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaaa                                                                   >  2:648290/1‑150 (MQ=255)
                                                                   aaaGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaaa                                                                   >  2:1361033/1‑150 (MQ=255)
                                                                   aaaGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaaaa                                                                  <  1:122943/151‑1 (MQ=255)
                                                                    aaGACAACATTCAGAGTGTGTTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGGGTAGGTGCGGTTTCGGTCGGGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAAt                                                                 <  1:1052132/151‑1 (MQ=255)
                                                                    aaGAAAACATTCACAGTGTGTTTATTTTTAACACATGGGGGCGAGCAATGTCATGACAGTGTAGGTGCGTTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAAt                                                                 <  1:1302503/151‑1 (MQ=255)
                                                                       acaacaTTCACAGTGTGGTTATTTGTTACACATAGGGGGGAGAAATGTCATGACAGTGTAGGTGCGGTTACTTTCGAAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATa                                                                <  2:359535/149‑1 (MQ=255)
                                                                       acaacaTTCACAGGGTGGTTATTTGTTACACATGGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACGGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATaa                                                               <  2:944248/150‑1 (MQ=255)
                                                                            aTTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcac                                                         >  1:453784/1‑151 (MQ=255)
                                                                                  aGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt                                                    >  1:527983/1‑150 (MQ=255)
                                                                                   gtgtGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTAt                                                  >  2:1383188/1‑151 (MQ=255)
                                                                                     gtgGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCt                                                >  1:1210716/1‑151 (MQ=255)
                                                                                       gtttATTTGTACCACAGGGGGGCGAAAAATGACATAACAGTGTGGGTGGGTTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTc                                              <  2:246646/149‑1 (MQ=255)
                                                                                           ttttttcaaaaatgggggCGAGTAATTTAATGACAGTTAGGGTGCGTTTACTGTCGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTAt                                                  <  2:864186/129‑1 (MQ=255)
                                                                                            tttGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                         >  2:493031/1‑151 (MQ=255)
                                                                                            tttGTTACACAGGGGGGCGAGAAATGTCAGGACGGTGGGGGTGGGTTTCCTTTCGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCt                                          <  1:451172/150‑1 (MQ=255)
                                                                                             tggTTACACATAGGGGAGAGAAAGGACAGGCGAGTGAGGGGGCGTTTACTGCCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                         <  2:456507/148‑1 (MQ=255)
                                                                                              tGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAt                                       >  1:1237881/1‑151 (MQ=255)
                                                                                                taacacAGGGGGGCGAGAAATGTCAGGACAGTGTGGGGGGGTTTACTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTa                                        <  2:297438/146‑1 (MQ=255)
                                                                                                taacaaATGGGGGCGAGAATTTCCATGACAGTGTGGGGGCGTTTCCTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTa                                        <  1:1225517/146‑1 (MQ=255)
                                                                                                   cacaTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt                                  >  1:578998/1‑151 (MQ=255)
                                                                                                     aaagggggCAAGAAATGTCAGGACAGTGGAGGTGGGTTTACTGTCGAAAAAAACAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGa                                    <  2:1148729/143‑1 (MQ=255)
                                                                                                      aTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCt                               >  1:276892/1‑151 (MQ=255)
                                                                                                         gggggaaaaaaTTTCATAACGGGGTAGGGGGGTTTACTGTGGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCt                               <  2:69601/139‑1 (MQ=255)
                                                                                                          gggCGAGAAATGTCAGGACGGTGAGGGTGGGTTTCCTGTCGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttttt                           <  1:121990/151‑1 (MQ=255)
                                                                                                            gcgagaaattaaaTGACGGGGGAGGGGGGGTTACTGCCGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttc                         <  1:946403/139‑1 (MQ=255)
                                                                                                              gAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctt                       >  2:591830/1‑151 (MQ=255)
                                                                                                              gAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctt                       >  2:177994/1‑151 (MQ=255)
                                                                                                              gAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                        >  2:1361217/1‑150 (MQ=255)
                                                                                                                  aaTGTCATGACAGTGTGGGGGGGTTTACTGTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                        <  1:963364/146‑1 (MQ=255)
                                                                                                                  aaTGTCAGGACAGTGTGGGTGGGGTTACTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttt                      <  1:556544/148‑1 (MQ=255)
                                                                                                                  aaTGTAATGACGGTGAAGGTGCGGTTACTGTCGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt                    <  2:510990/150‑1 (MQ=255)
                                                                                                                   aTGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttttc                   >  1:1419359/1‑150 (MQ=255)
                                                                                                                   aTGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttttc                   >  2:425856/1‑150 (MQ=255)
                                                                                                                      tCATGACAGTGTGGGTGCGTTTATTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCa                  <  1:981427/148‑1 (MQ=255)
                                                                                                                      ccATGAAATTGTGGGGGGGTTTATTGTCGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                        <  1:1361256/141‑1 (MQ=255)
                                                                                                                        aTGACGGTGTAGGGGGGTTTACTGTGGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt             <  1:1172637/151‑1 (MQ=255)
                                                                                                                                   ggggtgtttACTGCCGAAAAAAAAAAAAAAACCCGCCAATGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACg      <  1:1328433/140‑1 (MQ=255)
                                                                                                                                   ggTGGGGTTACTTTCGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGa  <  2:596684/151‑1 (MQ=255)
                                                                                                                                   ggTGCGGTTCCTGTAGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACg      <  2:192909/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                     
GCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGA  >  NC_000913/4287562‑4287843

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: