Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,046,0680TA58.8% 9.1 / 31.2 34intergenic (‑167/+100)ygfF/gcvPputative oxidoreductase YgfF/glycine decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (3/11);  new base A (19/1);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 1.34e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGT  >  NC_000913/3045928‑3046208
                                                                                                                                            |                                                                                                                                            
ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCct                                                                                                                                              <  2:767732/141‑1 (MQ=255)
 tGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCctc                                                                                                                                             <  2:1387158/141‑1 (MQ=255)
 tGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCctc                                                                                                                                             <  2:901114/141‑1 (MQ=255)
 tGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCctc                                                                                                                                             <  2:659714/141‑1 (MQ=255)
 tGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCctc                                                                                                                                             <  1:788182/141‑1 (MQ=255)
          tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAg                                                                                                                                    <  1:841843/141‑1 (MQ=255)
                aTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACTCCCCAAAGGGGCgc                                                                                                                              >  2:1359671/1‑140 (MQ=255)
                 ttGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTACCCACACCCAAAACCACAACCCAAACGGGCGAg                                                                                                                             >  1:1445187/1‑141 (MQ=38)
                               gCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGc                                                                                                               <  2:1961128/141‑1 (MQ=255)
                                 tCACTCTTCATTGTTGGCACGCTTTTCCCCCTCGCCCATTCAGGGAGAGGGCCGGGGTGTGGGTAAATACTCGCGCCAGTGCCGTCCTGTTCCCCTCACCCTAACCCACTCCCCACAGGGGCGAGTTGACTGACTGAGCAc                                                                                                             <  2:334323/141‑1 (MQ=25)
                                 tCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCAc                                                                                                             <  2:334324/141‑1 (MQ=255)
                                                                              agagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCACACCCGAAACCACAACCCCAAGGGGCGAGGGGACTGTCGGAGCACTTTTTTACTTTTTAAACTGACTAAAAATGCGCCGAAGCGccggtt                                                                >  1:215090/1‑137 (MQ=255)
                                                                              agagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCCCCCAAAACCACACCCCAAACGGGCGCGGGGGCTGTGTGAGGCCATTTGTTATTTTTGATCTGACTAAAAATGCGCCCAAGCGCCGTTa                                                                >  1:215101/1‑141 (MQ=255)
                                                                               gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCGCACCCGAACCCACACCCCCAAGGGGCGAGGGGACAGTCTGAGCACTTTTGTTATTTGTCATCTGAAAAAAAAGGCGCCGAAGCGCCTTTat                                                               >  1:383656/1‑140 (MQ=255)
                                                                               gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCGCACCCAAAACCACACCCCAAAAGGGCGGGGGGGCTGTCTGAGGACGTTTTTAAGTTGTAGTATGACTAAAAAAGCGGCGAAGCGCCgtttt                                                               >  1:1283280/1‑139 (MQ=255)
                                                                               gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTATGTGCACTTTTGTACTTTGTCATATGACTAAAAACGCGCCGACGCGCCTTTAg                                                               >  2:473904/1‑141 (MQ=25)
                                                                               gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAAACCACTCCCCAAAGGGGCGAGGGGACTGTCTGTGCACTTTTGTTAGTTTTCATCTGACTAAAAAGGCGCCGAAGCGCCGTTAg                                                               >  2:147144/1‑141 (MQ=25)
                                                                               gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCAAAACCACACCCCCAAGGGGCGAGGGGACTGTGTGAGCCCGTTTGTAATTTTTGATCTGACTAAAAATGCGCCGGACCGCCGTTAg                                                               >  1:1290348/1‑141 (MQ=255)
                                                                                           gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCTCACCCAAACCCACACCCCAAAGGGGCGAGGGGTCTGTGTGAGGCCTTTTTTAATTTGTGATATGAATAAAAATGCGCCGAAGCGCCGTTATAAAATATTGGaa                                                   >  1:1471127/1‑141 (MQ=255)
                                                                                           gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCACACCCAAACCCACTCCCCCAAGGGGCGCGGGGTCTGTCTGAGCACGTTTGTAATTTTTAAACATACTAAAAATGCGCCGAAGCGCCTTTATAAAAAATTCGaa                                                   >  1:1922309/1‑141 (MQ=255)
                                                                                           gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTGATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGaa                                                   >  2:709044/1‑141 (MQ=255)
                                                                                           gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCTCTCCCCCAAGGGGGGAGGGGCCGGGCTGAGCACAGTTGTACTTTTTGATCTGACTAAAAAGGCGCCGAAGCGCCGTTTGAAAATATTCGaa                                                   >  2:709006/1‑141 (MQ=255)
                                                                                           gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCGAACCCACAACCCCAATGGGCGAGGGGACTGTCTGAGGACTTTTGTAATTTGTCATCTGAATAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGaa                                                   >  2:708964/1‑141 (MQ=255)
                                                                                           gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGAACCCAAACCCACACCCCAAAGGGGCGAGGGGACAGTCTGATCAAGTTTTTAATTTGTCATAATTCAAAAAATGCGCCGAAGCGCCTTTAGAAAAAACTCGaa                                                   >  2:673401/1‑141 (MQ=255)
                                                                                           gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCTAAACCACTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACGTTTGTACTTTGTCAACTGACTAAAAATGCGCCGAAGCGCCGTTAGTAAATATTCGaa                                                   >  2:193801/1‑141 (MQ=37)
                                                                                           gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTACCCTCACCCTAACCCACTACCCAAAGGGGCGAGGGGAATGTCTGAGCACATTTTTACTTTGTCATCAGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTAGaa                                                   >  2:50818/1‑141 (MQ=255)
                                                                                           gAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCTAACCCACTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACGTTTGTACTTTGTCATCTGACTAAAAATGCGCCGAAGCGCCTTTAGAAAATAGTCGaa                                                   >  2:673343/1‑141 (MQ=38)
                                                                                            aGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCACACCCCCCAGGGGCGGGGGGACGGGCGGAGCCCGTTTGTAAATTTTCATCTGACAAAAAATGCGCCCAACCGCCGTTAGAAAAAAGTCGaaa                                                  >  1:1404875/1‑140 (MQ=255)
                                                                                               gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAg                                               <  1:1609888/141‑1 (MQ=255)
                                                                                                tAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGt                                              <  1:110352/141‑1 (MQ=255)
                                                                                                tAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGt                                              <  1:111437/141‑1 (MQ=255)
                                                                                                       tCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCACTCCCCCAAGGGGCGAGGGGACTGTCTGAGCACTTTTTTACTTTTTGATATGACTAAAAATGCGCCGAAGCGCCTTTAGAAAATAGTCGAATAAGTGAATTaa                                       >  2:1396982/1‑140 (MQ=37)
                                                                                                        cgcaCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACACCCCCAAGGGGCGAGGGGACTGTCTGTGCACGTTTGTAATTTGTCAAATTAAAAAAAAGTCGCCGAAGCGCCGTTAGAAAATATTAGAATAATTGAATTACt                                      >  1:457857/1‑141 (MQ=255)
                                                                                                         gcgcCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCGTTAGAAAATAGTAGAATCAGTGAATTACTg                                     >  2:867725/1‑141 (MQ=255)
                                                                                                                      ccTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc                        <  1:1059375/141‑1 (MQ=255)
                                                                                                                      ccTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc                        <  1:1058574/141‑1 (MQ=255)
                                                                                                                                   cccTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGc           <  1:1359671/141‑1 (MQ=255)
                                                                                                                                       aaCCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGaa       <  1:95162/141‑1 (MQ=255)
                                                                                                                                       aaCCCTCTCCCCAAAGCGGCGAGGGGACTGCCCGCGCACTAGAGAACTATCTCAACCGACTAAAAAGGCGCCGACCCGCCATGAGACAATAGTCGAATCAGTGAATCACTGGTATTCGCTAATCGGTCCGCACCAGCAGaa       <  2:43583/141‑1 (MQ=255)
                                                                                                                                            tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACCGGt  >  2:319603/1‑141 (MQ=255)
                                                                                                                                            tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGt  >  2:210637/1‑141 (MQ=255)
                                                                                                                                            tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGt  >  2:1225601/1‑141 (MQ=255)
                                                                                                                                            tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGt  >  1:1548686/1‑141 (MQ=255)
                                                                                                                                            tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCGTTAGAAAATAGTCGAATCAGTGAATTACTGGTAGTCGCTAATCGGTACGCAGGAGCAGAACAGGt  >  2:319370/1‑141 (MQ=255)
                                                                                                                                            tctcCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTGTGTCATCTGACTAAAAAGGCGCCGAATCGCCTTAAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGt  >  2:319602/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
TTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGT  >  NC_000913/3045928‑3046208

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