Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,994,415 G→A pseudogene (34/60 nt) pbl → pseudogene, peptidoglycan‑binding enzyme family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,994,4150GA98.7% 234.3 / NA 74pseudogene (34/60 nt)pblpseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (36/37);  minor base C (1/0);  total (37/37)
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

AACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGA  >  NC_000913/2994304‑2994550
                                                                                                               |                                                                                                                                       
aaCGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAg                                                                                                              <  1:1199720/139‑1 (MQ=255)
  cGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGcc                                                                                                            >  2:442683/1‑139 (MQ=255)
    tGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTa                                                                                                          >  1:1138009/1‑139 (MQ=255)
    tGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGAAAGCAATATGCCCAAAAATATATATTGTATATACTAGGCCTAATGAACTACACAATCGTAAGGCTTCAGCAAAATGAGCTTa                                                                                                          >  1:475558/1‑139 (MQ=255)
       cAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAAt                                                                                                                 <  1:1005714/129‑1 (MQ=255)
       cAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAAt                                                                                                                 >  2:1005714/1‑129 (MQ=255)
            ccTGGGGATTTATAAGGCGGATTTAGCAATTGATAACCAAAAAAAGCAACCGCAATATGCCGAAAAGTATATATTGTATACACCAGGCTTAATGAACTAAACAGTCGAAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  <  2:475558/139‑1 (MQ=255)
            ccTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  <  1:442683/139‑1 (MQ=255)
            ccTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  <  2:1101253/139‑1 (MQ=255)
                   aCTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCtttttttt                                                                                           <  2:847385/139‑1 (MQ=255)
                               ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                               <  1:1031286/139‑1 (MQ=255)
                               ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                               <  2:1138009/139‑1 (MQ=255)
                               ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                               >  2:1152508/1‑139 (MQ=255)
                                                   aaaaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  >  1:22540/1‑100 (MQ=255)
                                                   aaaaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  <  2:22540/100‑1 (MQ=255)
                                                   aaaaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGtt                                                           >  2:1340070/1‑139 (MQ=255)
                                                    aaaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  >  2:55677/1‑99 (MQ=255)
                                                    aaaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  <  1:55677/99‑1 (MQ=255)
                                                     aaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAg                                                                                                                       >  2:687938/1‑77 (MQ=255)
                                                     aaaaGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAg                                                                                                                       <  1:687938/77‑1 (MQ=255)
                                                           gaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                                   <  1:1340070/139‑1 (MQ=255)
                                                              aGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  >  2:93889/1‑89 (MQ=255)
                                                              aGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  <  1:93889/89‑1 (MQ=255)
                                                                                  atatTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAAt                                                                                                                 <  2:1076606/54‑1 (MQ=255)
                                                                                  atatTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAAt                                                                                                                 >  1:1076606/1‑54 (MQ=255)
                                                                                  atatTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTa                            <  1:170473/139‑1 (MQ=255)
                                                                                     ttGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCct                                                                                                    >  2:401501/1‑64 (MQ=255)
                                                                                     ttGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCct                                                                                                    <  1:401501/64‑1 (MQ=255)
                                                                                      tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCttt                                                                                                <  1:596712/67‑1 (MQ=255)
                                                                                      tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCttt                                                                                                >  2:596712/1‑67 (MQ=255)
                                                                                      tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAAttt                        <  1:1152508/139‑1 (MQ=255)
                                                                                       gTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACGCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATttc                       >  1:890087/1‑139 (MQ=255)
                                                                                        tatataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATttct                      >  1:132758/1‑139 (MQ=255)
                                                                                         atataCCATGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  <  2:885101/62‑1 (MQ=255)
                                                                                         atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  >  1:885101/1‑62 (MQ=255)
                                                                                            taCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCAt                  <  2:132758/139‑1 (MQ=255)
                                                                                             aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAg                                                                >  1:320144/1‑92 (MQ=255)
                                                                                             aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAg                                                                <  2:320144/92‑1 (MQ=255)
                                                                                             aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGAt                                                      >  2:1049051/1‑102 (MQ=255)
                                                                                             aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGAt                                                      <  1:1049051/102‑1 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCaaaa                                                                                                                  <  2:587721/37‑1 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCaaaa                                                                                                                  >  1:587721/1‑37 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  >  1:976619/1‑53 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  <  2:424734/53‑1 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  >  1:424734/1‑53 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  <  2:976619/53‑1 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                               <  2:57205/72‑1 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCaaa                                                                               >  1:57205/1‑72 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAacac                                                                            >  2:595205/1‑75 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAacac                                                                            <  1:595205/75‑1 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt                                                 <  1:481079/102‑1 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt                                                 >  2:481079/1‑102 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCa                                         >  2:1209435/1‑110 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCa                                         <  1:1209435/110‑1 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCAt                  >  1:136970/1‑133 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCAt                  <  2:136970/133‑1 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt            <  1:436063/139‑1 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt            >  1:782468/1‑139 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt            >  1:30482/1‑139 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt            >  1:2260/1‑139 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTATATAGTGCAAACACCCATACg                                                                     >  1:875963/1‑81 (MQ=255)
                                                                                                  gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTATATAGTGCAAACACCCATACg                                                                     <  2:875963/81‑1 (MQ=255)
                                                                                                      aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  <  1:279166/49‑1 (MQ=255)
                                                                                                      aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  <  1:934764/49‑1 (MQ=255)
                                                                                                      aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  >  2:934764/1‑49 (MQ=255)
                                                                                                      aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctct                                                                                                  >  2:279166/1‑49 (MQ=255)
                                                                                                       aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCa                                                                                 <  2:123164/65‑1 (MQ=255)
                                                                                                       aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCa                                                                                 >  1:123164/1‑65 (MQ=255)
                                                                                                       aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAATCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGGTATTGAgcc       >  2:795118/1‑137 (MQ=255)
                                                                                                       aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGtt                                                           <  1:156325/87‑1 (MQ=255)
                                                                                                       aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGtt                                                           >  2:156325/1‑87 (MQ=255)
                                                                                                            cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                            <  2:791786/97‑1 (MQ=255)
                                                                                                            cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGa  >  2:371227/1‑139 (MQ=255)
                                                                                                            cTAAACAATCGTAACGCTTCAGCAAAATGCGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                            >  1:791786/1‑97 (MQ=255)
                                                                                                               |                                                                                                                                       
AACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGA  >  NC_000913/2994304‑2994550

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