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,2500CG98.0% 163.5 / NA 51intergenic (‑349/+37)aroC/USA300HOU_1405chorismate synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (23/27);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCATATAATCATTAATTAAAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTC‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATAAAA  >  CP000730/1539119‑1539374
                                                                                                                                   |                                                                                                                             
atcatATAATCATTAATTAAAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTaaaa                                                                                                                     >  2:552773/1‑141 (MQ=255)
 tcatATAATCATTAATTAAAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTaaaaa                                                                                                                    <  1:791508/141‑1 (MQ=255)
 acatATAACCATTAACTAAAAATGATTAGCCGTTTTAACATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTaaaaa                                                                                                                    <  1:357192/140‑1 (MQ=255)
            ttaattaaAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTATTAAGTG‑CGTTTAAAAAAATTATGatat                                                                                                         >  2:429723/1‑141 (MQ=255)
             taattaaAAATGATTAGCAGTTTTTTGATAATCGTGTCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATAtt                                                                                                        >  1:987684/1‑141 (MQ=255)
              aagtaaAAATGATTAGTCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAATAACGAACTGTTCATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTa                                                                                                       <  1:33673/138‑1 (MQ=255)
                 taaAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTaaaa                                                                                                    >  1:770356/1‑141 (MQ=255)
                 taaAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTaaaa                                                                                                    >  2:43836/1‑141 (MQ=255)
                  aaaaaTGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAat                                                                                                   >  1:483635/1‑141 (MQ=255)
                      aTGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCCAATCCCGCAGTATTCACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATAc                                                                                               <  1:718944/141‑1 (MQ=255)
                      aTGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAAtataa                                                                                               >  2:950460/1‑140 (MQ=255)
                      aTGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATAc                                                                                               <  1:299839/141‑1 (MQ=255)
                      aTGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAAATGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATAc                                                                                               <  2:132079/141‑1 (MQ=255)
                          ttAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACatat                                                                                           <  2:662533/141‑1 (MQ=255)
                                gTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCACTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAg                                                                                     >  2:281922/1‑141 (MQ=255)
                                   ttattaTAATTGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAAGCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAgcgc                                                                                  <  2:947326/141‑1 (MQ=255)
                                     aaTATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATCCATATGTTAAGGGCtt                                                                                >  1:301817/1‑141 (MQ=255)
                                       tataATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACGAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTACGATATTAAAATATACATATGTTAAGCGCtttt                                                                              >  1:665209/1‑141 (MQ=255)
                                       tataATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCtttt                                                                              >  2:116958/1‑141 (MQ=255)
                                       tataATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATAAATATGTTAAGCGCtttt                                                                              >  2:466565/1‑141 (MQ=255)
                                       tataATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CATTTAAAAAAATTATGTTATTAAAATATACATATGTTAAGCGCtttt                                                                              >  2:157910/1‑141 (MQ=255)
                                       tataATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATAAAAAAGAACGGACTGTTAATTTTTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCtttt                                                                              >  2:288/1‑141 (MQ=255)
                                        ataATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGCTATTAAAATATACATATGTTAAGCGCTTTTg                                                                             <  1:368679/141‑1 (MQ=255)
                                        ataATCGTATCGCAACTTTAGCAAATAATATACTGACTACAGTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTg                                                                             <  2:141226/141‑1 (MQ=255)
                                           aTCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGCCACGGCT‑AAAAAAGTATGATGTTAACATAGACATATGTTAAGCGCTttttttt                                                                          >  1:225498/1‑140 (MQ=255)
                                              gTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTtatat                                                                       >  2:825485/1‑140 (MQ=255)
                                                  cGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTg                                                                   >  2:24066/1‑141 (MQ=255)
                                                    cAACTTTAGCAAATGATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGAc                                                                 >  2:674103/1‑141 (MQ=255)
                                                      acgttAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACtt                                                               <  2:947184/138‑1 (MQ=255)
                                                      aCTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACGGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACtt                                                               <  2:708340/141‑1 (MQ=255)
                                                        tttAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTACAATATACATATTTTAAGCGCTTTTGTTATACTTTGACTTCa                                                             <  1:885415/141‑1 (MQ=255)
                                                                     ttacttactACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCaa                                                <  2:887807/140‑1 (MQ=255)
                                                                      taCTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAAt                                               <  2:303759/141‑1 (MQ=255)
                                                                      taCTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTT‑AAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAAtt                                              >  2:717490/1‑141 (MQ=255)
                                                                        cTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATATAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTc                                             >  2:638035/1‑141 (MQ=255)
                                                                        cTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTc                                             <  1:619118/141‑1 (MQ=255)
                                                                            cTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCtttt                                         <  1:890309/141‑1 (MQ=255)
                                                                              aCATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTAc                                       <  2:483635/141‑1 (MQ=255)
                                                                                   tCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATa                                  >  1:529016/1‑141 (MQ=255)
                                                                                     ttAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACa                                <  1:393950/141‑1 (MQ=255)
                                                                                     ttAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACa                                <  1:281107/141‑1 (MQ=255)
                                                                                     ttAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACa                                <  1:786115/141‑1 (MQ=255)
                                                                                           cAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATaaaaa                          >  2:668677/1‑141 (MQ=255)
                                                                                            aaaaGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATaaaaaa                         <  1:552773/141‑1 (MQ=255)
                                                                                              aaGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAata                       <  1:616755/141‑1 (MQ=255)
                                                                                              aaGAACGAACTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACAAATGTTAAGCGATCTTTGTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAata                       <  1:982226/141‑1 (MQ=255)
                                                                                                    gaacTGTTAATTATTAATAAAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATg                 >  1:721031/1‑141 (MQ=255)
                                                                                                    gaacTGTTAATTATTAATAAAGTTAATCAATG‑CGGGTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATg                 >  1:986639/1‑141 (MQ=255)
                                                                                                     aacTGTTAATTATTAATAAAGTTAATCAGTG‑CGCTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATga                <  1:24066/141‑1 (MQ=255)
                                                                                                                   aataaAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATaaaa  <  1:429723/141‑1 (MQ=255)
                                                                                                                   aataaAGTTAATCAGTG‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAACATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATaaaa  <  1:643592/141‑1 (MQ=255)
                                                                                                                                   |                                                                                                                             
ATCATATAATCATTAATTAAAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTC‑CGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATAAAA  >  CP000730/1539119‑1539374

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