Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi2,383,9560AT22.2% 38.7 / 14.0 27intergenic (+32/+39)elaD/yfbKprotease, capable of cleaving an AMC‑ubiquitin model substrate/Von Willebrand factor domain putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base T (2/4);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-02
Rejected as consensus: Frequency below/above cutoff threshold.

CAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAC  >  NC_000913_3_pae_pgi/2383810‑2384094
                                                                                                                                                  |                                                                                                                                          
cAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAat                                                                                                                                          <  1:178064/149‑1 (MQ=255)
cAGAAAGTGTTATGTAAATAAAAATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGGCAGCAAat                                                                                                                                          <  2:218302/149‑1 (MQ=255)
 aGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCACAAGCTGAAGCAGCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAata                                                                                                                                         >  1:528107/1‑149 (MQ=255)
  aaaaTTGTAAATTAGAAAAAGATGGCGATCAAAATAACAATAAGTATAATAGAAACTTATAAATAAAACACAAACAAACAAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAAATAAGTCGACCACTCAACAActat                                                                                                                                        <  2:605518/148‑1 (MQ=255)
        gTTATTTAGATAAATATGGCGATACAAATGGCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAAACCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGAc                                                                                                                                  <  2:298601/149‑1 (MQ=255)
         ttttttAGATAAAGAGTGAAATGCAAAAATAAATCATTAGTATAGAAACTTAGAAATAAAACACTCACAACCCAAGAATAGAGCATCCGGCAAAAGAGTGAGTGAACATCATTGTTAATTAATGCGACCAGTCAGCAAATATGCTGACt                                                                                                                                 <  2:357436/146‑1 (MQ=255)
         ttATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACt                                                                                                                                 <  2:389408/149‑1 (MQ=255)
          tATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACTg                                                                                                                                <  1:459363/149‑1 (MQ=255)
                 aTAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAAc                                                                                                                         >  1:574516/1‑149 (MQ=255)
                      tatGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTGAACATCAAAGTTAATGAATACGACCAGTCATCAAATATGCTGACTAGTGGAACATAgg                                                                                                                    >  1:150645/1‑148 (MQ=255)
                          gCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCCTCCGGCAAAAGAGTCAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCCAACATAGtttcc                                                                                                                >  2:405539/1‑147 (MQ=255)
                                        aatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTCACATCATTGTTAATTAATTCGCCCAGTCAGCAAATATGCTAACTGGTCGAACATAGTTTGCTAAACAGTCATCac                                                                                                  >  2:599409/1‑149 (MQ=255)
                                                          ttAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTcaca                                                                                >  2:391724/1‑149 (MQ=255)
                                                                     acacTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAg                                                                     >  2:376242/1‑149 (MQ=255)
                                                                                ccAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCCCACCATCCGCCAGTTCAATCAGGc                                                          >  2:205107/1‑149 (MQ=255)
                                                                                  aaaaaTAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCatatt                                                                                                                                         <  1:94062/68‑2 (MQ=255)
                                                                                                      aaGAGTGAGTCAACATCATTGTTAATTAATTCGACCAGTCAGAAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAAca                                    <  2:232843/149‑2 (MQ=255)
                                                                                                                  aCATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTc                        <  1:341501/149‑1 (MQ=255)
                                                                                                                  aCATCATTGTAAAGTAAGTCGACCAGACAGCAAATATGCCGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTc                        <  1:150254/149‑1 (MQ=255)
                                                                                                                       atTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagtca                                                                                  <  1:265434/86‑1 (MQ=255)
                                                                                                                       atTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagtca                                                                                  >  2:265434/1‑86 (MQ=255)
                                                                                                                           ttaattaaTTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACGGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGc               >  2:365070/1‑149 (MQ=255)
                                                                                                                                   ttCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCCTCACTGACTGATGTCAGTCACACCAGCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGCATATTCACAATTTGCCTGCTCAg       >  2:293993/1‑149 (MQ=255)
                                                                                                                                    tCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCACCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGc      <  1:599409/149‑1 (MQ=255)
                                                                                                                                        ccAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAc  >  2:271962/1‑149 (MQ=255)
                                                                                                                                        ccAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAc  <  1:365079/149‑1 (MQ=255)
                                                                                                                                          aGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGtcatca                                                                                                   >  2:542751/1‑50 (MQ=21)
                                                                                                                                          aGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGtcatca                                                                                                   <  1:542751/50‑1 (MQ=21)
                                                                                                                                                  |                                                                                                                                          
CAGAAATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAC  >  NC_000913_3_pae_pgi/2383810‑2384094

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