Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,539,250 C→G intergenic (‑349/+37) aroC ← / ← USA300HOU_1405 chorismate synthase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,539,2500CG95.8% 68.6 / NA 24intergenic (‑349/+37)aroC/USA300HOU_1405chorismate synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (12/11);  minor base T (0/1);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCATCATATAATCATTAATTAAAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTCCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATAAAAATAATTTTTGT  >  CP000730/1539116‑1539385
                                                                                                                                      |                                                                                                                                       
tGCATCATATAATCATTAATTAAAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGAAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACAGTTAATTATTAATAAAGTTAATCAGTGCGTTTa                                                                                                                                   >  2:101978/1‑141 (MQ=255)
 acatcatATAATCATTAATTAAAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTaa                                                                                                                                  <  2:282247/140‑1 (MQ=255)
          aaTCATTAATTAAAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTAt                                                                                                                         >  1:97921/1‑141 (MQ=255)
           atcatttttaattaaTGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATg                                                                                                                        >  2:168399/14‑141 (MQ=255)
                       aaaTGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAatat                                                                                                            >  1:30358/1‑141 (MQ=255)
                                       tAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCt                                                                                            <  1:261002/141‑1 (MQ=255)
                                           ataATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTg                                                                                        <  2:244595/141‑1 (MQ=255)
                                             aaTCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGtt                                                                                      >  2:105676/1‑141 (MQ=255)
                                                  tATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACt                                                                                 >  2:102533/1‑141 (MQ=255)
                                                   aTCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACtt                                                                                >  1:341182/1‑141 (MQ=255)
                                                         aCTGTAGCAAGTTACATATTGATTACAAGTCTTAATATATATGATCATATTGTTAATTATTAATAAAGTTAATCAGTTCTTCTAAAAAAATTATTATATTAAAATATTCATACTTTAACCGTTTTTGTTTTACTTTGACtt                                                                          <  2:26317/141‑1 (MQ=255)
                                                               gCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCaaaaa                                                                    >  1:157558/1‑141 (MQ=255)
                                                                     aataTACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAAt                                                              <  2:259556/141‑1 (MQ=255)
                                                                          aCTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAAtt                                                         <  2:125937/141‑1 (MQ=255)
                                                                                              cAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCaa                                                                       <  2:57138/107‑1 (MQ=255)
                                                                                              cAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCaa                                                                       >  1:57138/1‑107 (MQ=255)
                                                                                                 aaGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAata                                  <  1:10657/141‑1 (MQ=255)
                                                                                                 aaGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAata                                  <  1:18022/141‑1 (MQ=255)
                                                                                                  aGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAatat                                 >  2:52995/1‑141 (MQ=255)
                                                                                                     acgaacTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATaaa                              <  2:297003/141‑1 (MQ=255)
                                                                                                          cTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATgaga                         >  2:196470/1‑141 (MQ=255)
                                                                                                                        taaAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATAAAaat           <  2:341182/141‑1 (MQ=255)
                                                                                                                           aGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATAAAaataat        <  2:143541/141‑1 (MQ=255)
                                                                                                                                 tCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATAAAAATAATTTTTGt  >  1:289748/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                       
TGCATCATATAATCATTAATTAAAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTCCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATAAAAATAATTTTTGT  >  CP000730/1539116‑1539385

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