Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,7850TA47.4% 31.0 / 23.1 38intergenic (‑399/+149)sufA/ydiHFe‑S cluster assembly protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (13/7);  new base A (0/18);  total (13/25)
Fisher's exact test for biased strand distribution p-value = 1.59e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTT‑AAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTC  >  NC_000913/1764637‑1764929
                                                                                                                                                     |                                                                                                                                                
tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTT‑AAGTCTGTGCAACACAGGCTTAAGGGCTTCTACACCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATa                                                                                                                                                 >  1:168137/1‑150 (MQ=255)
  agagATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTT‑AAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaa                                                                                                                                                >  1:661219/1‑149 (MQ=255)
   gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTT‑AAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGc                                                                                                                                             <  1:346560/151‑1 (MQ=255)
         ggTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTT‑AAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAccc                                                                                                                                       <  1:916520/151‑1 (MQ=255)
                    gCAGCTTTCAAAAAAAAGAAACCCCTT‑AAGTCTGGGCGACACAGGCTAAGGGGTTTCTCCCCCACCCGGCGCTTATCTCCGGCACTCTCAGTGGTTTAGCTCTGGAGGGGGCGATAAGAATAATCTCATAAATCTACCCCGCCGtttt                                                                                                                               <  1:233617/148‑1 (MQ=255)
                         gTTCAGAAATAAGAAAACCCTT‑AAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACaa                                                                                                                        >  2:347889/1‑150 (MQ=255)
                            cAGAAATAAGAAAACCCTT‑AAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTg                                                                                                                    >  2:607369/1‑151 (MQ=255)
                                 aaaaaaaaaCCCTTAAAGTCTGTGCGACACAGTTTTAAGGTTTTCTCCCCCACCCGGCGTTTATCTCCGGCACTCTCAGTGGTTTAGCTCTTGAGGGGGCAAAAAGAATAATCTCAAAAGGCAACCCCGCCGTTTTAACACAAACTGCGAt                                                                                                                <  2:561258/146‑1 (MQ=255)
                                    aGAAAACCCTT‑AAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGt                                                                                                            >  2:383229/1‑151 (MQ=255)
                                          cccTT‑AAGTCTGTGCGACACAGGCTTAAGGTTTTCTCCCCCATCCGGCGCTTATCCCCGGCACTCTCAGTGGCTTAGCTCTGGAAGGGGCGATAAGAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtattat                                                                                                       <  2:218886/150‑1 (MQ=255)
                                           ccTT‑AAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt                                                                                                     >  1:672990/1‑151 (MQ=255)
                                                aaGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg                                                                                                  >  2:783298/1‑150 (MQ=255)
                                                      gggCGACACGGGTTAAGGGGTTTCCCCCCCACCGGGCGTTTATCCCCGCCCCTCTCAGTGGCTTAGCTTTTGAGGGGGCAAAAAAAATAATCTAATAAAGTTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACa                                                                                              <  1:225992/146‑1 (MQ=255)
                                                               aGGATTAAGGGTTTCTCCCCCACCCGGCGATTATCTCGGGCTTTCCATGTGGTTTAGCTTTTAAAGGGGCGATAAAAAAAACCTAAAAACGCGATCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt                                                                                  <  2:771549/151‑1 (MQ=255)
                                                                ggCTTAAGGGTTTCTCCCCCACCCGGCTCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTAAAGGGGAAAAAAAAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCgg                                                                                   <  2:28069/149‑1 (MQ=255)
                                                                ggCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt                                                                                  >  1:570071/1‑150 (MQ=255)
                                                                 gCTAAAGGGTTGCTCCCCCACCCGGCGCTTTCCGCGGGCCCTCACAGTGGCTTACCTCTAGAAGGGGCGAAAAAAAAAATCTCAAAAAACAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAg                                                                                <  2:860320/151‑1 (MQ=255)
                                                                     aaGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATa                                                                             >  2:519577/1‑150 (MQ=255)
                                                                       ggttttCTCCCCCACCCGGCCCCTACCCCGGGCCCTCTCAGGGGTTTAGCTCTTGAAGGGGCGATAAGAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATaa                                                                            <  2:627119/146‑1 (MQ=255)
                                                                           ttcccccccACCCGGCGTTTATCCCCGGCACTCTCAGTGGTTTAGCTCTTGAGGGGGCGAAAAAAAAAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGt                                                                       <  2:592478/145‑1 (MQ=255)
                                                                                        ggCGCTTTTCTCCGGCCCCCTAAGGGGTTTGGCTTTAAAGGGGGAGAAAAAAATAATCTCAAAATAATACACCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaagaag                                                              <  2:101597/146‑1 (MQ=255)
                                                                                        ggCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCAt                                                         >  2:373507/1‑151 (MQ=255)
                                                                                               aTCTCCGCCACTCTCAGTGGTTTAGCTCTTAAAGGGGCAAAAAAAAAAATCACAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAAc                                                    <  2:943572/149‑1 (MQ=255)
                                                                                                 cgccGGCCCTCTCAGGGTTTTAGCTCTTGAGGGGGAAAAAAGAATAATCTAATAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCt                                                  <  2:640826/147‑1 (MQ=255)
                                                                                                           ctcGGTGGCTTAGCTTTTAAGGGGACAAAAAAAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAgg                                         <  1:623239/148‑1 (MQ=255)
                                                                                                              aGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATc                                   >  2:112430/1‑151 (MQ=255)
                                                                                                                 gtcTTAGCTCTTGAGGGGGGAATAAAAATATTCTAAAAAATCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAg                                <  1:934141/149‑1 (MQ=255)
                                                                                                                 ggCTTAGCTCTTAAAGGGGCAATAAGAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCa                                 <  2:688954/150‑1 (MQ=255)
                                                                                                                    ttgactctTGAAGGGGCAAAAAAAATAATCAAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAgg                               <  2:651168/145‑1 (MQ=255)
                                                                                                                      aGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCg                            >  2:177595/1‑150 (MQ=255)
                                                                                                                        ctttagAAGGGCCAATAAAAATAGCCGAATAACGCAATCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGt                              <  2:313905/141‑1 (MQ=255)
                                                                                                                         tgttGAAGGGGAGAAAATAATAATCACAAAAACCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCg                            <  1:243090/145‑1 (MQ=255)
                                                                                                                         tctTGAGGGGGCAAAAAAAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCg                            <  2:750888/147‑1 (MQ=255)
                                                                                                                              aaGGGGCAAAAAGAATAATCAAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTgg                      <  1:447950/148‑1 (MQ=255)
                                                                                                                                      aTAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATTGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCgg           >  1:96534/1‑151 (MQ=255)
                                                                                                                                        aaGAATAATCTCAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGc          <  2:334193/150‑1 (MQ=255)
                                                                                                                                          aaataatCTAAAAAACCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCAtt       <  2:62519/150‑1 (MQ=255)
                                                                                                                                                    aTAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTc  <  2:802104/146‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTT‑AAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTC  >  NC_000913/1764637‑1764929

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