Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,433,718 T→C intergenic (‑44/+273) pinR ← / ← ynaE Rac prophage; putative site‑specific recombinase/cold shock protein, Rac prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,433,7180TC87.2% 76.0 / 7.3 47intergenic (‑44/+273)pinR/ynaERac prophage; putative site‑specific recombinase/cold shock protein, Rac prophage
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (20/21);  total (23/24)
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

TCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTC  >  NC_000913/1433622‑1433766
                                                                                                |                                                
taggtggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAg                                               >  2:1316859/3‑100 (MQ=16)
tCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAAGTCTAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAg                                               <  1:1839971/100‑1 (MQ=255)
    tggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGaaaa                                           <  1:64845/100‑1 (MQ=17)
    tggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGaaaa                                           <  2:1042324/100‑1 (MQ=18)
     ggtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaat                                          <  2:200311/100‑1 (MQ=18)
      gtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaata                                         <  2:1166192/100‑1 (MQ=18)
      gtCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAaata                                         >  2:632557/1‑100 (MQ=18)
          gATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTa                                     >  2:1263008/1‑100 (MQ=18)
           aTCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATCATTa                                     >  2:1722216/1‑99 (MQ=17)
           aTCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTAt                                    >  2:826824/1‑100 (MQ=33)
           aTCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTAt                                    >  1:434260/1‑100 (MQ=18)
              cAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGat                                 >  2:1800250/1‑100 (MQ=18)
                 cGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGattat                              >  1:1397205/1‑100 (MQ=18)
                    tGATATCCGACAGTAAGCAAAAATTCGATACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTAtttt                           <  1:870435/100‑1 (MQ=12)
                     gATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTAttttt                          >  2:1772433/1‑100 (MQ=18)
                      atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTtg                         <  2:938673/100‑1 (MQ=32)
                      atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtg                         <  2:247256/100‑1 (MQ=17)
                      atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtg                         <  1:241678/100‑1 (MQ=17)
                      atatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtg                         <  2:1637207/100‑1 (MQ=18)
                       tatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTtgt                        >  2:591498/1‑100 (MQ=33)
                       tatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgt                        >  1:342466/1‑100 (MQ=17)
                       tatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgt                        <  1:177611/100‑1 (MQ=17)
                       tatCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgt                        >  2:711565/1‑100 (MQ=18)
                        atCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATTTCATAATGCATCCCAGAAAATAATTATGATTATTTTtgtg                       <  1:1201216/100‑1 (MQ=12)
                        atCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgtg                       <  2:305320/100‑1 (MQ=18)
                        atCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTtgtg                       <  1:849034/100‑1 (MQ=18)
                          ccGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCa                     <  2:552724/100‑1 (MQ=18)
                          ccGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCa                     <  2:299938/100‑1 (MQ=18)
                          ccGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAAAAAATAATTATGATTATTTTTGTGCa                     >  2:1215924/1‑100 (MQ=255)
                           cGACAGTAAGCAAAAATTCGAGACATAGTGCGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGc                      >  2:1234100/1‑98 (MQ=255)
                                gTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTtg               <  1:1459457/100‑1 (MQ=32)
                                  aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta             >  2:161517/1‑100 (MQ=18)
                                  aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta             <  2:373635/100‑1 (MQ=18)
                                  aaGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta             >  1:1428576/1‑100 (MQ=18)
                                  aaGCAAAAATTCGAGAAATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta             <  2:536498/100‑1 (MQ=12)
                                  aaGCAAAAATGCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta             <  2:1133497/100‑1 (MQ=12)
                                  aaGAAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgta             <  1:296051/100‑1 (MQ=16)
                                   aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat            <  2:1467493/100‑1 (MQ=18)
                                   aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtat            >  2:43264/1‑100 (MQ=18)
                                   aGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGATGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTtgta             <  1:83278/99‑1 (MQ=17)
                                    gCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatg           >  2:31438/1‑100 (MQ=18)
                                    gCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatg           >  2:1290776/1‑100 (MQ=18)
                                    gCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgt              >  1:890923/1‑97 (MQ=18)
                                     cAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatgt          >  2:637601/1‑100 (MQ=18)
                                      aaaaaTTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTtgtatgta         >  2:1450363/1‑100 (MQ=18)
                                       aaaaTTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTac        >  1:1145187/1‑100 (MQ=18)
                                             cGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCCCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTc  <  1:727152/100‑1 (MQ=11)
                                                                                                |                                                
TCGGTGGTCTGATCCAGCGTTGATATCCGACAGTAAGCAAAAATTCGAGACATAGTGAGACTCTATACGAAATTGGTGTTCATATCATAATGCATCTCAGAAAATAATTATGATTATTTTTGTGCATATTTGTATGTACACGTTC  >  NC_000913/1433622‑1433766

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: