Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007779_cfb1,097,5800CT24.1% 35.0 / 11.4 29intergenic (+329/+226)ycdU/serXinner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base C (7/15);  new base T (0/7);  total (7/22)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAC  >  NC_007779_cfb/1097334‑1097829
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         
aGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACg                                                                                                                                                                                                                                                        <  1:162882/250‑1 (MQ=255)
        tgCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGCTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGaaaaaaaa                                                                                                                                                                                                                                                <  1:1125535/250‑1 (MQ=255)
               tctCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTa                                                                                                                                                                                                                                         <  1:1247925/248‑1 (MQ=255)
                cttccgatctccTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTAc                                                                                                                                                                                                                                        <  1:29661/240‑1 (MQ=255)
                          ccTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGAATCGCTTCCCTCCCCCCGCGGGCCGCCCCCCCCCTGCGTTCTCGGTGTGTCCCCCCGGCTGTCGCCAGCTTTCCACCCCCCGCCGGGGCGTCTCCTCCCCCCCCCGCGTGCGAATTACACAAAAAAAAAcccctcacatatcgaa                                                                                                                                                                                                                                >  2:29661/1‑238 (MQ=9)
                           cTAACAGACGTGGGCGGGGGCCCCCTGTGTGTACAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCATCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAg                                                                                                                                                                                                                             <  1:1152499/250‑1 (MQ=255)
                                                                                  tACTGTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCg                                                                                                                                                                      <  1:58683/250‑1 (MQ=255)
                                                                                   aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCgg                                                                                                                                                                     <  1:976946/250‑1 (MQ=255)
                                                                                                 tcttctTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCTCTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGTTTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt                                                                                                                                                       <  1:1294305/250‑1 (MQ=255)
                                                                                                 tcttctTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt                                                                                                                                                       <  1:538526/250‑1 (MQ=255)
                                                                                                 gcttcggtAAAGAGGGCTGAGAGGGGAGGATTGACTTGTTCCGCACGCCCTGCGGGCATTTCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTATCGACCCCAGGGCGGGTCTTCTCATGCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACt                                                                                                                                                       <  2:1524811/243‑1 (MQ=17)
                                                                                                       ttAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTCGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCCTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTc                                                                                                                                                 <  1:1454395/250‑1 (MQ=255)
                                                                                                            tatGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGt                                                                                                                                            <  1:1552534/250‑1 (MQ=255)
                                                                                                                  ggTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTcc                                                                                                                                      <  1:870667/250‑1 (MQ=255)
                                                                                                                                   cTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGt                                                                                                                     <  1:379918/250‑1 (MQ=255)
                                                                                                                                          cgctcgcCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGGGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGaa                                                                                                              <  1:996623/250‑1 (MQ=255)
                                                                                                                                               gcCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATgg                                                                                                                                                                                                            >  2:646301/1‑151 (MQ=33)
                                                                                                                                               gcCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATgg                                                                                                                                                                                                            <  1:646301/151‑1 (MQ=33)
                                                                                                                                                cccTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGc                                                                                                                                                                                                                            >  2:961133/1‑134 (MQ=34)
                                                                                                                                                cccTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGc                                                                                                                                                                                                                            <  1:961133/134‑1 (MQ=34)
                                                                                                                                                               gCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATc                                                                                         <  1:182666/250‑1 (MQ=255)
                                                                                                                                                                                      cAACTGGCTGTCGCCATTTGTCGAACCCAGGTCGGGTCTTCTCATCCCCCCGGTGTGTGC‑AA‑TATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGa                                                                  <  1:486315/250‑1 (MQ=255)
                                                                                                                                                                                      cAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATAAGGCGGTGAgg                                                                                                                                                                                                    >  2:445626/1‑120 (MQ=25)
                                                                                                                                                                                       aaCTGGCTCTCGCCAGTTGTCGACCCCCGGTCGGGGCTTGTCATCCCCCCGGTTTTTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCGGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGaa                                                                 <  1:1433132/250‑1 (MQ=255)
                                                                                                                                                                                               gTCGCCAGTTGTCGAACCCCGGTC‑GGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGcc                                                        <  1:410510/250‑1 (MQ=255)
                                                                                                                                                                                                                cccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGa                                                                                                                                                                                                                                                       >  2:1215494/1‑43 (MQ=38)
                                                                                                                                                                                                                      cGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGaaaaaa                                                                                                                                                                                                                                                  >  2:564026/1‑42 (MQ=38)
                                                                                                                                                                                                                                      cccGGTGTGTGC‑AA‑TATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCGTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTT‑AATATGGCGGTGa                 <  1:1500312/250‑1 (MQ=14)
                                                                                                                                                                                                                                            gtgtgC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCttt                                                                                                                                                                                                                    >  2:974138/1‑50 (MQ=38)
                                                                                                                                                                                                                                             tgtgC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGctct                                                                                                                                                                                                                         >  2:1640644/1‑44 (MQ=38)
                                                                                                                                                                                                                                                      tACACGAAAAAAAAGCCGGTACTTTCGTACGGGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGAttgac    <  1:1744244/250‑4 (MQ=17)
                                                                                                                                                                                                                                                        cacGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGAttgactc  <  1:56295/250‑6 (MQ=11)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         
AGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGC‑AA‑TACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAC  >  NC_007779_cfb/1097334‑1097829

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