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

AAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGC  >  NC_000913/4287569‑4287855
                                                                                                                                                       |                                                                                                                                        
aaaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAAt                                                                                                                                          <  2:1093185/151‑1 (MQ=255)
aaaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAAt                                                                                                                                          <  1:87454/151‑1 (MQ=255)
aaaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAAt                                                                                                                                          >  1:471565/1‑151 (MQ=255)
 aaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATa                                                                                                                                         <  1:475353/151‑1 (MQ=255)
     gagGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATaaaaa                                                                                                                                     <  1:686121/151‑1 (MQ=255)
      agGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAAc                                                                                                                                    <  1:80018/151‑1 (MQ=255)
            aTACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCAt                                                                                                                              >  2:398470/1‑151 (MQ=255)
              aCTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTg                                                                                                                            >  2:108477/1‑151 (MQ=255)
                  aCTAACTCTAAAGTGTTATTTTACATGCACTTACAATTGATAAAAGACAACATTCCCAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGca                                                                                                                          <  2:378608/149‑1 (MQ=255)
                    tAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGcaaca                                                                                                                       >  1:962248/1‑150 (MQ=255)
                     aaCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAAt                                                                                                                     <  1:158631/151‑1 (MQ=255)
                         ctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGTGCGAGAAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATg                                                                                                                    <  1:1067500/148‑1 (MQ=255)
                           aaaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGgcgc                                                                                                                <  2:1017998/150‑1 (MQ=255)
                                gttATTTTACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTTACACAGAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGT‑TTACTGCCGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTc                                                                                                           <  2:226677/148‑1 (MQ=255)
                                ggTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTCACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCg                                                                                                          >  2:633232/1‑151 (MQ=255)
                                   aTTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGAt                                                                                                       >  1:322732/1‑151 (MQ=255)
                                    ttttACATGCACTTACAATAAATGAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGGGAGAAATGTAAGGACAGTGTAGGTGCGG‑TTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGAt                                                                                                       <  2:695983/150‑1 (MQ=255)
                                       tACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGa                                                                                                     >  1:940633/1‑149 (MQ=255)
                                        aCATGCACTTAAAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAg                                                                                                   <  2:405400/150‑1 (MQ=255)
                                                ttACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg                                                                                           >  1:276565/1‑150 (MQ=255)
                                                  aCAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACAAATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAc                                                                                         <  2:18587/150‑1 (MQ=255)
                                                    aaTTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg                                                                                           <  2:334317/146‑1 (MQ=255)
                                                    aaTTGATAAAAGACAACATTCACAGTGTGTTTATTTGTAACACATAGGGGCGAGAAATGTAATGACAGTGTAGGTGCGG‑TTACTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACa                                                                                        <  2:633629/149‑1 (MQ=255)
                                                         aTTAAAGACAACATTCCCAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATa                                                                                  <  1:19188/150‑1 (MQ=255)
                                                           tAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                >  2:555421/1‑150 (MQ=255)
                                                            aaaGACAACATTCACAGTTTGGTTATTTGTTACAAAGAGGGGCGAGCAATGTCATGACAGTGTGGGTGCGT‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaaa                                                                               <  1:970722/150‑1 (MQ=255)
                                                             aaGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAAt                                                                             >  1:564711/1‑151 (MQ=255)
                                                             aaGACAACATTCACAGTGTGGATATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGAGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                >  2:447957/1‑148 (MQ=255)
                                                               gACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATa                                                                            >  1:931830/1‑150 (MQ=255)
                                                                     aTTCACAGTGTGGTTATTTGTTCCACATAGGGGCGAGCAATGTCATGCCAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgca                                                                      <  2:211677/150‑1 (MQ=255)
                                                                     aTTCACAGTGGGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGAAAGTGTGGGTGCGG‑TTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcac                                                                     <  2:225916/151‑1 (MQ=255)
                                                                              gtgGTTTTTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGGGGGG‑TTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt                                                                <  2:441625/147‑1 (MQ=255)
                                                                              gggtttATTTTTTCCAAATAGGGGCGAGAAATGAAAGGACAGTGTGGGTGCGG‑TTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATc                                                             <  1:397543/146‑1 (MQ=255)
                                                                               tgGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGAATCtt                                                           >  2:824032/1‑151 (MQ=255)
                                                                                    aTTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCt                                                      >  1:109761/1‑151 (MQ=255)
                                                                                     ttttttACAAATGGGGGCGAGCAATGTCATGACGGTGTGGGGGGGT‑TTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCt                                                      <  2:539702/150‑1 (MQ=255)
                                                                                      tttttccacaTAGGGGCGAGAAATGTCAGGACGGTGTAGGTGCGT‑TTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                                     <  1:73423/144‑1 (MQ=255)
                                                                                      tttttACACAAAGGGGCGAGAAATGTCATGACGGGGTAGGTGCGG‑TTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                                     <  2:812790/147‑1 (MQ=255)
                                                                                      ttgtacaacaggggggTGAGCAATTTCAGGCCAGTGTGGGTGGGT‑TTTTTTTCGGAAAAAACAAAAAAACCCGCCATTCCAACAATGGCGCAATCCGGATGAAGCCCCTATTACAAGGATAAAAATAACCGCACGCCGTATCTTCAGCtt                                                     <  1:228470/138‑1 (MQ=255)
                                                                                       ttttACAAATAGGGGCGAGAAATTTCATGACGGGGTAGGGGGGT‑TTACTGTTGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                                     <  2:271245/147‑1 (MQ=255)
                                                                                         ttACACATGGGGGCGAGAAATGTCATGACGGTGTGGGTGCGT‑TTACTGTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAtt                                                  <  1:571002/150‑1 (MQ=255)
                                                                                         ttACAAATGGGGGCGAGCAATGTAATGACGGTGTGGGTGGGG‑TTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAt                                                   <  2:341131/149‑1 (MQ=255)
                                                                                           acacAGGGGGGAGAGAAATGCCATGACAGTGTGGGGGGGT‑TTACTGTGGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg                                                 <  1:352367/149‑1 (MQ=255)
                                                                                              aaTAGGGGCGGGAATTGTCGTAAAAGGGTAGGGGCGG‑TTTTTGTGGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg                                                 <  1:382148/145‑1 (MQ=255)
                                                                                               aTAGGGGCGAGCAATGACATGACAGTGTAGGTGCGG‑TTTCTGTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCt                                           <  2:339401/151‑1 (MQ=255)
                                                                                                tAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtt                                          >  2:1083078/1‑151 (MQ=255)
                                                                                                  ggggggAGAAATGTCATGAAGGGGTGGGTGGGT‑TTCTTTTGGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg                                                 <  1:91006/142‑1 (MQ=255)
                                                                                                  ggggCGAGAATGGTCATGACAGTGAAGGTGCGG‑TTACGGCCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTcc                                            <  1:531368/147‑1 (MQ=255)
                                                                                                   ggggaaGCAATGTCATGACAGTGTAGGGGGGG‑TTACTGTGGAAAAAAACAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCTTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttttt                                       <  1:182395/146‑1 (MQ=255)
                                                                                                   gggCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttttt                                       >  2:905733/1‑151 (MQ=255)
                                                                                                   gggCGAGAAATGTAATGACAGTGTGGGTGCGT‑TTACTGTCGTAAAAAAAAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCTTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttttt                                       <  1:412949/151‑1 (MQ=255)
                                                                                                    ggCGAGCAATGTCATGACAGTGTGGGTGCGT‑TTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttttt                                      <  1:261342/151‑1 (MQ=255)
                                                                                                      cGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                                    >  1:436027/1‑151 (MQ=255)
                                                                                                       gaaaaaTGTCATAACAGTGTGGGGGGGG‑TTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttc                                     <  1:633140/145‑1 (MQ=255)
                                                                                                         gCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttt                                 >  1:86513/1‑151 (MQ=255)
                                                                                                           aaTTTCATAACAG‑GTAGGGGGGGTTTCCTGGCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttc                                     <  1:1082946/145‑1 (MQ=255)
                                                                                                           aaTGTCATGACGGTGTAGGTGGGT‑TTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttttc                               <  2:417568/151‑1 (MQ=255)
                                                                                                            aTGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttttc                               >  1:796687/1‑150 (MQ=255)
                                                                                                                 aTGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt                         >  1:23885/1‑151 (MQ=255)
                                                                                                                    gaagggagggTGCGT‑TTACTGTCGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTc                          <  1:847851/139‑1 (MQ=255)
                                                                                                                    aCAGTGTGGGTGGGT‑TTATTGTCGTGAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCa                              <  1:480577/143‑1 (MQ=255)
                                                                                                                                 gg‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa          >  1:796545/1‑150 (MQ=255)
                                                                                                                                  g‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaaa         >  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)
                                                                                                                                                       |                                                                                                                                        
AAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGG‑TTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGC  >  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: ‑