Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 575,915 T→A pseudogene (911/1040 nt) nmpC ← DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913575,9150TA97.4% 259.2 / ‑1.2 76pseudogene (911/1040 nt)nmpCDLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (31/43);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

TTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCA  >  NC_000913/575818‑576010
                                                                                                 |                                                                                               
tttGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACAGATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAAcc                                                                                               <  2:664695/100‑1 (MQ=255)
   gTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTa                                                                                            <  2:507160/100‑1 (MQ=255)
   gTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTa                                                                                            >  1:604925/1‑100 (MQ=255)
    tAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTAc                                                                                           >  1:72159/1‑100 (MQ=255)
     aaTCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTAcc                                                                                          <  2:781846/100‑1 (MQ=255)
      aTCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCt                                                                                         >  2:280101/1‑100 (MQ=255)
        cAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCttt                                                                                       <  2:249511/100‑1 (MQ=255)
        cAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCttt                                                                                       >  1:348858/1‑100 (MQ=255)
           cGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGa                                                                                    <  1:122720/100‑1 (MQ=255)
                gTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCa                                                                               <  2:816407/100‑1 (MQ=255)
                 tAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAg                                                                              <  2:1394189/100‑1 (MQ=255)
                 tAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAg                                                                              <  2:196374/100‑1 (MQ=255)
                 tAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAg                                                                              <  1:971220/100‑1 (MQ=255)
                  aGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAgg                                                                             <  1:1573452/100‑1 (MQ=255)
                   gACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGt                                                                            >  1:199801/1‑100 (MQ=255)
                     cATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTaa                                                                          <  2:468642/100‑1 (MQ=255)
                      atatTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAg                                                                         >  1:1489/1‑100 (MQ=255)
                       tatTTTTGTTCAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCTAAGAACCTTACCTTTAGATTGCAGGTAAGc                                                                        <  2:1031195/100‑1 (MQ=255)
                        atTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCa                                                                       <  1:285140/100‑1 (MQ=255)
                        atTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCa                                                                       <  1:1098163/100‑1 (MQ=255)
                        atTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCa                                                                       >  2:311705/1‑100 (MQ=255)
                         ttttttttGTTGTTATAGGTTGCACATACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCaa                                                                      >  2:23879/1‑100 (MQ=255)
                          ttttGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAac                                                                     <  2:367296/100‑1 (MQ=255)
                            ttgttgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAacgg                                                                   <  1:609243/100‑1 (MQ=255)
                             tgttgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAacgga                                                                  <  2:1264511/100‑1 (MQ=255)
                              gttgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAacggac                                                                 <  2:541435/100‑1 (MQ=255)
                               ttgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAacggacg                                                                <  2:604922/100‑1 (MQ=255)
                                tgAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAacggacgg                                                               >  2:795686/1‑100 (MQ=255)
                                  aaGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACggacggac                                                             >  1:1571694/1‑100 (MQ=255)
                                   aGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGgacggacg                                                            <  2:30704/100‑1 (MQ=255)
                                   aGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGgacggacg                                                            >  1:58622/1‑100 (MQ=255)
                                   aGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGgacggacg                                                            >  2:221591/1‑100 (MQ=255)
                                   aGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGgacggacg                                                            <  2:1454890/100‑1 (MQ=255)
                                    gTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTTATCGCCCCATACCCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGc                                                           <  2:756689/100‑1 (MQ=255)
                                    gTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGgacggacgn                                                           <  2:756645/100‑2 (MQ=255)
                                     tAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCa                                                          <  1:601279/100‑1 (MQ=255)
                                     tAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCa                                                          <  1:495104/100‑1 (MQ=255)
                                        tAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGAc                                                       <  1:795691/100‑1 (MQ=255)
                                           gTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCg                                                     <  2:1023233/99‑1 (MQ=255)
                                                 ccTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGa                                              >  2:434435/1‑100 (MQ=255)
                                                 ccTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGa                                              <  2:348857/100‑1 (MQ=255)
                                                   tACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAAc                                            >  2:1554968/1‑100 (MQ=255)
                                                   tACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAAc                                            <  2:1010970/100‑1 (MQ=255)
                                                   tACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAAc                                            >  1:539251/1‑100 (MQ=255)
                                                   tACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAAc                                            <  1:641201/100‑1 (MQ=255)
                                                      aTCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGa                                         >  1:1448316/1‑100 (MQ=255)
                                                      aTCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGa                                         >  1:1448339/1‑100 (MQ=255)
                                                          aCATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATAtt                                     >  2:1529100/1‑100 (MQ=255)
                                                           cATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATtg                                    >  2:1061761/1‑100 (MQ=255)
                                                           cATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATtg                                    <  1:609032/100‑1 (MQ=255)
                                                            atatTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATtgt                                   <  1:431130/100‑1 (MQ=255)
                                                            atatTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATtgt                                   <  1:55386/100‑1 (MQ=255)
                                                             tatTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATtgtg                                  >  1:1199095/1‑100 (MQ=255)
                                                              atTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGc                                 >  2:1394615/1‑100 (MQ=255)
                                                              atTTGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGc                                 <  1:23879/100‑1 (MQ=255)
                                                                 tGACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAAc                              <  2:752690/100‑1 (MQ=255)
                                                                  gACTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTAAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACa                             >  2:857605/1‑100 (MQ=255)
                                                                    cTAAGTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGc                           >  2:893580/1‑100 (MQ=255)
                                                                        gTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCg                       >  1:952930/1‑100 (MQ=255)
                                                                        gTCCTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCg                       <  1:1020011/100‑1 (MQ=255)
                                                                           cTGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAg                    <  1:528334/100‑1 (MQ=255)
                                                                            tGATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGCTATTGTGCAACAGCTTCGAGGt                   >  1:723784/1‑100 (MQ=255)
                                                                             gATCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGtt                  >  1:1345397/1‑100 (MQ=255)
                                                                              aTCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGttt                 >  2:1158434/1‑100 (MQ=255)
                                                                              aTCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGttt                 <  1:969773/100‑1 (MQ=255)
                                                                               tCGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGtttt                <  1:1059571/100‑1 (MQ=255)
                                                                                cGCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTg               <  2:58622/100‑1 (MQ=255)
                                                                                 gCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTgg              >  2:1460484/1‑100 (MQ=255)
                                                                                 gCCCCATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTgg              <  2:1460080/100‑1 (MQ=255)
                                                                                    ccATACTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCt           <  2:883956/100‑1 (MQ=255)
                                                                                       tACTCCAAGAACCTTACCTTTAGATTGCAGGTAGGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTGTTGGGCTtta        >  2:1035143/1‑100 (MQ=255)
                                                                                        aCTCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttat       >  2:897395/1‑100 (MQ=255)
                                                                                          tCCAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTttatta     <  1:273053/100‑1 (MQ=255)
                                                                                           ccAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAg    >  1:541966/1‑100 (MQ=255)
                                                                                           ccAAGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAg    >  1:1496111/1‑100 (MQ=255)
                                                                                             aaGAACCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCa  >  2:1162427/1‑100 (MQ=255)
                                                                                                 |                                                                                               
TTTGTAATCAACGAAAGTAGACATATTTTTGTTGAAGTAATAGGTTGCACCTACATCAACATATTTGACTAAGTCCTGATCGCCCCATACTCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCA  >  NC_000913/575818‑576010

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