Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,381,542 A→G 7.7% pseudogene (722/969 nt) ycjV → putative ABC transporter ATP‑binding protein YcjV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,381,5420AG 7.7% 151.5 / 3.3 52pseudogene (722/969 nt)ycjVputative ABC transporter ATP‑binding protein YcjV
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base G (2/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CCCACGATCAGACCGAAGCGATGACCATGGCGACGCGGATTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATGGGAATACGACCGGAAGATATTCATCCGG  >  NC_000913/1381395‑1381691
                                                                                                                                                   |                                                                                                                                                     
cccACGATCAGACCGAAGCGATGACCATGGCGACGCGGATTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAGtt                                                                                                                                                     >  1:224586/1‑150 (MQ=255)
cccACGATCAGACCGAAGCGATGACCATGGCGACGCGGATTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAGtt                                                                                                                                                     >  1:219363/1‑150 (MQ=255)
 ccACGATCAGACCGAAGCGATGACCATGGCGACGCGGATTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAAttt                                                                                                                                                    <  2:671905/150‑1 (MQ=255)
 ccACGATCAGACCGAAGCGATGACCATGGCGACGCGGATTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAAttt                                                                                                                                                    <  2:649849/150‑1 (MQ=255)
  cACGATCAGACCGAAGCGATGACCATGGCGACGCGGATTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAAtttt                                                                                                                                                   >  2:29489/1‑150 (MQ=255)
  cACGATCAGACCGAAGCGATGACCATGGCGACGCGGATTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAAtttt                                                                                                                                                   >  2:1218777/1‑150 (MQ=255)
              gAAGCGATGACCATGGCGACGCGGATTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCAcg                                                                                                                                       >  1:596182/1‑150 (MQ=255)
                  cGATGACCATGGCGACGCGGATTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCAcgatcg                                                                                                                                   <  1:620495/150‑1 (MQ=255)
                               gACGCGGATTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAAtt                                                                                                                      <  1:477342/150‑1 (MQ=255)
                                      aTTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACg                                                                                                               >  2:854920/1‑150 (MQ=255)
                                      aTTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACg                                                                                                               >  1:19339/1‑150 (MQ=255)
                                         gtgATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGaa                                                                                                            >  1:109006/1‑150 (MQ=255)
                                          tgATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGaaa                                                                                                           >  1:1195400/1‑150 (MQ=255)
                                           gATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTGGCGGCACGATCGATGGCGATAAATTAGTTACGGAAAc                                                                                                          >  1:952973/1‑150 (MQ=255)
                                            aTTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACg                                                                                                         >  2:1409698/1‑150 (MQ=255)
                                            aTTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACg                                                                                                         >  2:168185/1‑150 (MQ=255)
                                             ttATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGc                                                                                                        <  2:249342/150‑1 (MQ=255)
                                                  aaaGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTaaa                                                                                                   >  1:95549/1‑150 (MQ=255)
                                                       cGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAAc                                                                                              <  1:938793/150‑1 (MQ=255)
                                                       cGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAAc                                                                                              <  2:769146/150‑1 (MQ=255)
                                                                tCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGa                                                                                     <  1:132593/150‑1 (MQ=255)
                                                                  aGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAag                                                                                   <  2:95553/150‑1 (MQ=255)
                                                                    cAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAagag                                                                                 >  1:296797/1‑150 (MQ=255)
                                                                         aGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAAtt                                                                            <  1:835304/150‑1 (MQ=255)
                                                                               gcCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCATCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGt                                                                      >  2:686312/1‑150 (MQ=255)
                                                                                  taaaaCCGTTTATAACCAACCCGCTAATATGTTTTTTTCCTGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCt                                                                   <  1:1218777/149‑1 (MQ=255)
                                                                                  gAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCt                                                                   <  1:29489/150‑1 (MQ=255)
                                                                                    aaaCCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAGTTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGa                                                                 <  2:224586/150‑1 (MQ=255)
                                                                                    aaaCCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAGTTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGa                                                                 <  2:219363/150‑1 (MQ=255)
                                                                                       ccGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGaaaa                                                              <  1:174682/150‑1 (MQ=255)
                                                                                        cGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAac                                                             <  1:1102556/150‑1 (MQ=255)
                                                                                            tataACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAg                                                         <  2:435727/150‑1 (MQ=255)
                                                                                            tataACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAg                                                         >  1:972377/1‑150 (MQ=255)
                                                                                            tataACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAg                                                         >  2:1430680/1‑150 (MQ=255)
                                                                                              taACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTAGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACAAAGGa                                                       >  1:928815/1‑150 (MQ=255)
                                                                                                     cccGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTg                                                <  1:1081614/150‑1 (MQ=255)
                                                                                                     cccGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTg                                                <  2:592127/150‑1 (MQ=255)
                                                                                                       cgcgAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCa                                              <  2:109006/150‑1 (MQ=255)
                                                                                                            atatGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAgg                                         <  1:318768/150‑2 (MQ=255)
                                                                                                               tgtttgttTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGcaca                                      >  1:1412978/1‑150 (MQ=255)
                                                                                                                  ttgtttCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCg                                   >  1:1267183/1‑150 (MQ=255)
                                                                                                                       tCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATg                              >  2:277473/1‑150 (MQ=255)
                                                                                                                        ccGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATgg                             <  2:1332999/150‑1 (MQ=255)
                                                                                                                         cGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAACCCCATCGTGATggg                            >  1:1051602/1‑150 (MQ=255)
                                                                                                                          ggATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCCGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATGGGa                           <  2:436666/150‑1 (MQ=255)
                                                                                                                                aTTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATGGGAATACGa                     >  1:997153/1‑150 (MQ=255)
                                                                                                                                    gATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATGGGAATACGACCgg                 >  1:1371513/1‑150 (MQ=255)
                                                                                                                                       cACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATGGGAATACGACCGGAAg              <  1:168185/150‑1 (MQ=255)
                                                                                                                                        aCCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATGGGAATACGACCGGAAGa             >  1:515939/1‑150 (MQ=255)
                                                                                                                                         ccAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATGGGAATACGACCGGAAGat            >  2:410759/1‑150 (MQ=255)
                                                                                                                                              gATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATGGGAATACGACCGGAAGATATTCa       >  1:1083183/1‑150 (MQ=255)
                                                                                                                                                 gAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATGGGAATACGACCGGAAGATATTCATcc    <  2:95549/150‑1 (MQ=255)
                                                                                                                                                  aaTTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTGCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGAGGGGAATACGACCGGAAGATATTCATCCg   <  2:296797/150‑1 (MQ=255)
                                                                                                                                                   aTTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATGGGAATACGACCGGAAGATATTCATCCgg  <  1:686312/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCCACGATCAGACCGAAGCGATGACCATGGCGACGCGGATTGTGATTATGAAAGACGGGATTGTTCAGCAAGTAGGTGCGCCGAAAACCGTTTATAACCAACCCGCGAATATGTTTGTTTCCGGATTTATTGGATCACCAGCGATGAATTTTATTCGCGGCACGATCGATGGCGATAAATTCGTTACGGAAACGCTTAAATTAACCATTCCCGAAGAGAAATTAGCGGTTCTGAAAACACAGGAAAGTTTGCATAAGCCCATCGTGATGGGAATACGACCGGAAGATATTCATCCGG  >  NC_000913/1381395‑1381691

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