Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 584,262 (A)7→8 coding (582/750 nt) appY → global transcriptional activator; DLP12 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913584,2551.A89.8% 192.3 / 15.5 59coding (575/750 nt)appYglobal transcriptional activator; DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base . (1/5);  new base A (29/24);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTA  >  NC_000913/584156‑584347
                                                                                                    |                                                                                              
cGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGc                                                                                                 >  2:1415944/1‑100 (MQ=255)
cGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGc                                                                                                 >  2:1265278/1‑100 (MQ=255)
cGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGc                                                                                                 <  1:99518/100‑1 (MQ=255)
 gaaattaattTAAACAAGTAAAATTTAAAAAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGc                                                                                                 <  2:85590/97‑1 (MQ=255)
         tttATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaaaac                                                                                       >  1:1471387/1‑98 (MQ=255)
            ataCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGAGACTAGGATGAGGTATGCaaaaaaaactca                                                                                    >  2:949281/1‑95 (MQ=255)
              aCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATa                                                                                 >  2:1532352/1‑100 (MQ=255)
              aCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATa                                                                                 <  1:1284043/100‑1 (MQ=255)
                 aGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTAAAATATAAAGAGTTAATAGAATGAGGTATGCaaaaaaaactca                                                                                    >  2:1425690/1‑90 (MQ=255)
                       agagaaataaAAAAAAGATTAAGGGATGAAGGAACGTCATAGACTGAAATATTGAGAGATACTATGATAAAAAAAAAAAAAAAAA‑‑CTCATAACTTCAAAc                                                                        <  2:1365342/95‑1 (MQ=255)
                       aGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAAc                                                                        <  2:36778/100‑1 (MQ=255)
                            aataaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTa                                                                   >  1:1178302/1‑100 (MQ=255)
                             aaaaaaaaaaaaTTAACAGATAAACAAACCTCATTTACAAAAATATTCAGATATACCAGGATAAAATAAACAAAAAAAA‑‑CTCATAACTTCAAACTCTTAt                                                                  <  2:596891/98‑1 (MQ=255)
                              taaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAA‑‑CTCATAACTTCAAACTCTTATTc                                                                <  2:492438/100‑1 (MQ=255)
                               aaaaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTc                                                                >  2:197600/1‑100 (MQ=255)
                                  aaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAAAAAAAAAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTAt                                                             <  2:1289377/100‑1 (MQ=255)
                                   aaaaGATTAAGGGATGAAGGAACGTCATTTACGGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATc                                                            <  1:377993/100‑1 (MQ=255)
                                    aaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCa                                                           <  1:988990/100‑1 (MQ=255)
                                    aaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAAGAAGAAAAATACAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCa                                                           <  1:918943/100‑1 (MQ=255)
                                     aaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCaa                                                          <  1:235955/100‑1 (MQ=255)
                                        aTTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGt                                                       >  1:125197/1‑100 (MQ=255)
                                         ttAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTc                                                      >  2:66644/1‑100 (MQ=255)
                                         ttAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTc                                                      >  2:1179929/1‑100 (MQ=255)
                                         ttAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCg                                                     <  1:590418/100‑1 (MQ=255)
                                           aaGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTa                                                   <  1:1082927/100‑1 (MQ=255)
                                             gggATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAg                                                  <  1:555149/100‑1 (MQ=255)
                                              ggATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGc                                                 >  2:110020/1‑100 (MQ=255)
                                              ggATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGc                                                 >  1:713866/1‑100 (MQ=255)
                                              ggATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGAAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGc                                                 <  2:1472811/100‑1 (MQ=255)
                                               gATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGcc                                                >  2:1352090/1‑100 (MQ=255)
                                                 tGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCa                                              >  1:546342/1‑100 (MQ=255)
                                                  gAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCTATGTCGTAGCCCAg                                             >  2:1468416/1‑100 (MQ=255)
                                                  gAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAg                                             >  2:339205/1‑100 (MQ=255)
                                                   aaGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCACACTCTTATTCTATCAATGTCGTAGCCCAGa                                            >  2:369576/1‑100 (MQ=255)
                                                       aaCGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtg                                        >  1:904777/1‑100 (MQ=255)
                                                       aaCGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtg                                        >  1:772668/1‑100 (MQ=255)
                                                        aCGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtgt                                       >  1:649785/1‑100 (MQ=255)
                                                        aCGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtgt                                       <  1:470837/100‑1 (MQ=255)
                                                          gttcttcaCTGAAATATTGAGAGATACTAGAATGTAGTGAGAAAAAAAATCTCATAACCTCAAGCTCTTATTCTAACAATGTCGTAGCCCAGAAAtgtg                                      <  2:1416633/92‑1 (MQ=255)
                                                          gTCATTTACTGAAATATTGAGATATACTAGTATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATTTCGTAGCCCAGAAAtgt                                       <  2:201091/98‑1 (MQ=255)
                                                          gTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTgg                                     <  1:308110/100‑1 (MQ=255)
                                                           tCATTTACTGAAATATTGAGATATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGc                                    <  1:1378452/100‑1 (MQ=255)
                                                           tCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGc                                    <  1:197600/100‑1 (MQ=255)
                                                            cATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCt                                   <  2:1089816/99‑1 (MQ=255)
                                                             aTTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCta                                  <  1:1523540/100‑1 (MQ=255)
                                                              tttACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCtat                                 >  1:745964/1‑100 (MQ=255)
                                                                  cTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAAc                              >  1:1223460/1‑99 (MQ=255)
                                                                   tGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAg                            <  1:1461370/100‑1 (MQ=255)
                                                                    gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                           >  1:863022/1‑100 (MQ=255)
                                                                    gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACCTAAAAACTCCAAACCTTAATTCTATAATTTATTTACCCAGAAAAGTgggatt                                 >  1:478869/1‑90 (MQ=255)
                                                                     aaaTATTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTa                          <  1:1485958/100‑1 (MQ=255)
                                                                       atatTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACt                        <  2:1564187/100‑1 (MQ=255)
                                                                         atTGAGAGATACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttc                      >  1:790192/1‑100 (MQ=255)
                                                                            gagagaTACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcata                   >  1:800557/1‑100 (MQ=255)
                                                                               agaTACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatattt                >  1:489252/1‑100 (MQ=255)
                                                                               agaTACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatattt                >  1:280997/1‑100 (MQ=255)
                                                                                gaTACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttc               <  2:1356763/100‑1 (MQ=255)
                                                                                  tACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcat             >  2:1196258/1‑100 (MQ=255)
                                                                                  tACTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttca              <  1:1265283/99‑1 (MQ=255)
                                                                                   aCTAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcata            <  1:803157/100‑1 (MQ=255)
                                                                                     tAGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAtg          <  1:1223084/100‑1 (MQ=255)
                                                                                      aGGATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAtgt         >  1:18515/1‑100 (MQ=255)
                                                                                        gATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGc       >  2:570748/1‑100 (MQ=255)
                                                                                        gATGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGc       >  2:570761/1‑100 (MQ=255)
                                                                                          tGAGGTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCAt     >  2:1134590/1‑100 (MQ=255)
                                                                                             ggTATGCAAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTa  >  2:1407170/1‑100 (MQ=255)
                                                                                                    |                                                                                              
CGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAA‑‑CTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTA  >  NC_000913/584156‑584347

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: