Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,209,4200CG51.8% 11.1 / 36.7 28pseudogene (200/501 nt)stfEpseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (5/8);  new base G (14/0);  total (20/8)
Fisher's exact test for biased strand distribution p-value = 5.80e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: Frequency below/above cutoff threshold.

CCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGC  >  NC_000913/1209357‑1209487
                                                               |                                                                   
caaTGGGCATGAACCTCCGTTGAATACCCCTGAATACTTAAAATAGAGCGACCTGTATCCATCCcccg                                                                 >  2:353415/3‑68 (MQ=255)
  aTGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATcgcccgcc                                                               >  1:534089/1‑68 (MQ=255)
      gCATGATCCTCCGTTGCGTACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCGCCCGCCCGTc                                                           >  2:224776/1‑68 (MQ=255)
           aTCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATccc                                                      <  2:608009/68‑1 (MQ=255)
             cctccGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCGCCCGCCCGTCATCCCAg                                                    >  2:884299/1‑68 (MQ=255)
             cctccGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCGCCCGCCCGTCATCCCAg                                                    >  1:786780/1‑68 (MQ=255)
             cctccGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAg                                                    >  2:940108/1‑68 (MQ=255)
                  gTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACg                                               <  2:461478/68‑1 (MQ=255)
                  gTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACg                                               >  1:912970/1‑68 (MQ=255)
                   ttGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCGCCCGCCCGTCATCCCAGCCACGa                                              >  1:1059409/1‑68 (MQ=255)
                       aTACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCGCCCGCCCGTCATCCCGGCCACGAATaa                                          >  2:520474/1‑68 (MQ=255)
                            ccTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTca                                     <  1:41487/68‑1 (MQ=255)
                                aaTACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTcaccac                                 <  2:27921/68‑1 (MQ=255)
                                    cTTAAAATAGAGCGACCTGTATCAATCGCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTaa                             >  1:525621/1‑68 (MQ=255)
                                    cTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTaa                             <  2:291283/68‑1 (MQ=255)
                                          aTAGAGCGACCTGTATCGATCGCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAgg                       >  2:932134/1‑68 (MQ=255)
                                          aTAGAGCGACCTGTATCAATCGCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAgg                       >  2:448127/1‑68 (MQ=255)
                                          aTAGAGCGACCTGTATCAATCGCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAgg                       >  2:494597/1‑68 (MQ=255)
                                          aTAGAGCGACCTGTATCAATCGCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAgg                       >  2:750998/1‑68 (MQ=255)
                                           tAGAGCGACCTGTATCGATCCCCCGCCCGACATCCCAGCCACGAATAAACTCCCCACGTAAATCAGGc                      >  2:452881/1‑68 (MQ=255)
                                           tAGAGCGACCTGTATCAATCGCCCGCCCGTCATCCCAGCCGCGAATAAACTCACCGCGTGAATCAGGc                      >  2:397347/1‑68 (MQ=255)
                                           tAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTcaccac                                 >  2:323360/1‑57 (MQ=255)
                                           tAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTcaccac                                 <  1:323360/57‑1 (MQ=255)
                                                cGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAattt                 <  2:124648/68‑1 (MQ=255)
                                                 gACCTGTATCAATCGCCCGCCCGTCATCCCAGCCACGAATAAACTCACCGCGTAAATCAGGCTattta                >  1:777815/1‑68 (MQ=255)
                                                 gACCTGTATCAATCGCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAattta                >  1:163684/1‑68 (MQ=255)
                                                        aTCAATCGCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCg         >  1:274951/1‑68 (MQ=255)
                                                         tCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCgg        <  2:589194/68‑1 (MQ=255)
                                                         acAATCACCCGCACGTCATCCCAGCCACGAAAAAACTCACCACGTAAATCAGGCAATTTATTTGACgg        >  2:318787/2‑68 (MQ=255)
                                                               ccccGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAACTTATTTGTCGGATAAGc  <  2:131289/68‑1 (MQ=255)
                                                               |                                                                   
CCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGC  >  NC_000913/1209357‑1209487

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