Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913448,5820AG23.3% 51.6 / 12.9 30F27L (TTC→CTC) ‡cyoCcytochrome o ubiquinol oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base A (13/10);  new base G (3/4);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCAC  >  NC_000913/448444‑448730
                                                                                                                                          |                                                                                                                                                    
tttCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCTTCGGATTTTggat                                                                                                                                            <  5:7792/149‑3 (MQ=255)
       cAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCTTCGGATTTTgga                                                                                                                                             >  1:126872/1‑140 (MQ=255)
       cAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCTTCGGATTTTgga                                                                                                                                             <  2:126872/141‑2 (MQ=255)
        aGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCTTCGGATTTTgg                                                                                                                                              <  4:94771/139‑1 (MQ=255)
        aGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCTTCGGATTTTgg                                                                                                                                              >  3:94771/1‑139 (MQ=255)
                      ttCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCTTCGGATTTTgga                                                                                                                                             >  2:115821/1‑125 (MQ=255)
                      ttCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCTTCGGATTTTgga                                                                                                                                             <  1:115821/126‑2 (MQ=255)
                         gAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGtt                                                                                                                   >  4:100010/1‑149 (MQ=255)
                              tGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTgc                                                                                                              <  1:52424/149‑1 (MQ=255)
                                tCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGaa                                                                                                                                                      >  2:174968/1‑107 (MQ=255)
                                       cTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGAtttt                                                                                                                                                <  6:72002/106‑1 (MQ=255)
                                       cTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGAtttt                                                                                                                                                >  5:72002/1‑106 (MQ=255)
                                                              ttctggtgAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCaaa                                                                                      >  3:107974/8‑141 (MQ=255)
                                                              ttctggtgAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCaaa                                                                                      <  4:107974/134‑1 (MQ=255)
                                                                             tAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCa                                                                                >  3:96465/1‑132 (MQ=255)
                                                                             tAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCa                                                                                <  4:96465/132‑1 (MQ=255)
                                                                                gTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGt                                                                              >  3:6013/1‑131 (MQ=255)
                                                                                gTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGt                                                                              <  4:6013/131‑1 (MQ=255)
                                                                                 tAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGc                                                           >  3:38717/1‑149 (MQ=255)
                                                                                    cAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCt                                                         >  4:4336/1‑148 (MQ=255)
                                                                                     aaaCAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATc                                                                                                                                 >  8:94982/1‑75 (MQ=255)
                                                                                     aaaCAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATc                                                                                                                                 <  7:94982/75‑1 (MQ=255)
                                                                                          aGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAAt                                                  >  6:20805/1‑149 (MQ=255)
                                                                                           gATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAAtc                                                 <  3:100010/149‑1 (MQ=255)
                                                                                                       aaTGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAg                                                                 <  3:49745/121‑1 (MQ=255)
                                                                                                       aaTGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGATCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAg                                                                 >  4:49745/1‑121 (MQ=255)
                                                                                                                        ggTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTc                    <  5:20805/149‑1 (MQ=255)
                                                                                                                           aGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGAt                 >  1:73034/1‑149 (MQ=255)
                                                                                                                            gATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGAtt                >  4:75038/1‑149 (MQ=255)
                                                                                                                                 aaaaTCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTa                                                       >  4:4458/1‑105 (MQ=255)
                                                                                                                                 aaaaTCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTa                                                       <  3:4458/105‑1 (MQ=255)
                                                                                                                                          aGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCAc  <  2:73034/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
TTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCAC  >  NC_000913/448444‑448730

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