Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*U000962,947,4650T.23.1% 66.7 / 16.7 27intergenic (+2/‑32)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base T (7/13);  new base . (2/4);  total (9/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAT  >  U00096/2947318‑2947575
                                                                                                                                                    |                                                                                                              
ctaaTGCGTGATCGGGTTATACAAAAGGGTGCATGAATACTCGAGCGGAGTATTGTCCGCTTACTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAtc                                                                                                               <  2:964106/147‑2 (MQ=17)
 taaTGCGTGATCGGGTTAAAATTAAGGCTGAATGAAATCTCGAGCGGAGTATAGTGCGCTTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAtca                                                                                                              <  2:964126/148‑3 (MQ=18)
  aaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaa                                                                                                             <  2:964107/149‑1 (MQ=255)
   aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTatt                                                                                                            >  2:964165/1‑147 (MQ=255)
    tCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaat                                                                                                           >  2:964152/1‑148 (MQ=255)
          aTCGGGGTAAAATAAAGCTTGCATGAAAACTCTATCGGAGTATATTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaatt                                                                                                          <  2:964108/144‑3 (MQ=18)
          aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgatg                                                                                                    >  2:964151/1‑149 (MQ=255)
           tCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTTCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                         >  1:964108/1‑144 (MQ=255)
           tCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAt                                                                                                                >  1:964114/1‑137 (MQ=255)
               ggTAAAAAACAAGGTTGCAAGAAAACGCGATCGTAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGt                                                                                                <  2:964109/149‑1 (MQ=17)
                ttAAAAACACGGTTCCATGAAAACCCGAGCGGAGTATAGTGCGAATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaa                                                                                              <  2:964123/148‑1 (MQ=17)
                ttAAAAACACGCTTGCATGAAAAATCGTGCGTAGTATTGTGCGCTTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaa                                                                                              <  2:964144/148‑1 (MQ=14)
                gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTa                                                                                               <  1:964111/149‑1 (MQ=255)
                gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAACCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTa                                                                                               <  1:964110/149‑1 (MQ=255)
                 tAAAATTAAGGTTGCATGAAAACGCGAGCGTAGTATAGTGCGCCTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaa                                                                                             <  2:964114/149‑1 (MQ=17)
                 tAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaa                                                                                              <  1:964113/149‑1 (MQ=255)
                     aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaagca                                                                                         >  2:964141/1‑145 (MQ=255)
                      aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAg                                                                                         <  1:964116/149‑1 (MQ=255)
                       aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGt                                                                                        <  2:964117/149‑1 (MQ=255)
                       aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGt                                                                                        <  2:964118/149‑1 (MQ=255)
                         aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTa                                                                                       <  1:964119/149‑1 (MQ=255)
                         aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTa                                                                                       <  1:964121/149‑1 (MQ=255)
                           ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAATTTTGAGGAGGTAAAGCAGTACg                                                                                     >  1:964143/1‑149 (MQ=255)
                            gTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAATTTTTATAAGTTAAAGCAGTAcag                                                                                    >  1:964148/1‑147 (MQ=18)
                               gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGa                                                                                 <  2:964122/149‑1 (MQ=255)
                                  tGAAAACGCGAGCGGAGTATAGTGAGCATCCACGGACGCGGGGTGTAGCATCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCTGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGActc                                                                              <  2:964125/149‑3 (MQ=255)
                                        cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAt                                                                                                                >  1:964149/1‑109 (MQ=255)
                                             gCGGAGTATAGTGCGCTTCCACGGACGAGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCGCGCAACCAAATAAAATTTGATGAAGTAATGCAGTACGGTGACGCGGGGTGGagca                                                                   <  2:964128/149‑1 (MQ=17)
                                                                                                                                                ccAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAt  >  1:964185/1‑115 (MQ=255)
                                                                                                                                                 cAATTAAAATTTGTTGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtcgt                              >  2:964135/1‑86 (MQ=255)
                                                                                                                                                 cAATTAAAATTTGTTGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtcgt                              <  1:964135/86‑1 (MQ=255)
                                                                                                                                                    |                                                                                                              
CAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAT  >  U00096/2947318‑2947575

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