Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,572,630 (T)6→7 coding (217/1518 nt) yjiT → putative uncharacterized protein YjiT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,572,6241.T98.0% 170.7 / ‑3.5 49coding (211/1518 nt)yjiTputative uncharacterized protein YjiT
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (25/23);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

ACCGGATGGCCGTCCTCTTTATGAATATCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCGTACCGAAAGGGATGGA  >  NC_000913/4572494‑4572754
                                                                                                                                    |                                                                                                                                  
aCCGGATGGCCGTCCTCTTTATGAATATCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCtttttttgc                                                                                                                           >  2:254264/1‑137 (MQ=255)
 ccGGATGGCCGTCCTCTTTATGAATATCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCtttttttgct                                                                                                                          <  1:845126/140‑5 (MQ=255)
  cGGATGGCCGTCCTCTTTATGAATATCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCtttttttgctc                                                                                                                         >  1:664532/1‑135 (MQ=255)
         ccGTCCTCTTTATGAATATCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTgg                                                                                                                   <  2:34280/140‑1 (MQ=255)
              ctctTTATGAATATCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCg                                                                                                              <  1:405584/140‑1 (MQ=255)
              ctctTTATGAATATCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGACTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCg                                                                                                              <  2:813009/140‑1 (MQ=255)
               tctTTATGAATATCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGc                                                                                                             <  2:286427/140‑1 (MQ=255)
               tctTTATGAATATCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGc                                                                                                             <  1:518586/140‑1 (MQ=255)
                       aaTATCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTa                                                                                                     >  1:730237/1‑140 (MQ=255)
                         tatCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATg                                                                                                   <  1:881505/140‑1 (MQ=255)
                          atCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCTTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGa                                                                                                  >  2:864509/1‑140 (MQ=255)
                           tCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGaa                                                                                                 >  1:721820/1‑140 (MQ=255)
                             aTGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACg                                                                                               >  2:292954/1‑140 (MQ=255)
                              tGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACgg                                                                                              >  2:421797/1‑140 (MQ=255)
                                cAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCa                                                                                            <  1:421797/140‑1 (MQ=255)
                                cAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCa                                                                                            <  2:721820/140‑1 (MQ=255)
                                   ccAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATg                                                                                         >  2:117660/1‑140 (MQ=255)
                                   ccAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATg                                                                                         >  1:875797/1‑140 (MQ=255)
                                      atgatgAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCgg                                                                                      >  2:888488/1‑140 (MQ=255)
                                       tgatgaATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGa                                                                                     >  1:688482/1‑140 (MQ=255)
                                         atgaATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATg                                                                                   <  2:178324/140‑1 (MQ=255)
                                             atatACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACg                                                                               <  1:349943/140‑1 (MQ=255)
                                             atatACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACg                                                                               <  2:664532/140‑1 (MQ=255)
                                               ataCCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTg                                                                             <  1:877079/140‑1 (MQ=255)
                                               ataCCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTg                                                                             <  1:292954/140‑1 (MQ=255)
                                                taCCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTgg                                                                            <  1:888488/140‑1 (MQ=255)
                                                     aaTTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTgg                                                                            <  2:46272/135‑1 (MQ=255)
                                                     aaTTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTgg                                                                            >  1:46272/1‑135 (MQ=255)
                                                     aaTTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGAccc                                                                       <  1:182976/140‑1 (MQ=255)
                                                             cAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCtatatata                                                               >  2:616925/1‑140 (MQ=255)
                                                                cTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATaaaa                                                            >  2:814661/1‑140 (MQ=255)
                                                                      cGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAg                                                      >  2:596905/1‑140 (MQ=255)
                                                                      cGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAg                                                      <  1:254264/140‑1 (MQ=255)
                                                                        tGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAggg                                                    >  1:813912/1‑140 (MQ=255)
                                                                                   aaTCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCAc                                         <  2:813912/140‑1 (MQ=255)
                                                                                          aTCTAATATATGTAATAGGGACTTTGCTGCTTTTTTTGTTCTTTTTTTTCTGGGAGTGGTACCGGAGGGATTATGAACGGCAAGGGGGAGGGAGGTGGGACCTTATATAAAAAAAATTAGGTATTTCATTCACCGCcaca                                  <  2:43946/140‑1 (MQ=255)
                                                                                           tCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACag                                 >  2:814008/1‑140 (MQ=255)
                                                                                             tAATATATGTAATAGGGACTTTGCTGCTTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACagag                               <  1:864509/140‑1 (MQ=255)
                                                                                             tAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTTTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACagag                               <  2:875797/140‑1 (MQ=255)
                                                                                             tAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACagag                               <  1:596905/140‑1 (MQ=255)
                                                                                               atatatGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCt                             <  1:616925/140‑1 (MQ=255)
                                                                                                     gTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTAc                       >  2:360081/1‑140 (MQ=255)
                                                                                                         tAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATc                   >  2:48204/1‑140 (MQ=255)
                                                                                                          aGGGACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCg                  >  1:118682/1‑140 (MQ=255)
                                                                                                            ggACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCGTa                >  2:900171/1‑140 (MQ=255)
                                                                                                            ggACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCGTa                >  2:668492/1‑140 (MQ=255)
                                                                                                            ggACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCGTa                >  1:918235/1‑140 (MQ=255)
                                                                                                            ggACTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCGTa                >  1:843835/1‑140 (MQ=255)
                                                                                                              aCTTTGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCGTAcc              >  2:505823/1‑140 (MQ=255)
                                                                                                                tttGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCGTACCGa            >  2:753666/1‑140 (MQ=255)
                                                                                                                tttGCTGCCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCGTACCGa            >  1:476301/1‑140 (MQ=255)
                                                                                                                      gcCTTTTTTGTTTTTTTTTTGCTCTGGGTGGTCCCGCAGGGATTATGAACGGCAATGCGGAGGGACGGGGGCCCCTATAAAAAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCGTACCGAAAGGGa      <  2:688482/140‑1 (MQ=255)
                                                                                                                      gcCTGTTTTGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCGTACCGAAAGGGa      <  2:118682/140‑1 (MQ=255)
                                                                                                                          tttttGTTCTTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCGTACCGAAAGGGATGGa  <  1:48204/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                  
ACCGGATGGCCGTCCTCTTTATGAATATCATGCAACCAATGATGAATATACCCAATTAACGCAGCTACTCCGTGCAGTCGGTCAATCACAATCTAATATATGTAATAGGGACTTTGCTGCCTGTTTTGTTCTTTTTTGCTCTGAGTGGTACCGCAGGGATTATGAACGGCAATGCGGATGGACGTGGGACCCTATATATAAAAAAATAGGGATTTCATTCACCGCCACAGAGCTTGGTACTATCGTACCGAAAGGGATGGA  >  NC_000913/4572494‑4572754

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