Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077794,291,3000GA27.5% 129.3 / 28.1 69intergenic (‑339/‑54)acs/nrfAacetate‑‑CoA ligase/ammonia‑forming nitrite reductase cytochrome c552 subunit
Reads supporting (aligned to +/- strand):  ref base G (42/8);  new base A (0/19);  total (42/27)
Fisher's exact test for biased strand distribution p-value = 4.80e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTG  >  NC_007779/4291051‑4291549
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         
aGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTg                                                                                                                                                                                                                                                           <  1:1096083/250‑1 (MQ=255)
 gATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGa                                                                                                                                                                                                                                                          <  2:568697/250‑1 (MQ=255)
    cACGCAAACGTAGAATTTTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAACGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGCGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCTTGACAGTGTAGGTGCGGTTACTGTCGTGcaac                                                                                                                                                                                                                                                       >  2:416280/1‑249 (MQ=255)
    cACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGaaaa                                                                                                                                                                                                                                                       >  1:1087841/1‑250 (MQ=255)
           aaGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                                                                                                                                                                                                                                                <  1:1074856/250‑1 (MQ=255)
             gTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaaaa                                                                                                                                                                                                                                              <  2:702946/250‑1 (MQ=255)
                     gtgCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCAt                                                                                                                                                                                                                                      >  2:133974/1‑250 (MQ=255)
                      tgCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCAtt                                                                                                                                                                                                                                     >  2:877980/1‑250 (MQ=255)
                           aTAATCGGCAGGTAAATTTTTAAAGGTAAGCTACAAAAGTTGATTAATTCTTTGAGGAACATTCAGTTATGGATGCGGTTAAAAAAGAGGAATATATTAAATAACTCTAAAGTGGTATTTTACATACACTAACAATTAATAAAAGACATCATTCGCATTGTGGTTATTTGTAACACAAAGGGCAAAGAAATGTCAGACGTGTGTAGGTGCTGTTACGATGTTAAAAAAAAATAAAAACCGCCATTGcaac                                                                                                                                                                                                                                <  2:895785/250‑1 (MQ=255)
                           aTAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGcaac                                                                                                                                                                                                                                >  1:602281/1‑250 (MQ=255)
                                 ggCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGgc                                                                                                                                                                                                                          <  1:783391/250‑1 (MQ=255)
                                  gCAGGGTAATTTTTGAAGGTCAGAAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTAAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGgcg                                                                                                                                                                                                                         <  1:856176/250‑1 (MQ=255)
                                               tGAAGTTAAGCAAAAAAAGTTGATAAATTTTTTGAGGAACATGCAGTTATGCATGTTGTTAAAAAAGAGGAAAATACTGACAAACTATAAAGTGTTATTTTACATGCACTTAAAATTGATAAAAAACAACATTCAAAGTGTGTTTTTTTGTTACAAATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGTTTACTGTCGTAAAAAAAAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGaa                                                                                                                                                                                                            <  1:766337/250‑1 (MQ=255)
                                                     tCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCt                                                                                                                                                                                                      >  1:1073961/1‑250 (MQ=255)
                                                           aaaaaaGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCATTTACAATTGATTAAAGAAAACATTCACAGTGTGGTTATTTTTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACa                                                                                                                                                                                                <  1:206384/248‑1 (MQ=255)
                                                                          ttCTTTGAGGAACATGCAGTTATGCATGTTTTTAAAAAAGAGGAAGATACTGATTAACTCTAAAGTGTTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTTTTTGTTACACATAGGGGCGAGCAATTTCATGACAGTGTAGGTGCGTTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAac                                                                                                                                                                                 <  2:733882/250‑1 (MQ=255)
                                                                          ttCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAac                                                                                                                                                                                 >  1:142859/1‑250 (MQ=255)
                                                                          ttCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAac                                                                                                                                                                                 >  1:1022091/1‑250 (MQ=255)
                                                                           tttttGAGAAACATGCAGTTATGCTTGTTTTAAAAAAAAAGGAAGATACTGACTAACTAAAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACTTGCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTAATGACAGTGGAGGTGCGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacg                                                                                                                                                                                <  2:813146/248‑1 (MQ=255)
                                                                            cTTTGAGAAACATGAAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTTACACAGGGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgc                                                                                                                                                                               <  2:1062552/250‑1 (MQ=255)
                                                                             tttGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgca                                                                                                                                                                              >  2:866423/1‑250 (MQ=255)
                                                                             tttGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgca                                                                                                                                                                              >  1:577403/1‑250 (MQ=255)
                                                                              ttGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcac                                                                                                                                                                             >  1:123842/1‑250 (MQ=255)
                                                                              ttGAGGAACATGCAGTTAGGCATGCTGTTGAAAAAGAGAAAAATATTGACTAATTCTAAAGTGGTATTTTACATGCATTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcac                                                                                                                                                                             <  2:972409/250‑1 (MQ=255)
                                                                                 agaaaCATGCAGTTTTGCATGCTGTTGAAAAAGAGAAAGATACTTACTAACTCTAAAGGGGTATTTTACATGCCTTTACAATTAATTAAAGACAACATTCACAGTGTGTTTTTTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGcc                                                                                                                                                                          <  1:355709/247‑1 (MQ=255)
                                                                                 aGGAACATGCATTTATGCATGCTTTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGTTTTTTTACATGCATTTACAATTGATTAAAGACACCATTCACATTGTGTTTTTTTTTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTCCTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGcc                                                                                                                                                                          <  1:523239/250‑1 (MQ=255)
                                                                                 aGGAACATGAAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTAAAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGcc                                                                                                                                                                          <  1:114147/250‑1 (MQ=255)
                                                                                  ggAACATGAAGTTATGCATGCTGTTGAAAAAAAGGAAGATACTAACTAACTCTAAGGTGGTATTTTACATGCATTTAAAATTGATAAAAGACAACATTCACAGTGTGTTTATTTGTAACACATAGGGGCGAGCAATGTAATGACAGTGTAGGGGCGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCg                                                                                                                                                                         <  2:892413/250‑1 (MQ=255)
                                                                                     aCAGGCAGTTATGCATGTTGTGAAAAAAGAGAAAGATACTGACTAACTAAAAAGTGGTATTTTACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGAAATTTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTAt                                                                                                                                                                      <  2:594520/250‑1 (MQ=255)
                                                                                      cATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATc                                                                                                                                                                     >  1:389826/1‑250 (MQ=255)
                                                                                       atacAGTAATGCATGTTGTTAAAAAAAAGGACGACACTTACGAACTCTAAAGATGTATTTTAGATCCACTTACAATTTATAAAGGACAACATTCAGTGTGTGGTTTTTTGTTAACAATGGGGGCGGGCAATATCATGCCATAGTGGGTAAGGTTCTTGGCTAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCt                                                                                                                                                                    <  2:835082/247‑1 (MQ=255)
                                                                                       aTGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCt                                                                                                                                                                    >  2:279379/1‑250 (MQ=255)
                                                                                          catttAGGCATGTGGTTAAAAAAGAGGAAGATACTGACTAACTCTAAGGTGTTATTTTACATCCATTTAAAATTGATAAAAAACAACATACAGGGTGTGTTTATTTGTTAAAAATAGGGGCGAAAAATGCAATGACGTTGTAGGTGCGTTTCCTGTGGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCa                                                                                                                                                                 <  2:926801/247‑1 (MQ=255)
                                                                                            tttATGCATGTGTGTAAAAAAGAGGAAGAACCTGACTAATTAAAAGGGGGTATTTAAAATGCATTTAAAATGGATAAAAGAAAACATGAACAGGGTGGTTATTGGTTACACAGGGGGGCGAGAAATGTCATGACAGTGTGGGGGCGTTTACTGTCGGAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGc                                                                                                                                                               <  2:197167/249‑1 (MQ=255)
                                                                                                  catgcTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg                                                                                                                                                         >  1:542784/1‑250 (MQ=255)
                                                                                                    tgcTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAt                                                                                                                                                       >  2:785849/1‑250 (MQ=255)
                                                                                                    tgcTGTTAAAAAAGAGGAAGATACTGAATAACTCTAAAGTGGTATTTTACATGCATTTAAAATTGATAAAAGAAAAAATTCACAGTGTGTTTATTTTTTACACATGGGGGCGAGAAATGTCATGACAGTGTAGGGGCGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAt                                                                                                                                                       <  1:133974/250‑1 (MQ=255)
                                                                                                       ctttaaaaaaaGAGAAATATACTGATTACCTTTAAAGTGTTTTTTTACATGCCTTTACAATTGATAAAAAACAACATTCACAGTGTGTTTTTTTTTTACAAATGGGGGCGACAAATGTCAGGACAGTGTAGGTGGGTTTCTTGTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTc                                                                                                                                                     <  1:95692/244‑1 (MQ=255)
                                                                                                           agAAAAAGAGGAAAATTCTCACAAACTTAAAGGTTGTTTTTTAATTGCTCTTAAATTTAATAAAAGAAAACTTTCACATGGTGTTTTTTTTTTAAAAAAGGGGGCGACAAATTTCATGACGGTGTGGGTGGGTTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttt                                                                                                                                                 <  2:894862/249‑1 (MQ=255)
                                                                                                            gAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttt                                                                                                                                                >  1:714486/1‑250 (MQ=255)
                                                                                                              aaaaGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttttt                                                                                                                                              >  1:511247/1‑250 (MQ=255)
                                                                                                                aaGAGGAAAATACTTACTAAACAAAAAGTGGTTTTTTAAATGAACTTAAAATTGAAAAAAAAAAACATTCAAAGGGGTGTTATTTTTAACACAGGGGGGCAAACAATTACATGACAGTGTGGGTGCGGTTACTGACGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                                                                                                                                            <  2:714486/250‑1 (MQ=255)
                                                                                                                     gAAAATACTGACTAACTCAAAGGTGGTATTTTACATGCACTTAAAATTGATAAAAGACAAAATTCACGGTGTGTTTTTTTTTTACACATGGGGGCAAAAAATGACATGACAGTGTAGGGGGGTTTACTGTAGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttttc                                                                                                                                       <  2:669482/250‑1 (MQ=255)
                                                                                                                      aaaatattaacaaaCTCAAAATTTGTATTTTACATCCATTTAAAATTAATAAAACAAAAATTTCACTGTGTTTTTTTTTTTTACAAATGGGGGCGAAAAATTTCATGACGGTGTAGGTGCGTTTACTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCa                                                                                                                                      <  2:1087841/238‑1 (MQ=255)
                                                                                                                         tatataacaatctcaAAAGTGGTTTTTTACATCCATTTAAAATTAATAAAAAAAAAAATCCAGAGGGTGTTTATTTTTACAAAAAGGGGGCGACAAATTAAATACCAGTGAAGGGGGGTTTCTGTTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACtt                                                                                                                                   <  2:1123469/239‑1 (MQ=255)
                                                                                                                          tACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTc                                                                                                                                  >  1:180307/1‑250 (MQ=255)
                                                                                                                           aCTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCt                                                                                                                                 >  2:118737/1‑250 (MQ=255)
                                                                                                                           aCTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTATTTTTCACTTCt                                                                                                                                 >  2:321059/1‑250 (MQ=255)
                                                                                                                             tGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGt                                                                                                                               >  2:347211/1‑250 (MQ=255)
                                                                                                                              gACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGtt                                                                                                                              >  1:728190/1‑250 (MQ=255)
                                                                                                                              gACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGtt                                                                                                                              >  1:709483/1‑250 (MQ=255)
                                                                                                                                   aCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGc                                                                                                                         >  1:939354/1‑250 (MQ=255)
                                                                                                                                    ctctAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCt                                                                                                                        >  2:715214/1‑250 (MQ=255)
                                                                                                                                        aaaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAAc                                                                                                                    >  2:1011167/1‑250 (MQ=255)
                                                                                                                                        aaaGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAAc                                                                                                                    >  1:807/1‑250 (MQ=255)
                                                                                                                                                aTTTTAAATGCATTTAAAATGAATAAAAGACAACATTCACAGTGTGTTTTTTTGTTACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGAAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgct                                                                                                            <  1:359486/250‑1 (MQ=255)
                                                                                                                                                   ttACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgct                                                                                                         >  1:521972/1‑250 (MQ=255)
                                                                                                                                                     aCATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTcc                                                                                                       >  1:67944/1‑250 (MQ=255)
                                                                                                                                                     aCATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTcc                                                                                                       >  2:1098189/1‑250 (MQ=255)
                                                                                                                                                                          aaGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGa                                                                                  >  2:947121/1‑250 (MQ=255)
                                                                                                                                                                            gACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAg                                                                                >  1:506525/1‑250 (MQ=255)
                                                                                                                                                                                     cacaGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATg                                                                       >  1:196529/1‑250 (MQ=255)
                                                                                                                                                                                        gggggtgTTATTTGTTCAAAAGAGGGGAAAAAAATGACATAACGGTGGAGGGGGGGTGACGGAAGAAAAAAAAAAAAAAAACCGCAATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGaaa                                                                    <  2:389826/244‑1 (MQ=255)
                                                                                                                                                                                              gTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTg                                                              >  1:207891/1‑250 (MQ=255)
                                                                                                                                                                                              gTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTg                                                              >  1:208552/1‑250 (MQ=255)
                                                                                                                                                                                                       tACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCgcagc                                                     >  1:105057/1‑250 (MQ=255)
                                                                                                                                                                                                          acaTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATc                                                  >  1:424626/1‑250 (MQ=255)
                                                                                                                                                                                                          acaTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATc                                                  >  2:379563/1‑250 (MQ=255)
                                                                                                                                                                                                           caTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATcc                                                 >  2:242150/1‑250 (MQ=255)
                                                                                                                                                                                                              tgggggCGAGAATTGTCATGACAGTGTGGGGGAGTTTACTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCg                                               <  1:877980/248‑1 (MQ=255)
                                                                                                                                                                                                                    cGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAat                                         >  2:71813/1‑250 (MQ=255)
                                                                                                                                                                                                                     gaacAAAGTCATAACAGTGTGGGTGCGGTTAGTGGAGTAAAAACAAAAAAAAACCGCCATTCCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACCGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAata                                        <  2:915451/247‑1 (MQ=255)
                                                                                                                                                                                                                                             tGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCgg                 >  1:328456/1‑250 (MQ=255)
                                                                                                                                                                                                                                               cGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGaa               >  2:389520/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         tCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCACTCCTGGAAAGCCACCTCGGAACAGTCagag      >  2:907284/1‑250 (MQ=255)
                                                                                                                                                                                                                                                             gAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCgtg  >  1:681274/1‑250 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         
AGATCACGCAAAAGTAGAATTGTGCAAATAAACGGCAGGGTAATTTTTGAAGGTCAGCAACAAAAGTTGATTAATTCTTTGAGGAACATGCAGTTATGCATGCTGTTGAAAAAGAGGAAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGGAAGCGAAGAATGAAACCTTTGCCCCGCAGCATCCCGATCAATATCTCTCCTGGAAAGCCACCTCGGAACAGTCAGAGCGTG  >  NC_007779/4291051‑4291549

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