Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,613,762 T→C intergenic (‑31/‑37) scpA ← / → USA300HOU_1497 possible segregation and condensation protein A/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,613,7620TC96.5% 94.1 / NA 29intergenic (‑31/‑37)scpA/USA300HOU_1497possible segregation and condensation protein A/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (17/11);  total (18/11)
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

TCAAATTTTTGGATAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCATGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTAGAAGATGCATGGAAAGCTGAAAAT  >  CP000730/1613637‑1613897
                                                                                                                             |                                                                                                                                       
tCAAATTTTTGGATAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAAtt                                                                                                                   >  1:49225/1‑148 (MQ=255)
tCAAATTTTTGGATAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAAtt                                                                                                                   >  1:283981/1‑148 (MQ=255)
  aaaTTTTTGGATAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTgg                                                                                                                 >  1:292643/1‑148 (MQ=255)
           gATAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAAc                                                                                                        <  1:238875/148‑1 (MQ=255)
             tAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCATGAAAATAACTTACTTTTAAATTGGGGTGATAACaa                                                                                                      >  2:363535/1‑148 (MQ=255)
             tAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACaa                                                                                                      >  2:251894/1‑148 (MQ=255)
             tAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACaa                                                                                                      >  1:355315/1‑148 (MQ=255)
             tAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACaa                                                                                                      >  2:114724/1‑148 (MQ=255)
               aGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACaaaa                                                                                                    <  1:121959/148‑1 (MQ=255)
                  tGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGc                                                                                                 >  2:308667/1‑148 (MQ=255)
                  tGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGc                                                                                                 >  2:47320/1‑148 (MQ=255)
                        aaTAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAgg                                                                                           >  1:118909/1‑148 (MQ=255)
                               cTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGAt                                                                                    >  2:81600/1‑148 (MQ=255)
                                   tGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAAt                                                                                >  1:45385/1‑148 (MQ=255)
                                      tCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTc                                                                             <  1:139678/148‑1 (MQ=255)
                                       ccATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCt                                                                            <  1:309300/148‑1 (MQ=255)
                                         aTTGAGAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCtat                                                                          <  2:269435/148‑1 (MQ=255)
                                         aTTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCtat                                                                          <  1:304388/148‑1 (MQ=255)
                                         aTTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCtat                                                                          <  2:61846/148‑1 (MQ=255)
                                         aTTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCtat                                                                          <  2:264284/148‑1 (MQ=255)
                                             aaaGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTAAAGCACGAAAATAACTTACTTTTAAAATAAAAAAAAAAAAAAATGCAACAGGCACTGATAAATTTCtattata                                                                      <  2:347368/148‑2 (MQ=255)
                                                             cTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCa                                                      >  2:213215/1‑148 (MQ=255)
                                                                               ccTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCa                                    <  2:52880/148‑1 (MQ=255)
                                                                                 tCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATg                                  >  1:203514/1‑148 (MQ=255)
                                                                                          gTCATATGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTag                         >  2:33494/1‑148 (MQ=255)
                                                                                             atatGTGCATAATTCTCTAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTagaag                      >  1:287180/1‑148 (MQ=255)
                                                                                                           ctctAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTAGAAGATGCATGGAAAGCt        >  1:219658/1‑148 (MQ=255)
                                                                                                           ctctAAAATATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTAGAAGATGCATGGAAAGCt        >  1:234427/1‑148 (MQ=255)
                                                                                                                 aaTATTATAGCACGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTAGAAGATGCATGGAAAGCTGAAAAt  <  1:378681/148‑1 (MQ=255)
                                                                                                                             |                                                                                                                                       
TCAAATTTTTGGATAAGATGCAGCAATAAATCTAATGGTCCATTGAAAGCATCTAATTTAACTTCATACATAATATCTACCTCGTATTGCGTCATATGTGCATAATTCTCTAAAATATTATAGCATGAAAATAACTTACTTTTAAATTGGGGTGATAACAAAATGCAACAGGCACTGATAAATTTCTATTATCAATTTCATACTAAGCAACATTATTTCCTATGTCATGACATTTTAGAAGATGCATGGAAAGCTGAAAAT  >  CP000730/1613637‑1613897

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