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

GAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGG  >  NC_000913/1764640‑1764932
                                                                                                                                                |                                                                                                                                                      
gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGc                                                                                                                                                 <  2:461117/151‑1 (MQ=255)
  gaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTa                                                                                                                                               <  2:449456/151‑1 (MQ=255)
         aaaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑Cccgcc                                                                                                                                       <  1:648639/151‑1 (MQ=255)
           aaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGt                                                                                                                                     >  1:252252/1‑151 (MQ=255)
                 gCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAac                                                                                                                               >  1:702127/1‑151 (MQ=255)
                    gCGTACAAAAATAAAAAACCCCTTAAGTCTGTGCGACACGGGTTTAAGGGTTTCTCCCCCATCCGGCGCTTATCCCCGGCCCTCTCAGTGGTTTAGCTCTTGAG‑GGGGCGATAAGAATAATCTCATAAAGCTAC‑CCCGCCGTTTTAAcaca                                                                                                                            <  1:804224/151‑1 (MQ=255)
                      gTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACaa                                                                                                                           >  2:328215/1‑150 (MQ=255)
                          aaaaaTAAGAAACCCCTTAAGTCTGTGCGACACAGGTTAAAGGTTTTCTCCCCCACCCGGCGTTTATCTCCGGCCCTCCCAGTGGCTTAGCTTTTAAA‑GGGGCGATAAGAATAATCTAAAAAAGCTAC‑CCCGCCGTTTTAACACAAACTg                                                                                                                       <  2:213763/148‑1 (MQ=255)
                              aTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGAtt                                                                                                                  >  1:515778/1‑151 (MQ=255)
                                   aaacCCTTTAGGTTTGTGCGACACAGGCTTAGGGGTTTCTCCCCCATCCGGCCCTTATCTCCGGCCCTCTCAGTGGTTTAGCTTTTGAG‑GGGGCAATAAAAAAAATCTAAAAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGta                                                                                                              <  1:769120/150‑1 (MQ=255)
                                    aaaCCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGtatt                                                                                                            >  2:444818/1‑151 (MQ=255)
                                       ccttaaaGTTTGTGCGACACAGCCTAAGGGGTTTCTACCCCACCCGGCCTTTATCCCCGCCCCTCTCAGTGGCTGAGCTTTTAAA‑GGGGAAAAAAAAATAATCTAAAAAACCAAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTattatt                                                                                                         <  2:186995/146‑1 (MQ=255)
                                                       acacggaattAGGGGTTTCTCCCCCATCGGGCGTTTATCTCCGGCCCTCGCAGTGGCTTAGCTTTGGAA‑GGGGGAAAAAGAATAATCTAATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAg                                                                                          <  1:229951/142‑1 (MQ=255)
                                                          caGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGTGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCgg                                                                                      >  2:143797/1‑151 (MQ=255)
                                                            ggCTTAAGGTTTTCTCCCCCACCCGGCGCTTATCCCCGGCCCTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAAAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                                                                        <  1:81962/147‑1 (MQ=255)
                                                                  aGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATa                                                                                >  2:613835/1‑149 (MQ=255)
                                                                       ttCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGt                                                                          >  2:441376/1‑150 (MQ=255)
                                                                           cccccACCGGGCTCTTATCCCCGCCCCTCTCAGTGGTTTAGCTTTTGAG‑GGGGCGAAAAAAATAATCTAAAAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTa                                                                      <  1:253420/149‑1 (MQ=255)
                                                                            ccccACCCGTCACTAATCCCCGCCACCCTCAGTGGCTGAGCTTTAGAA‑GGGGCGTAAAAAATAATCAAAAAAAGTTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAACATCAGGCGGTAGATAAGCa                                                                            <  2:615729/143‑1 (MQ=255)
                                                                            ccccACCCGGCGCTTATCCCCGCCCCCCTCAGTGGTTTAGCTTTGGAG‑GGGGCAATAAGAATAATCTAAAAAAGCTAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTa                                                                      <  1:198327/149‑1 (MQ=255)
                                                                                    ggCGCTTATCCCGGGCACTCTCAGGGGCTTAGTTCTAGAG‑GGGGCGAAAAAAAAAATCTAAAAAACCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTc                                                              <  1:431937/149‑1 (MQ=255)
                                                                                         ttATCTCCGCCCCCCCCGGTGG‑TTACCTTTTGGAGGGGGCAAAAAAAAAAATCAAAAAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCg                                                          <  2:774463/148‑1 (MQ=255)
                                                                                             cgccGGCTCTCTCAGTGGTTTAGCTCTTGAA‑GGGGCAATAAAAATATTCTCAAAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTgg                                                   <  1:214225/149‑1 (MQ=255)
                                                                                             ccccGGCACTCTCAGTGGCTTAGCTTTTGAG‑GGGGCAATAAAAATAATCTAAAAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCt                                                     <  1:734206/147‑1 (MQ=255)
                                                                                             ccccGCCACTCTCAGTGGTTTAGCT‑TTGAAGGGGGCAAAAAAAATAATCTCA‑AAAACTAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGa                                                         <  2:520573/143‑1 (MQ=255)
                                                                                               ccGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGAc                                                 >  2:13847/1‑151 (MQ=255)
                                                                                                      tctcGGTGGTTTAGCTGGTGAA‑GGGGAAAAAAAAATAATCTAAAAAAGCAAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAg                                             <  2:503847/148‑1 (MQ=255)
                                                                                                           gTGGTTTAGCTCTTGAA‑GGGTCAAAAAAAATAATCTAATAAACCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATcc                                     <  2:257994/151‑1 (MQ=255)
                                                                                                             ggCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATcc                                     >  2:405694/1‑149 (MQ=255)
                                                                                                                ttAGCTCTTGAG‑GGGGCAAAAAGAATAATCTAAAAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGt                                 <  1:483134/150‑1 (MQ=255)
                                                                                                                  aGCTCTTGAG‑GGGGAGAAAAGAATACACTCATAACACCAT‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGa                              <  1:594182/151‑1 (MQ=255)
                                                                                                                  aGCTCTTGAA‑GGGGCAAAAAAAATAATCAAATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCg                               <  1:689160/150‑1 (MQ=255)
                                                                                                                         gAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCg                       >  2:233220/1‑151 (MQ=255)
                                                                                                                             ggggCGAAAAAAATAATCTAAAAAAGATAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGaa                     <  2:66666/150‑1 (MQ=255)
                                                                                                                               ggcaaaaaGAATAATCTAAAAAAGCAAC‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACtt                  <  1:258614/145‑1 (MQ=255)
                                                                                                                               ggCAATAAAAATAATCTAAAAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACtt                  <  1:142545/151‑1 (MQ=255)
                                                                                                                                     aaGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGc             >  2:803574/1‑150 (MQ=255)
                                                                                                                                     aaGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGc             >  2:185457/1‑150 (MQ=255)
                                                                                                                                               tcATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAgg  >  1:420307/1‑151 (MQ=255)
                                                                                                                                                cATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAgg  >  2:233698/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAA‑GGGGCGATAAGAATAATCTCATAAAGCTAA‑CCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGG  >  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: ‑