Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913225,4670GT21.6% 50.4 / 16.1 37intergenic (+10/‑33)ileV/alaVtRNA‑Ile/tRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base G (17/12);  new base T (8/0);  total (25/12)
Fisher's exact test for biased strand distribution p-value = 3.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTACCT  >  NC_000913/225438‑225616
                             |                                                                                                                                                     
aaGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTgg                               >  1:969095/1‑150 (MQ=11)
aaGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTgg                               <  2:1550129/150‑1 (MQ=11)
 aGTCCACTCAGGCCTACCAAATTTGCACGGCACATTTGAAGAGGTTTTAACTACATGGTATGGGGCTAAAGCTCAGCTGGGAGCGCGCCTGCTTTGCACGCAGGCGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCGCTGTAGTGGt                              >  1:354824/1‑150 (MQ=11)
 aGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGt                              <  2:1892314/150‑1 (MQ=31)
 aGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGt                              >  2:1708777/1‑150 (MQ=31)
 aGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGGTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGt                              >  2:1617541/1‑150 (MQ=14)
  gTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGTGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGtt                             <  1:1343879/150‑1 (MQ=17)
   tCCACTCAGGCCTAACAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTta                            <  2:1889073/150‑1 (MQ=17)
    ccACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaa                           >  1:1931185/1‑150 (MQ=32)
      aCTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGTCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaaat                         <  2:1391284/150‑1 (MQ=17)
      aCTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaaat                         <  2:955395/150‑1 (MQ=32)
      aCTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTtaaat                         <  1:1560111/150‑1 (MQ=32)
           ggCCTTCGAAATTTGCAAGGCAAATTTGAAGATGTTTTAACTACCTGCTTTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGAAGTGGTTAAATaaaaa                    <  2:402946/150‑1 (MQ=11)
           ggCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATaaaaa                    >  1:435404/1‑150 (MQ=32)
            gCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAATTaaaaaa                   >  1:1905343/1‑150 (MQ=17)
            gCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATaaaaaa                   >  1:2082082/1‑150 (MQ=32)
            cctcttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGCAGGTCTGCCGATCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATaaaaaa                   >  2:1053427/7‑150 (MQ=11)
             ccTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAAt                  >  1:563890/1‑150 (MQ=32)
               cttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATAc                >  2:626525/4‑150 (MQ=255)
               cttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATAc                >  2:68590/4‑150 (MQ=255)
                ttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACt               >  2:1016178/3‑150 (MQ=12)
                ttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACt               >  1:807775/3‑150 (MQ=12)
                ttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACt               >  1:925155/3‑150 (MQ=12)
                ttcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACt               >  1:739179/3‑150 (MQ=12)
                 tcAAATTTGCCGTGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACtt              >  1:2139405/2‑150 (MQ=12)
                 ccAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACtt              >  1:1853732/1‑150 (MQ=32)
                  cAAATTTGCACGGCAATTTTGAAGAGTTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCATGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc             <  1:848544/150‑1 (MQ=12)
                  cAAATTTGCACGGCAAATTTGAAGTGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATTGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc             >  1:1372237/1‑150 (MQ=17)
                  cAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGATGTCTTCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc             >  2:1364554/1‑150 (MQ=17)
                  cAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTc             >  1:949496/1‑150 (MQ=32)
                   aaaTTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCa            >  1:476446/1‑150 (MQ=32)
                      tttGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCagag         <  1:1289256/150‑1 (MQ=32)
                       ttGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAgt        >  2:570610/1‑150 (MQ=32)
                        tGCCCGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACTCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATCCTTCAGAgtg       >  1:142712/1‑150 (MQ=12)
                          cACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCAGCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTa     <  1:1195106/150‑1 (MQ=17)
                          cACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCAGGGAGAGCGCCTGCTTTGCACGCATGATGTCTTCGGTTCGATCCCGCATTGCTACACCAGCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTa     >  1:1349534/1‑150 (MQ=11)
                           aCGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTAc    <  1:1898822/150‑1 (MQ=32)
                             ggCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTACCt  >  1:278939/1‑150 (MQ=32)
                             ggCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCAACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTAcca  >  2:1039589/1‑149 (MQ=17)
                             |                                                                                                                                                     
AAGTCCACTCAGGCCTACCAAATTTGCACGGCAAATTTGAAGAGGTTTTAACTACATGTTATGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCTTTGCACGCAGGAGGTCTGCGGTTCGATCCCGCATAGCTCCACCATCTCTGTAGTGGTTAAATAAAAAATACTTCAGAGTGTACCT  >  NC_000913/225438‑225616

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