Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,413,571 T→C intergenic (‑106/‑293) acrS ← / → acrE acrAB operon transcriptional repressor/cytoplasmic membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,413,5710TC96.0% 84.4 / NA 25intergenic (‑106/‑293)acrS/acrEacrAB operon transcriptional repressor/cytoplasmic membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (10/14);  minor base A (1/0);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATA  >  NC_000913/3413508‑3413632
                                                               |                                                             
aaaaaaGCGAAGGGTAATCTATCAACTAATGTGTATATATAAGAGAGGCTAATATTCAAATGCCATaa                                                           >  2:182490/1‑68 (MQ=255)
   aTAGCGAAGGTTAATCTATCACCTAATGCGTATTTATACGAGAGGCTAATATTGAGTCGCCATAAATc                                                        <  2:1138763/68‑1 (MQ=255)
      gCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCCATAAATCGtt                                                     <  1:663118/68‑1 (MQ=255)
       cGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCCATAAATCGTta                                                    <  2:274782/68‑1 (MQ=255)
          aGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCCATAAATCGTtaaat                                                 <  2:933849/68‑1 (MQ=255)
            tttAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCCATAAATCGTtaaatta                                               <  1:91053/67‑3 (MQ=255)
            gTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTATTATTGAGTTGCCATAAATCGTtaaataa                                               <  2:700510/68‑1 (MQ=255)
                  cTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCCATAAATCGTtaaataa                                               <  1:327561/62‑1 (MQ=255)
                  cTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCCATAAATCGTtaaataa                                               >  2:327561/1‑62 (MQ=255)
                     tCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCCATAAATCGTTAAATAAATAatatat                                      <  2:198684/68‑1 (MQ=255)
                       aCCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCCATAAATCGTTAAATAAATAatatatat                                    >  2:100391/1‑68 (MQ=255)
                        ccTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCCATAAATCGTTAAATAAATAATATAtatt                                   <  2:615615/68‑1 (MQ=255)
                            aTGTGTATTTATACGAGAGGCTAATATTGAGTTGCCATAAATCGTTAAATAAATAATATATATTAttt                               <  2:1115040/68‑1 (MQ=255)
                              gtgtATTTATACGAGAGGCTAATATTGAGTTGCAATAAATCGATAAATAAATAATATATATTAATTAc                             >  2:351840/1‑68 (MQ=255)
                                gtATTTATACGAGAGGCTAAAATTGAGTTGCCATAAATCGTTAAATAAATAATATATATTATTTACCt                           <  2:125299/68‑1 (MQ=255)
                                  aTTTATACGAGAGGCTAATATTGAGTTGCCATAAATCGTTAaataaata                                            >  2:423749/1‑49 (MQ=255)
                                  aTTTATACGAGAGGCTAATATTGAGTTGCCATAAATCGTTAaataaata                                            <  1:423749/49‑1 (MQ=255)
                                           agagGCTAATATTGAGTTGCCATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTc                >  1:778116/1‑68 (MQ=255)
                                                   atatTGAGTTGCCATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATc        <  1:215698/68‑1 (MQ=255)
                                                     atTGAGATGCCATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCaa      >  1:38522/1‑68 (MQ=255)
                                                      ttGAGTTGCCATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAat     >  1:272562/1‑68 (MQ=255)
                                                         aGTTGCCATAAGTCGTTAAATAAATAATATATATTATATACCTAAGATACATTCACTACATCAatata  >  2:1126126/1‑68 (MQ=255)
                                                         aGTTGCCATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAatata  >  2:1126127/1‑68 (MQ=255)
                                                         aGTTGCCATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAatat   >  2:836081/1‑67 (MQ=255)
                                                         aGTTGCCATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAatat   <  1:836084/67‑1 (MQ=255)
                                                               |                                                             
AAAATAGCGAAGGTTAATCTATCACCTAATGTGTATTTATACGAGAGGCTAATATTGAGTTGCTATAAATCGTTAAATAAATAATATATATTATTTACCTAAGATACATTCACTACATCAATATA  >  NC_000913/3413508‑3413632

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