Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,287,7190TA44.1% 75.2 / 47.9 69intergenic (‑348/‑45)acs/nrfAacetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base T (30/8);  new base A (0/30);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 4.67e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTG  >  NC_000913/4287569‑4287869
                                                                                                                                                      |                                                                                                                                                       
aaaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                                                                                                                                         >  1:822252/1‑151 (MQ=255)
 aacagagGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                                                                                                                                         >  1:834697/4‑150 (MQ=255)
 aaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATa                                                                                                                                                        <  1:180782/151‑1 (MQ=255)
    agagGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaaaa                                                                                                                                                     <  2:568908/151‑1 (MQ=255)
      agGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAAc                                                                                                                                                  <  2:1033666/151‑1 (MQ=255)
           gATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGcc                                                                                                                                              >  2:139888/1‑150 (MQ=255)
            aTACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCa                                                                                                                                             <  1:1167402/150‑1 (MQ=255)
                   cTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGcaac                                                                                                                                      >  1:355737/1‑150 (MQ=255)
                     aaCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAAt                                                                                                                                   >  1:43649/1‑151 (MQ=255)
                     aaCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAAt                                                                                                                                   >  1:933581/1‑151 (MQ=255)
                      aCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATg                                                                                                                                  >  1:285507/1‑151 (MQ=255)
                      aCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATg                                                                                                                                  >  1:285522/1‑151 (MQ=255)
                          tAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGgcgc                                                                                                                              >  1:514693/1‑151 (MQ=255)
                           aaaGGGGTTTTTTACATGCACTTAAAATTGATAAAAGACAACATCCACAGTGTGGTTTTTTGTTACACATGGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCa                                                                                                                             <  1:752697/151‑1 (MQ=255)
                            aaGGGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCaa                                                                                                                            <  2:928129/151‑1 (MQ=255)
                               tGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTc                                                                                                                         >  1:707617/1‑151 (MQ=255)
                                ggTATTTTACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGTTTTTTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGGGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCg                                                                                                                        <  2:1005239/151‑1 (MQ=255)
                                  tATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCgg                                                                                                                       >  2:143442/1‑150 (MQ=255)
                                  tATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCgg                                                                                                                       >  2:529364/1‑150 (MQ=255)
                                            gCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGcccc                                                                                                             <  2:1099379/150‑1 (MQ=255)
                                            gCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCt                                                                                                            <  2:153258/151‑1 (MQ=255)
                                                ttACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg                                                                                                         >  2:201426/1‑150 (MQ=255)
                                                  aCAATTGATTAAAGACAACATTCACAGTGTGTTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGTTTACTGTCGTGAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATg                                                                                                         <  1:793057/148‑1 (MQ=255)
                                                  aCAATTGATTAAAGACAACATTCACAGTGTGTTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGa                                                                                                        <  1:417395/149‑1 (MQ=255)
                                                    aaTTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACa                                                                                                      >  2:783775/1‑149 (MQ=255)
                                                         aTTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTAAGGACAGTGTAGGTGCGGTTACTGTCGTAAAAAAAAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaa                                                                                               <  1:294234/151‑1 (MQ=255)
                                                               gACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATa                                                                                          >  1:105361/1‑150 (MQ=255)
                                                                  aacaTTCACAGTGTGTTTTTTTGTTACACATGGGGGCAAGAAATGTCATGACGTGGTAGGTGCGTTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacg                                                                                      <  1:916563/151‑1 (MQ=255)
                                                                  aacaTTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacg                                                                                      >  2:621003/1‑151 (MQ=255)
                                                                  aacaTTCACAGTGTGGTTATTTGTTACACAAAGGGGCGAGCAATGTCATGACAGTGTGGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacg                                                                                      <  1:281444/151‑1 (MQ=255)
                                                                    caTTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgc                                                                                     >  2:206998/1‑150 (MQ=255)
                                                                     aTTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcac                                                                                   >  1:724157/1‑151 (MQ=255)
                                                                     aTTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGATACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgca                                                                                    >  1:171661/1‑150 (MQ=255)
                                                                        cacaGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGcc                                                                                >  2:249396/1‑151 (MQ=255)
                                                                           aGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt                                                                              >  1:718803/1‑150 (MQ=255)
                                                                             tgtgTTTTTTTTTTACACAAGGGGGCGAGAAATGTAATAACAGTGTGGGTGCGTTTACTGCCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt                                                                              <  2:1140422/148‑1 (MQ=255)
                                                                             tgtgGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATc                                                                           >  2:225370/1‑151 (MQ=255)
                                                                              gtgtTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTGGGTGCGTTTACTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTAt                                                                            <  2:1093088/149‑1 (MQ=255)
                                                                              gtgGTTTTTTTTTACACATAGGGGCGAACATTGTCATGCCAGTGTGGGTGCGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTAt                                                                            <  1:725051/149‑1 (MQ=255)
                                                                                  tttttttttACACATGGGGGCGAGAATTGTCATCACATTTTGGGTGGGTTTTCTGTCGAAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTc                                                                        <  2:613550/146‑1 (MQ=255)
                                                                                   tttttGTTACAAAAGGGGGCGAAAAATGCAAGGAAGGTGTGGGTGGGTTTACGGTGGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAg                                                                      <  1:381244/148‑1 (MQ=255)
                                                                                     tttGTTACAAATGGGGGCAAGAAATGTAAGGAAGGGGTAGGTGCGTTTACTGTAGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTc                                                                        <  1:537636/146‑1 (MQ=255)
                                                                                            cacaTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt                                                            >  1:273537/1‑151 (MQ=255)
                                                                                               aagggggggAGAAATGTCATGACGGTGTGGGGGGGTTTACTGTGGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTcc                                                          <  1:1000440/147‑1 (MQ=255)
                                                                                                  ggggCGAGCAATGTCATGACAGTGTGGGTGGGGTTACTGTGGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCt                                                         <  2:990545/148‑1 (MQ=255)
                                                                                                   gggCGAGAAATGTCAGGACAGTGTGGGGGCGTTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttttt                                                     <  2:43651/151‑1 (MQ=255)
                                                                                                   gggCGACAAATGTCATGCCAGTGTAGGGGGGTTTACTGGCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttttt                                                     <  2:328333/151‑1 (MQ=255)
                                                                                                      cGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                                                  >  2:891278/1‑151 (MQ=255)
                                                                                                       gAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctt                                                 >  1:191967/1‑151 (MQ=255)
                                                                                                       gAGCAATGTCATGACAGTGTAGGCGCGGTTACTGTCGTGCAAAACAATAAAACCCGCCAGTGCAACAATGGCGCAATTCGGATGAAGCACCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctt                                                 >  1:1131788/1‑151 (MQ=255)
                                                                                                        agaaaTGTCAGGACGGTGTGGGTGGGGTTACTTTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttt                                                <  1:42027/148‑1 (MQ=255)
                                                                                                        aacagtGTCATGAAAGTGTAGGTGCGTTTCCTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                                                  <  2:1064878/144‑1 (MQ=255)
                                                                                                             tgtaagaaCAGTGTGGGTGGGTTTACTGTGGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCa                                            <  1:344142/143‑1 (MQ=255)
                                                                                                                 aaaaCAGTGTAGGTGGGGTTACGGTGGTAAAAAAAAATAAAAACCCCCATCGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAACCGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACt                                          <  2:59946/145‑1 (MQ=255)
                                                                                                                    aCAGTGTGGGGGCGTTTCCTGTCGGAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTc                                        <  2:404806/147‑1 (MQ=255)
                                                                                                                       gtgtGGGTGGGTTTATTGTCGGAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCa                                  <  2:621315/150‑1 (MQ=255)
                                                                                                                         gtggtggCGGTTATTGTTGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTg                                      <  1:642377/138‑1 (MQ=255)
                                                                                                                         gtgggggggtttccTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCAc                                 <  2:451866/136‑1 (MQ=255)
                                                                                                                            ggtggagttatattCGTAAAAAAAAACAAAAACCCGCCATTGCAAGAATGGCGCAATTCGTATGAAGCCCCTATCAAAGGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCTTTTTTTCTTTTTCACTTCTGTTc                                   <  1:664053/132‑1 (MQ=255)
                                                                                                                            ggggCGTTTACTGTCGGAAAAAACAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTg                                      <  2:832782/138‑1 (MQ=255)
                                                                                                                            ggTGGGTTTACTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTc                                        <  1:687494/139‑1 (MQ=255)
                                                                                                                            ggTGCGTTTACTGTCGTAAAAAACAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGa                            <  2:1125613/151‑1 (MQ=255)
                                                                                                                                 ggTTACTGTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACaa                        <  2:4007/150‑1 (MQ=255)
                                                                                                                                    tACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACgg                    >  2:28162/1‑151 (MQ=255)
                                                                                                                                     actggcgaaaaaaaaaaaaaaaaCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACg                     <  2:982105/140‑1 (MQ=255)
                                                                                                                                     aCTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgc                   >  2:1177225/1‑151 (MQ=255)
                                                                                                                                      cTGTCGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACgg                    <  1:162786/149‑1 (MQ=255)
                                                                                                                                       tgtgttaaaaaaaaaaaaaACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgct                  <  1:206572/143‑1 (MQ=255)
                                                                                                                                        gtggtaaaaaaaaaaaaaaaCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACgg                    <  2:912908/141‑1 (MQ=255)
                                                                                                                                        gtcgaaaaaaaaaaaaaaaaCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctg                 <  1:31383/144‑1 (MQ=255)
                                                                                                                                         tCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgct               >  1:248562/1‑151 (MQ=255)
                                                                                                                                             aaaaaaaaaaaaaaaCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgct               <  1:972243/146‑1 (MQ=255)
                                                                                                                                             aaaaaaCAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgct               <  2:1179504/146‑1 (MQ=255)
                                                                                                                                              aaaaaaaaaaaaACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgc                <  1:987399/145‑1 (MQ=255)
                                                                                                                                              aaaaaaaaaaaaACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTcc             <  2:448971/148‑1 (MQ=255)
                                                                                                                                                aaaCAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCa         >  2:558526/1‑150 (MQ=255)
                                                                                                                                                  aCAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCaaa       >  1:484627/1‑150 (MQ=255)
                                                                                                                                                     aTAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCt   >  1:704832/1‑151 (MQ=255)
                                                                                                                                                      tAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTg  >  1:83578/1‑151 (MQ=255)
                                                                                                                                                      tAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATATTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTg  >  1:254126/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
AAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTG  >  NC_000913/4287569‑4287869

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