Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,7830CA44.7% 53.3 / 16.9 38intergenic (‑397/+151)sufA/ydiHFe‑S cluster assembly protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (13/8);  new base A (0/17);  total (13/25)
Fisher's exact test for biased strand distribution p-value = 6.20e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGG  >  NC_000913/1764640‑1764932
                                                                                                                                                |                                                                                                                                                      
gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGc                                                                                                                                                 <  2:461117/151‑1 (MQ=255)
  gaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTa                                                                                                                                               <  2:449456/151‑1 (MQ=255)
         aaaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgcc                                                                                                                                       <  1:648639/151‑1 (MQ=255)
           aaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGt                                                                                                                                     >  1:252252/1‑151 (MQ=255)
                 gCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAac                                                                                                                               >  1:702127/1‑151 (MQ=255)
                    gCGTACAAAAATAAAAAACCCCTTAAGTCTGTGCGACACGGGTTTAAGGGTTTCTCCCCCATCCGGCGCTTATCCCCGGCCCTCTCAGTGGTTTAGCTCTTGAGGGGGCGATAAGAATAATCTCATAAAGCTACCCCGCCGTTTTAAcaca                                                                                                                            <  1:804224/151‑1 (MQ=255)
                      gTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACaa                                                                                                                           >  2:328215/1‑150 (MQ=255)
                          aaaaaTAAGAAACCCCTTAAGTCTGTGCGACACAGGTTAAAGGTTTTCTCCCCCACCCGGCGTTTATCTCCGGCCCTCCCAGTGGCTTAGCTTTTAAAGGGGCGATAAGAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTg                                                                                                                       <  2:213763/148‑1 (MQ=255)
                              aTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAtt                                                                                                                  >  1:515778/1‑151 (MQ=255)
                                   aaacCCTTTAGGTTTGTGCGACACAGGCTTAGGGGTTTCTCCCCCATCCGGCCCTTATCTCCGGCCCTCTCAGTGGTTTAGCTTTTGAGGGGGCAATAAAAAAAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGta                                                                                                              <  1:769120/150‑1 (MQ=255)
                                    aaaCCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtatt                                                                                                            >  2:444818/1‑151 (MQ=255)
                                       ccttaaaGTTTGTGCGACACAGCCTAAGGGGTTTCTACCCCACCCGGCCTTTATCCCCGCCCCTCTCAGTGGCTGAGCTTTTAAAGGGGAAAAAAAAATAATCTAAAAAACCAAACCCGCCGTTTTAACACAAACTGCGATTAGTattatt                                                                                                         <  2:186995/146‑1 (MQ=255)
                                                       acacggaattAGGGGTTTCTCCCCCATCGGGCGTTTATCTCCGGCCCTCGCAGTGGCTTAGCTTTGGAAGGGGGAAAAAGAATAATCTAATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAg                                                                                          <  1:229951/142‑1 (MQ=255)
                                                          caGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGTGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCgg                                                                                      >  2:143797/1‑151 (MQ=255)
                                                            ggCTTAAGGTTTTCTCCCCCACCCGGCGCTTATCCCCGGCCCTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAAAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                                                                        <  1:81962/147‑1 (MQ=255)
                                                                  aGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATa                                                                                >  2:613835/1‑149 (MQ=255)
                                                                       ttCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGt                                                                          >  2:441376/1‑150 (MQ=255)
                                                                           cccccACCGGGCTCTTATCCCCGCCCCTCTCAGTGGTTTAGCTTTTGAGGGGGCGAAAAAAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTa                                                                      <  1:253420/149‑1 (MQ=255)
                                                                            ccccACCCGTCACTAATCCCCGCCACCCTCAGTGGCTGAGCTTTAGAAGGGGCGTAAAAAATAATCAAAAAAAGTTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAACATCAGGCGGTAGATAAGCa                                                                            <  2:615729/143‑1 (MQ=255)
                                                                            ccccACCCGGCGCTTATCCCCGCCCCCCTCAGTGGTTTAGCTTTGGAGGGGGCAATAAGAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTa                                                                      <  1:198327/149‑1 (MQ=255)
                                                                                    ggCGCTTATCCCGGGCACTCTCAGGGGCTTAGTTCTAGAGGGGGCGAAAAAAAAAATCTAAAAAACCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTc                                                              <  1:431937/149‑1 (MQ=255)
                                                                                         ttATCTCCGCCCCCCCCGGTGGTTACCTTTTGGAGGGGGCAAAAAAAAAAATCAAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCg                                                          <  2:774463/148‑1 (MQ=255)
                                                                                             cgccGGCTCTCTCAGTGGTTTAGCTCTTGAAGGGGCAATAAAAATATTCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTgg                                                   <  1:214225/149‑1 (MQ=255)
                                                                                             ccccGGCACTCTCAGTGGCTTAGCTTTTGAGGGGGCAATAAAAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCt                                                     <  1:734206/147‑1 (MQ=255)
                                                                                             ccccGCCACTCTCAGTGGTTTAGCTTTGAAGGGGGCAAAAAAAATAATCTCAAAAACTAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGa                                                         <  2:520573/143‑1 (MQ=255)
                                                                                               ccGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGAc                                                 >  2:13847/1‑151 (MQ=255)
                                                                                                      tctcGGTGGTTTAGCTGGTGAAGGGGAAAAAAAAATAATCTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAg                                             <  2:503847/148‑1 (MQ=255)
                                                                                                           gTGGTTTAGCTCTTGAAGGGTCAAAAAAAATAATCTAATAAACCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATcc                                     <  2:257994/151‑1 (MQ=255)
                                                                                                             ggCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATcc                                     >  2:405694/1‑149 (MQ=255)
                                                                                                                ttAGCTCTTGAGGGGGCAAAAAGAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGt                                 <  1:483134/150‑1 (MQ=255)
                                                                                                                  aGCTCTTGAGGGGGAGAAAAGAATACACTCATAACACCATCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGa                              <  1:594182/151‑1 (MQ=255)
                                                                                                                  aGCTCTTGAAGGGGCAAAAAAAATAATCAAATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCg                               <  1:689160/150‑1 (MQ=255)
                                                                                                                         gAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCg                       >  2:233220/1‑151 (MQ=255)
                                                                                                                             ggggCGAAAAAAATAATCTAAAAAAGATACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGaa                     <  2:66666/150‑1 (MQ=255)
                                                                                                                               ggcaaaaaGAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACtt                  <  1:258614/145‑1 (MQ=255)
                                                                                                                               ggCAATAAAAATAATCTAAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACtt                  <  1:142545/151‑1 (MQ=255)
                                                                                                                                     aaGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGc             >  2:803574/1‑150 (MQ=255)
                                                                                                                                     aaGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGc             >  2:185457/1‑150 (MQ=255)
                                                                                                                                               tcATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAgg  >  1:420307/1‑151 (MQ=255)
                                                                                                                                                cATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAgg  >  2:233698/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGG  >  NC_000913/1764640‑1764932

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