Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913448,5790CT22.6% 56.8 / 13.4 31G28R (GGA→AGA) cyoCcytochrome o ubiquinol oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base C (14/10);  new base T (3/4);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATT  >  NC_000913/448444‑448716
                                                                                                                                       |                                                                                                                                         
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)
                                                                                                                                       |                                                                                                                                         
TTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATT  >  NC_000913/448444‑448716

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: ‑