Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,287,7190TA42.2% 71.9 / 37.1 64intergenic (‑348/‑45)acs/nrfAacetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base T (23/14);  new base A (0/27);  total (23/41)
Fisher's exact test for biased strand distribution p-value = 4.30e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGC  >  NC_000913/4287569‑4287855
                                                                                                                                                       |                                                                                                                                        
aaaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                                                                                                                          <  2:1093185/151‑1 (MQ=255)
aaaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                                                                                                                          <  1:87454/151‑1 (MQ=255)
aaaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                                                                                                                          >  1:471565/1‑151 (MQ=255)
 aaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATa                                                                                                                                         <  1:475353/151‑1 (MQ=255)
     gagGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaaaaa                                                                                                                                     <  1:686121/151‑1 (MQ=255)
      agGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAAc                                                                                                                                    <  1:80018/151‑1 (MQ=255)
            aTACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCAt                                                                                                                              >  2:398470/1‑151 (MQ=255)
              aCTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTg                                                                                                                            >  2:108477/1‑151 (MQ=255)
                  aCTAACTCTAAAGTGTTATTTTACATGCACTTACAATTGATAAAAGACAACATTCCCAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGca                                                                                                                          <  2:378608/149‑1 (MQ=255)
                    tAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGcaaca                                                                                                                       >  1:962248/1‑150 (MQ=255)
                     aaCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAAt                                                                                                                     <  1:158631/151‑1 (MQ=255)
                         ctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGTGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATg                                                                                                                    <  1:1067500/148‑1 (MQ=255)
                           aaaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGgcgc                                                                                                                <  2:1017998/150‑1 (MQ=255)
                                gttATTTTACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTTACACAGAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGTTTACTGCCGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTc                                                                                                           <  2:226677/148‑1 (MQ=255)
                                ggTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTCACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCg                                                                                                          >  2:633232/1‑151 (MQ=255)
                                   aTTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGAt                                                                                                       >  1:322732/1‑151 (MQ=255)
                                    ttttACATGCACTTACAATAAATGAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGGGAGAAATGTAAGGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGAt                                                                                                       <  2:695983/150‑1 (MQ=255)
                                       tACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGa                                                                                                     >  1:940633/1‑149 (MQ=255)
                                        aCATGCACTTAAAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAg                                                                                                   <  2:405400/150‑1 (MQ=255)
                                                ttACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg                                                                                           >  1:276565/1‑150 (MQ=255)
                                                  aCAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACAAATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAc                                                                                         <  2:18587/150‑1 (MQ=255)
                                                    aaTTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg                                                                                           <  2:334317/146‑1 (MQ=255)
                                                    aaTTGATAAAAGACAACATTCACAGTGTGTTTATTTGTAACACATAGGGGCGAGAAATGTAATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACa                                                                                        <  2:633629/149‑1 (MQ=255)
                                                         aTTAAAGACAACATTCCCAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATa                                                                                  <  1:19188/150‑1 (MQ=255)
                                                           tAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                >  2:555421/1‑150 (MQ=255)
                                                            aaaGACAACATTCACAGTTTGGTTATTTGTTACAAAGAGGGGCGAGCAATGTCATGACAGTGTGGGTGCGTTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaaa                                                                               <  1:970722/150‑1 (MQ=255)
                                                             aaGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAAt                                                                             >  1:564711/1‑151 (MQ=255)
                                                             aaGACAACATTCACAGTGTGGATATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGAGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                >  2:447957/1‑148 (MQ=255)
                                                               gACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATa                                                                            >  1:931830/1‑150 (MQ=255)
                                                                     aTTCACAGTGTGGTTATTTGTTCCACATAGGGGCGAGCAATGTCATGCCAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgca                                                                      <  2:211677/150‑1 (MQ=255)
                                                                     aTTCACAGTGGGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGAAAGTGTGGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcac                                                                     <  2:225916/151‑1 (MQ=255)
                                                                              gtgGTTTTTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGGGGGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt                                                                <  2:441625/147‑1 (MQ=255)
                                                                              gggtttATTTTTTCCAAATAGGGGCGAGAAATGAAAGGACAGTGTGGGTGCGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATc                                                             <  1:397543/146‑1 (MQ=255)
                                                                               tgGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGAATCtt                                                           >  2:824032/1‑151 (MQ=255)
                                                                                    aTTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCt                                                      >  1:109761/1‑151 (MQ=255)
                                                                                     ttttttACAAATGGGGGCGAGCAATGTCATGACGGTGTGGGGGGGTTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCt                                                      <  2:539702/150‑1 (MQ=255)
                                                                                      tttttccacaTAGGGGCGAGAAATGTCAGGACGGTGTAGGTGCGTTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                                     <  1:73423/144‑1 (MQ=255)
                                                                                      tttttACACAAAGGGGCGAGAAATGTCATGACGGGGTAGGTGCGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                                     <  2:812790/147‑1 (MQ=255)
                                                                                      ttgtacaacaggggggTGAGCAATTTCAGGCCAGTGTGGGTGGGTTTTTTTTCGGAAAAAACAAAAAAACCCGCCATTCCAACAATGGCGCAATCCGGATGAAGCCCCTATTACAAGGATAAAAATAACCGCACGCCGTATCTTCAGCtt                                                     <  1:228470/138‑1 (MQ=255)
                                                                                       ttttACAAATAGGGGCGAGAAATTTCATGACGGGGTAGGGGGGTTTACTGTTGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                                     <  2:271245/147‑1 (MQ=255)
                                                                                         ttACACATGGGGGCGAGAAATGTCATGACGGTGTGGGTGCGTTTACTGTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAtt                                                  <  1:571002/150‑1 (MQ=255)
                                                                                         ttACAAATGGGGGCGAGCAATGTAATGACGGTGTGGGTGGGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAt                                                   <  2:341131/149‑1 (MQ=255)
                                                                                           acacAGGGGGGAGAGAAATGCCATGACAGTGTGGGGGGGTTTACTGTGGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg                                                 <  1:352367/149‑1 (MQ=255)
                                                                                              aaTAGGGGCGGGAATTGTCGTAAAAGGGTAGGGGCGGTTTTTGTGGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg                                                 <  1:382148/145‑1 (MQ=255)
                                                                                               aTAGGGGCGAGCAATGACATGACAGTGTAGGTGCGGTTTCTGTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCt                                           <  2:339401/151‑1 (MQ=255)
                                                                                                tAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtt                                          >  2:1083078/1‑151 (MQ=255)
                                                                                                  ggggggAGAAATGTCATGAAGGGGTGGGTGGGTTTCTTTTGGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg                                                 <  1:91006/142‑1 (MQ=255)
                                                                                                  ggggCGAGAATGGTCATGACAGTGAAGGTGCGGTTACGGCCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTcc                                            <  1:531368/147‑1 (MQ=255)
                                                                                                   ggggaaGCAATGTCATGACAGTGTAGGGGGGGTTACTGTGGAAAAAAACAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCTTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttttt                                       <  1:182395/146‑1 (MQ=255)
                                                                                                   gggCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttttt                                       >  2:905733/1‑151 (MQ=255)
                                                                                                   gggCGAGAAATGTAATGACAGTGTGGGTGCGTTTACTGTCGTAAAAAAAAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCTTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttttt                                       <  1:412949/151‑1 (MQ=255)
                                                                                                    ggCGAGCAATGTCATGACAGTGTGGGTGCGTTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttttt                                      <  1:261342/151‑1 (MQ=255)
                                                                                                      cGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                                    >  1:436027/1‑151 (MQ=255)
                                                                                                       gaaaaaTGTCATAACAGTGTGGGGGGGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttc                                     <  1:633140/145‑1 (MQ=255)
                                                                                                         gCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttt                                 >  1:86513/1‑151 (MQ=255)
                                                                                                           aaTTTCATAACAGGTAGGGGGGGTTTCCTGGCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttc                                     <  1:1082946/145‑1 (MQ=255)
                                                                                                           aaTGTCATGACGGTGTAGGTGGGTTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttttc                               <  2:417568/151‑1 (MQ=255)
                                                                                                            aTGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttttc                               >  1:796687/1‑150 (MQ=255)
                                                                                                                 aTGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt                         >  1:23885/1‑151 (MQ=255)
                                                                                                                    gaagggagggTGCGTTTACTGTCGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTc                          <  1:847851/139‑1 (MQ=255)
                                                                                                                    aCAGTGTGGGTGGGTTTATTGTCGTGAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCa                              <  1:480577/143‑1 (MQ=255)
                                                                                                                                 ggTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa          >  1:796545/1‑150 (MQ=255)
                                                                                                                                  gTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaaa         >  2:739690/1‑150 (MQ=255)
                                                                                                                                      aCTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACg       <  2:404618/149‑1 (MQ=255)
                                                                                                                                       cTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgct    >  1:958867/1‑151 (MQ=255)
                                                                                                                                         gTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgc  >  2:636951/1‑151 (MQ=255)
                                                                                                                                           cGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgc  <  2:538076/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                        
AAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGC  >  NC_000913/4287569‑4287855

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: