Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,0140AC29.2% 45.2 / 14.0 24intergenic (‑32/+571)ompF/asnSouter membrane porin F/asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (17/0);  new base C (0/7);  total (17/7)
Fisher's exact test for biased strand distribution p-value = 2.89e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGT  >  NC_000913/986950‑987078
                                                                |                                                                
gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca                                                                  >  1:3035639/1‑65 (MQ=255)
gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca                                                                  >  1:996640/1‑65 (MQ=255)
gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca                                                                  >  1:3979444/1‑65 (MQ=255)
gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca                                                                  >  1:3347543/1‑65 (MQ=255)
gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca                                                                  >  1:215795/1‑65 (MQ=255)
gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca                                                                  >  1:2100459/1‑65 (MQ=255)
gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca                                                                  >  1:1970384/1‑65 (MQ=255)
gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca                                                                  >  1:1801520/1‑65 (MQ=255)
gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca                                                                  >  1:1610156/1‑65 (MQ=255)
gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca                                                                  >  1:1170050/1‑65 (MQ=255)
gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca                                                                  >  1:1324085/1‑65 (MQ=255)
 aCGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac                                                                 <  1:3217454/65‑1 (MQ=255)
 aCGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac                                                                 <  1:1893149/65‑1 (MQ=255)
 aCGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac                                                                 <  1:337775/65‑1 (MQ=255)
     tCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc                                                             >  1:1288771/1‑65 (MQ=255)
     tCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc                                                             >  1:344334/1‑65 (MQ=255)
     tCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc                                                             >  1:484222/1‑65 (MQ=255)
       aCTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACACCTGCCa                                                           >  1:2292161/1‑65 (MQ=255)
       aCTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCa                                                           >  1:2474441/1‑65 (MQ=255)
       aCTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCa                                                           >  1:3411354/1‑65 (MQ=255)
           ccAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGc                                                       >  1:1372917/1‑65 (MQ=255)
             aaaaTTTTGCCCTTCATCATTTTTTTTTCCCCTCTTTTTTTTTTTTATCCCCCCTGCCACTGCCg                                                     <  1:2214021/63‑1 (MQ=255)
               aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTc                                                   >  1:3524726/1‑65 (MQ=255)
               aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTc                                                   >  1:865334/1‑65 (MQ=255)
                     gcgcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAg                                             >  1:349850/1‑65 (MQ=255)
                      cgcTTCATTTTTATTTATTACCCTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGt                                            <  1:34891/65‑1 (MQ=255)
                      cgcTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATCCCCCCTGCCACTGCCGTCAATAAGt                                            <  1:1676552/65‑1 (MQ=255)
                                    ttattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaaa                              >  1:836740/1‑65 (MQ=255)
                                         aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTAc                         >  1:3374614/1‑65 (MQ=255)
                                          cccTCATGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg                        <  1:2784433/65‑1 (MQ=255)
                                          cccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg                        >  1:3737350/1‑65 (MQ=255)
                                          cccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg                        >  1:2277365/1‑65 (MQ=255)
                                           ccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg                       >  1:1108057/1‑65 (MQ=255)
                                           ccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg                       >  1:2216484/1‑65 (MQ=255)
                                           ccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg                       >  1:3614178/1‑65 (MQ=255)
                                             tCATGGTTTTTTTTAGGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa                     <  1:3850968/65‑1 (MQ=255)
                                                 ttttttttttATCACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt                 <  1:2716303/63‑1 (MQ=255)
                                                 tgTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt                 <  1:3999771/64‑1 (MQ=255)
                                                 gtttttttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt                 <  1:3712255/63‑1 (MQ=255)
                                                 gtttttttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt                 <  1:1891469/63‑1 (MQ=255)
                                                 gtttttttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt                 <  1:111989/63‑1 (MQ=255)
                                                 gtttttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt                 <  1:3234100/63‑1 (MQ=255)
                                                 gtttttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt                 <  1:146488/63‑1 (MQ=255)
                                                 gtttttttttAGGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt                 <  1:3941153/63‑1 (MQ=255)
                                                 ggTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTGACGGAACTAt                 <  1:3027592/65‑1 (MQ=255)
                                                   ttttttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg               <  1:393357/65‑1 (MQ=255)
                                                     ttttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgat             >  1:875881/1‑65 (MQ=255)
                                                         ttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga         <  1:2533364/65‑1 (MQ=255)
                                                         ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga         >  1:995196/1‑65 (MQ=255)
                                                          tATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag        <  1:3496678/65‑1 (MQ=255)
                                                                acCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgt  >  1:2070541/1‑65 (MQ=255)
                                                                |                                                                
GACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGT  >  NC_000913/986950‑987078

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