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,4160GA96.1% 157.3 / ‑0.6 51coding (657/657 nt)pphBserine/threonine‑specific protein phosphatase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/27);  minor base . (1/1);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

AACAAATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATA  >  NC_000913/2860283‑2860551
                                                                                                                                     |                                                                                                                                       
aacaaATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCa                                                                                                                                    <  2:1166624/139‑1 (MQ=255)
  caaATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  <  2:848075/139‑1 (MQ=255)
                         tataTTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAAc                                                                                                          <  2:264332/139‑1 (MQ=255)
                          atatTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAAcc                                                                                                         <  1:662836/139‑1 (MQ=255)
                            atTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                >  2:1609614/1‑115 (MQ=255)
                            atTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                <  1:1609614/115‑1 (MQ=255)
                            atTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCg                                                                                                       >  1:874914/1‑139 (MQ=255)
                                ggACATATGATGTTTGATAACATTCAGAGGTTCGATAACCAGATTTATAATGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCAgtgt                                                                                                   >  1:1081330/1‑139 (MQ=255)
                                     tatGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCCTGCTCTTCTGTGAAGCCTGa                                                                                                               <  2:1135516/122‑1 (MQ=255)
                                     tatGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGa                                                                                                               >  1:1135516/1‑122 (MQ=255)
                                       tgatgTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCa                                                                                            >  1:681684/1‑139 (MQ=255)
                                            tttGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  >  1:49966/1‑97 (MQ=255)
                                            tttGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  <  2:49966/97‑1 (MQ=255)
                                                atatCATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGcc                                                                                        >  1:1076965/1‑134 (MQ=255)
                                                aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGt                                                                                                             >  1:1468664/1‑113 (MQ=255)
                                                aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGt                                                                                                             <  2:1468664/113‑1 (MQ=255)
                                                aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGcc                                                                                        <  2:1076965/134‑1 (MQ=255)
                                                aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCattat                                                                                   >  2:1414939/1‑139 (MQ=255)
                                                             gTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                                      <  1:727707/139‑1 (MQ=255)
                                                              ttCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAAtctc                                                                                                                                     <  2:282904/76‑1 (MQ=255)
                                                              ttCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAAtctc                                                                                                                                     >  1:282904/1‑76 (MQ=255)
                                                               tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                >  2:242486/1‑80 (MQ=255)
                                                               tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                                <  1:242486/80‑1 (MQ=255)
                                                               tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCAt                                                                                                                 >  2:514806/1‑94 (MQ=255)
                                                               tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCAt                                                                                                                 <  1:514806/94‑1 (MQ=255)
                                                               tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGgtgt                                                                                                   <  1:1094465/108‑1 (MQ=255)
                                                               tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGgtgt                                                                                                   >  2:1094465/1‑108 (MQ=255)
                                                               tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGc                                                                        <  2:1635620/135‑1 (MQ=255)
                                                               tCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGc                                                                        >  1:1635620/1‑135 (MQ=255)
                                                               tCGATAACCAGCTTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCAGTGTTATCGCAGGCCATTATTCATTTTTCGCTAcc                                                                    <  2:1081330/139‑1 (MQ=255)
                                                                         gATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAAccc                                                                                                        >  2:923624/1‑93 (MQ=255)
                                                                         gATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAAccc                                                                                                        <  1:923624/93‑1 (MQ=255)
                                                                         gATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAagaga                                                          <  2:681684/139‑1 (MQ=255)
                                                                          aTTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCtg                                                                                                                         >  2:649418/1‑75 (MQ=255)
                                                                          aTTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCtg                                                                                                                         <  1:649418/75‑1 (MQ=255)
                                                                            ttATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       <  2:1584246/139‑1 (MQ=255)
                                                                                    aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  <  2:1364450/57‑1 (MQ=255)
                                                                                    aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATg                                                                                                                                  >  1:1364450/1‑57 (MQ=255)
                                                                                    aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGt                                                                                                             >  1:1589933/1‑77 (MQ=255)
                                                                                    aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGt                                                                                                             <  2:1589933/77‑1 (MQ=255)
                                                                                    aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAAccc                                                                                                        >  1:722514/1‑82 (MQ=255)
                                                                                    aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAAccc                                                                                                        <  2:722514/82‑1 (MQ=255)
                                                                                    aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTAcc                                                                    >  1:161205/1‑118 (MQ=255)
                                                                                    aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTAcc                                                                    <  2:161205/118‑1 (MQ=255)
                                                                                     tACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCa                                                                                       <  1:618330/98‑1 (MQ=255)
                                                                                     tACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCa                                                                                       >  2:618330/1‑98 (MQ=255)
                                                                                                       cgggcggCTGTCATTTTATAAAATAAAGTATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       >  1:1280170/1‑111 (MQ=255)
                                                                                                       cgggcggCTGTCATTTTATAAAATAAAGTATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGAt                                                       <  2:1280170/111‑1 (MQ=255)
                                                                                                         ggcggcTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAAcc                                                                                                         <  2:1087038/60‑1 (MQ=255)
                                                                                                         ggcggcTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAAcc                                                                                                         >  1:1087038/1‑60 (MQ=255)
                                                                                                                         taaaataaaGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt          <  2:350501/139‑1 (MQ=255)
                                                                                                                                 agtaatCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCAGGCAAACCGTCGTCTGtaata  >  1:1438823/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                       
AACAAATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCTGTAATA  >  NC_000913/2860283‑2860551

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