Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,860,416 G→A *219* (TAG→TAA pphB → serine/threonine‑specific protein phosphatase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,860,4160GA98.6% 220.2 / ‑4.4 70*219* (TAG→TAApphBserine/threonine‑specific protein phosphatase 2
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (32/37);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

AACAAATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATA  >  NC_000913/2860283‑2860551
                                                                                                                                     |                                                                                                                                       
aacaaATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCGGTCATTTTATAAAATAAAGTAATCTCa                                                                                                                                    <  1:1200591/139‑1 (MQ=255)
  caaATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  >  2:1611958/1‑139 (MQ=255)
     aTAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCa                                                                                                                                    <  2:361896/134‑1 (MQ=255)
     aTAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCa                                                                                                                                    >  1:361896/1‑134 (MQ=255)
      tAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGctct                                                                                                                              >  1:1168653/1‑137 (MQ=255)
       aaaCGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGctctt                                                                                                                             <  2:1595476/139‑4 (MQ=255)
          cGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCtg                                                                                                                         <  2:72388/139‑1 (MQ=255)
           ggCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCtg                                                                                                                         <  2:1247631/139‑1 (MQ=255)
           ggCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCtgt                                                                                                                        <  2:169332/139‑1 (MQ=255)
                     attTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAg                                                                                                              <  2:47332/139‑1 (MQ=255)
                      ttttATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGt                                                                                                             <  2:1168653/139‑1 (MQ=255)
                          atatTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAAcc                                                                                                         <  1:1611958/139‑1 (MQ=255)
                            atTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAg                                                                                                              >  1:506328/1‑132 (MQ=255)
                            atTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAg                                                                                                              <  2:506328/132‑1 (MQ=255)
                                     tatGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTg                                                                                              >  1:1462342/1‑139 (MQ=255)
                                      atgatgTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAg                                                                                                                    <  1:134443/116‑1 (MQ=255)
                                      atgatgTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAg                                                                                                                    >  2:134443/1‑116 (MQ=255)
                                      atgatgTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGc                                                                                             >  1:116658/1‑139 (MQ=255)
                                            tttGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGgtg                                                                                                    >  2:760659/1‑126 (MQ=255)
                                            tttGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGgtg                                                                                                    <  1:760659/126‑1 (MQ=255)
                                                aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCattat                                                                                   >  2:122354/1‑139 (MQ=255)
                                                aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCattat                                                                                   >  1:525764/1‑139 (MQ=255)
                                                aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCattat                                                                                   >  2:302965/1‑139 (MQ=255)
                                                           aCGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGc                                                                        <  2:759961/139‑1 (MQ=255)
                                                            cGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCt                                                                       <  2:116658/139‑1 (MQ=255)
                                                             gTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                                      <  2:1462342/139‑1 (MQ=255)
                                                               tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCa                                                                                                                                    <  2:579588/76‑1 (MQ=255)
                                                               tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCa                                                                                                                                    >  1:579588/1‑76 (MQ=255)
                                                               tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                <  2:1590976/80‑1 (MQ=255)
                                                               tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                >  1:1590976/1‑80 (MQ=255)
                                                                cGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  <  2:1177995/77‑1 (MQ=255)
                                                                cGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGACGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCa                                                                   >  2:162169/1‑139 (MQ=255)
                                                                cGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCa                                                                   <  1:162169/139‑1 (MQ=255)
                                                                cGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCCCATg                                                                                                                                  >  1:1177995/1‑77 (MQ=255)
                                                                   tAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCa                                                                <  1:302965/139‑1 (MQ=255)
                                                                          aTTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                >  2:128351/1‑69 (MQ=255)
                                                                          aTTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                <  1:128351/69‑1 (MQ=255)
                                                                          aTTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAAccc                                                                                                        >  2:1126282/1‑92 (MQ=255)
                                                                          aTTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAAccc                                                                                                        <  1:1126282/92‑1 (MQ=255)
                                                                                    aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  >  2:427907/1‑57 (MQ=255)
                                                                                    aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  <  1:427907/57‑1 (MQ=255)
                                                                                    aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTc                                               >  2:794897/1‑139 (MQ=255)
                                                                                     tACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                                      <  1:597915/115‑1 (MQ=255)
                                                                                     tACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                                      >  2:597915/1‑115 (MQ=255)
                                                                                            tCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAgcgc                                       <  2:1017773/139‑1 (MQ=255)
                                                                                              gCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGgtgt                                                                                                   >  1:710331/1‑77 (MQ=255)
                                                                                              gCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGgtgt                                                                                                   <  2:710331/77‑1 (MQ=255)
                                                                                               cctaaCAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  <  2:1077928/43‑1 (MQ=255)
                                                                                               ccGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  >  1:940097/1‑46 (MQ=255)
                                                                                               ccGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  >  1:1077928/1‑46 (MQ=255)
                                                                                               ccGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  <  2:940097/46‑1 (MQ=255)
                                                                                               ccGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCa                                                                >  1:830208/1‑111 (MQ=255)
                                                                                               ccGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCa                                                                <  2:830208/111‑1 (MQ=255)
                                                                                                  aaCAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGc                                                                                                                   >  1:265444/1‑57 (MQ=255)
                                                                                                  aaCAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGc                                                                                                                   <  2:265444/57‑1 (MQ=255)
                                                                                                       cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCtgtg                                                                                                                       >  2:1075047/1‑48 (MQ=255)
                                                                                                       cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCtgtg                                                                                                                       <  1:1075047/48‑1 (MQ=255)
                                                                                                       cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                                      <  1:1072208/97‑1 (MQ=255)
                                                                                                       cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                                      >  2:1072208/1‑97 (MQ=255)
                                                                                                       cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCAc                            >  1:1135873/1‑139 (MQ=255)
                                                                                                         ggcggcTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTg                                                                                              >  2:1602304/1‑71 (MQ=255)
                                                                                                         ggcggcTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTg                                                                                              <  1:1602304/71‑1 (MQ=255)
                                                                                                              cTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCAct                           >  2:606733/1‑133 (MQ=255)
                                                                                                              cTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCAct                           <  1:606733/133‑1 (MQ=255)
                                                                                                                 tCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTa                                                                                                 <  2:602897/60‑1 (MQ=255)
                                                                                                                 tCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTa                                                                                                 >  1:602897/1‑60 (MQ=255)
                                                                                                                 tCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGc                                                                 <  2:1486039/92‑1 (MQ=255)
                                                                                                                 tCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGc                                                                 >  1:1486039/1‑92 (MQ=255)
                                                                                                                         taaaataaaGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt          <  1:794897/139‑1 (MQ=255)
                                                                                                                              taaaGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGta     <  2:349175/139‑1 (MQ=255)
                                                                                                                                 agtaatCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGtaata  >  1:1535166/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                       
AACAAATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATA  >  NC_000913/2860283‑2860551

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: