Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,769,673 (A)8→7 pseudogene (1/111 nt) yibV → pseudogene, putative cryptic Rhs immunity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,769,6660A.98.3% 266.3 / ‑3.5 60pseudogene (322/327 nt)yibVpseudogene, putative cryptic Rhs immunity protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (33/26);  minor base T (0/1);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAGATGCTGGACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTAAAA  >  NC_000913/3769524‑3769803
                                                                                                                                              |                                                                                                                                         
aaTAGATGCTGGACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATa                                                                                                                                           <  1:163007/143‑1 (MQ=255)
          ggACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATaaaaaaatcat                                                                                                                                 >  1:293733/1‑139 (MQ=255)
           gACTTTAAGAATCATATACATGCTAATTATTAAATAAGGTAAAGATCATAAGACCATAACTGAAAAAACTAAATCCACTACAACAAACCGAAGATTGAATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCAtt                                                                                                                               <  2:392029/143‑1 (MQ=255)
              tttCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGat                                                                                                                            <  2:293733/143‑1 (MQ=255)
                     aTCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACATGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAAtt                                                                                                                     <  2:250864/143‑1 (MQ=255)
                      tCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTc                                                                                                                    <  1:1023291/143‑1 (MQ=255)
                      tCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTc                                                                                                                    >  1:699501/1‑143 (MQ=255)
                       ccTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCa                                                                                                                   >  1:116763/1‑143 (MQ=255)
                       acTATACATGCTATTTATTAAGTAAGAGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCa                                                                                                                   <  2:591811/142‑1 (MQ=255)
                        cTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAg                                                                                                                  >  2:372204/1‑143 (MQ=255)
                                gCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATACACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATAtt                                                                                                          >  1:99942/1‑143 (MQ=255)
                                 cTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATTAAAGATAAGAATACAAAAATCATTGATATAAATTCAGTTAATAttt                                                                                                         >  2:420323/1‑143 (MQ=255)
                                 cTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATAttt                                                                                                         >  2:417517/1‑143 (MQ=255)
                                   atttattAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGt                                                                                                       >  2:261272/1‑143 (MQ=255)
                                       attaAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGtg                                                                                                   <  1:964623/143‑1 (MQ=255)
                                          aagtaagGGAAAGATCATAAGACCATATCTGAAAAATCTAAATGCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAGGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGtgatg                                                                                                >  2:784984/1‑143 (MQ=255)
                                           agtaagGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGtgatga                                                                                               >  2:746162/1‑143 (MQ=255)
                                                 ggAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATAGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATGAAGTTAATATTTGTAGTGATGATaaaaa                                                                                         >  2:738888/1‑143 (MQ=255)
                                                   aaaGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGCTTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATCAGAAAAata                                                                                       >  2:715353/1‑143 (MQ=255)
                                                   aaaGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAata                                                                                       >  2:275130/1‑143 (MQ=255)
                                                          aTAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAAtt                                                                                <  1:975722/143‑1 (MQ=255)
                                                                    atCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAAt                                                                      >  1:582241/1‑143 (MQ=255)
                                                                      cggataaaTCTAAAACCGGATCAGCTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAg                                                                    <  1:818435/138‑1 (MQ=255)
                                                                              gcTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACtatta                                                            <  1:268276/142‑1 (MQ=255)
                                                                                tAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATg                                                          >  1:476971/1‑143 (MQ=255)
                                                                                tAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTAAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTATTATTATg                                                          >  2:186487/1‑143 (MQ=255)
                                                                                  aaTCCACTACAACATACCGATGAATGAATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGc                                                        <  2:854414/143‑1 (MQ=255)
                                                                                    tCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCa                                                      >  2:292161/1‑143 (MQ=255)
                                                                                    tCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCa                                                      >  2:92532/1‑143 (MQ=255)
                                                                                         tCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCAATtgat                                                 >  2:866329/1‑143 (MQ=255)
                                                                                           cAACTGGCCCCCGATTGCACTGAAAGAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGTCCATtgatga                                               <  1:68698/143‑1 (MQ=255)
                                                                                            aaCTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTgatgat                                              >  1:747884/1‑143 (MQ=255)
                                                                                             gcTTGCCGCTGATTGCAGTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGAtt                                             <  1:356410/142‑1 (MQ=255)
                                                                                             aCTTGCCGCTGATTGGATTGAAACAGTAAATAAAATAAAAGAGAAGAATAAAAAAATCATTGATATAAAGTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGAtt                                             <  1:227971/143‑1 (MQ=255)
                                                                                             aCTTGCCACTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGAtt                                             <  2:286543/143‑1 (MQ=255)
                                                                                                 gCCGCTGATTGCGTTGCCCCAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCtat                                         <  1:996399/143‑1 (MQ=255)
                                                                                                 gCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGATAATAGTACTATTATGGCCATTGATGATTCtat                                         <  1:261272/143‑1 (MQ=255)
                                                                                                  ccGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCtata                                        <  2:610707/143‑1 (MQ=255)
                                                                                                   cgatGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATaa                                       <  2:191829/140‑1 (MQ=255)
                                                                                                   cGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATaa                                       >  2:48922/1‑143 (MQ=255)
                                                                                                     cTGATTGGCTTGAAACAGGAAACAAAAGACAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAACTTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAg                                     <  2:892567/143‑1 (MQ=255)
                                                                                                     cTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAg                                     >  2:131263/1‑143 (MQ=255)
                                                                                                         ttGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCa                                 >  2:948786/1‑143 (MQ=255)
                                                                                                         ttGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCa                                 >  2:260005/1‑143 (MQ=255)
                                                                                                             aTTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATATGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGa                             >  2:208151/1‑143 (MQ=255)
                                                                                                             aTTGAAAAAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGa                             >  2:184388/1‑143 (MQ=255)
                                                                                                              ttaaaaCATAAAAAAAAATAAAAGTAAAAAATAAAAAAATCATTAATATAAATTCATTTAATATTTGTAGTGATTATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGat                            <  2:608214/140‑1 (MQ=255)
                                                                                                              ttGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGat                            <  1:666677/143‑1 (MQ=255)
                                                                                                                gAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGatga                          >  1:513543/1‑143 (MQ=255)
                                                                                                                     aGTAAATAAAATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGAt                     >  1:294365/1‑143 (MQ=255)
                                                                                                                       taaataaaATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGCTGCTTCTATAAAGTGCATCGATGAATGCTaa                   >  1:262977/1‑143 (MQ=255)
                                                                                                                       taaataaaATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATaa                   >  2:272921/1‑143 (MQ=255)
                                                                                                                       taaataaaATAAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATaa                   <  1:340488/143‑1 (MQ=255)
                                                                                                                        aaataaaataAAAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAAc                  <  2:726398/143‑1 (MQ=255)
                                                                                                                        aaagaaaagaaaagaTAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAAc                  <  2:690528/143‑1 (MQ=255)
                                                                                                                          ataaaataaaAGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACtt                >  2:434450/1‑143 (MQ=255)
                                                                                                                              aataaaaGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAAtt            >  1:1043159/1‑143 (MQ=255)
                                                                                                                                taaaaGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAAtttt          >  1:97595/1‑143 (MQ=255)
                                                                                                                                 aaaaGATAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATttta         <  2:116763/143‑1 (MQ=255)
                                                                                                                                       tAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTaaa   <  1:598676/143‑1 (MQ=255)
                                                                                                                                       tAAGAATAAAAAAATCATTGATATAAATTCAGTTAATATCTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCAAAAAAGTGCATCGATGAATGATAACTTAATTTTATTTaaa   <  2:1039312/143‑1 (MQ=255)
                                                                                                                                        aaGAATAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTaaaa  >  1:753005/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                         
AATAGATGCTGGACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTAAAA  >  NC_000913/3769524‑3769803

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