Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,632,837 A→G intergenic (‑552/‑53) ydfJ ← / → ynfT putative transporter YdfJ/Qin prophage; protein YnfT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,632,8370AG97.8% 104.2 / ‑4.3 46intergenic (‑552/‑53)ydfJ/ynfTputative transporter YdfJ/Qin prophage; protein YnfT
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (21/24);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

CTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAATTAGGGATTGGCATA  >  NC_000913/1632710‑1632973
                                                                                                                               |                                                                                                                                        
ctACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGt                                                                                                                              >  2:87282/1‑140 (MQ=25)
ctACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGt                                                                                                                              <  1:692412/140‑1 (MQ=18)
  aCATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCa                                                                                                                            >  1:576346/1‑140 (MQ=18)
   cATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCaa                                                                                                                           <  2:768684/140‑1 (MQ=17)
   cATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCaa                                                                                                                           <  1:651070/140‑1 (MQ=14)
     tAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACt                                                                                                                         <  2:309420/140‑1 (MQ=14)
     tAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACt                                                                                                                         <  2:327165/140‑1 (MQ=14)
        gATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACtata                                                                                                                      >  1:80550/1‑140 (MQ=18)
        gATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACtata                                                                                                                      <  1:3639/140‑1 (MQ=14)
         aTAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTatat                                                                                                                     >  1:67418/1‑140 (MQ=18)
         aTAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTatat                                                                                                                     >  1:559446/1‑140 (MQ=18)
          tAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATc                                                                                                                    <  1:87282/140‑1 (MQ=14)
          tAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATc                                                                                                                    <  2:194133/140‑1 (MQ=17)
          tAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATc                                                                                                                    <  2:751222/140‑1 (MQ=17)
           aGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCa                                                                                                                   >  1:162174/1‑140 (MQ=18)
             tttGCAGTTGCCGCTTCAGCTTGCGCCATAACCCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACtata                                                                                                                      <  2:467109/135‑1 (MQ=17)
             tttGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACtata                                                                                                                      >  1:467109/1‑135 (MQ=17)
             tttGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCata                                                                                                                 >  1:38113/1‑140 (MQ=18)
             tttGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCata                                                                                                                 >  1:537270/1‑140 (MQ=17)
                  aGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATaa                                                                                                            <  2:692150/140‑1 (MQ=14)
                  aGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATaa                                                                                                            <  1:534107/140‑1 (MQ=14)
                    ttGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAaca                                                                                                          <  2:330562/140‑1 (MQ=14)
                       ccGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCatatat                                                                                                              <  2:270595/133‑1 (MQ=14)
                       ccGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCatatat                                                                                                              >  1:270595/1‑133 (MQ=18)
                       ccGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACAt                                                                                                       >  1:511216/1‑140 (MQ=18)
                                                                       aGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTcc                                                       <  2:629598/140‑1 (MQ=18)
                                                                       aGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTcc                                                       <  2:162174/140‑1 (MQ=18)
                                                                            tCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATc                                                  >  1:735378/1‑140 (MQ=14)
                                                                             ccTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACAtctct                                                         >  1:303263/1‑132 (MQ=17)
                                                                             ccTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACAtctct                                                         <  2:303263/132‑1 (MQ=17)
                                                                               tGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCata                                               >  1:172591/1‑140 (MQ=17)
                                                                                      cTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCa                                        >  1:521429/1‑140 (MQ=14)
                                                                                        gAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCAtt                                      >  1:434540/1‑140 (MQ=17)
                                                                                                aaTGTTTATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCa                              >  1:494727/1‑140 (MQ=17)
                                                                                                    tttATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGcc                          <  2:110846/140‑1 (MQ=21)
                                                                                                    tttATCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGcc                          <  2:735378/140‑1 (MQ=21)
                                                                                                       aTCTACGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCta                       >  1:320890/1‑140 (MQ=18)
                                                                                                            cGAACTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATaaa                  >  1:266060/1‑140 (MQ=21)
                                                                                                              aaCTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAAt                <  1:689908/140‑1 (MQ=37)
                                                                                                              aaCTTTAAGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAAt                >  1:365443/1‑140 (MQ=18)
                                                                                                                  ttAAGAACACCCAGGATAAAAATTGTCAACTATATC‑‑ATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAATTAGGGa          <  1:483514/140‑1 (MQ=255)
                                                                                                                    aaGAACACCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAATTAGGGa          <  1:761659/140‑1 (MQ=37)
                                                                                                                         cacCCAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAATTAGGGATTGGc     >  2:644356/1‑140 (MQ=18)
                                                                                                                           cccAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAATTAGGGATTGGCat   >  2:107598/1‑140 (MQ=40)
                                                                                                                            ccAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAATTAGGGATTGGCata  <  2:434540/140‑1 (MQ=37)
                                                                                                                            ccAGGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAATTAGGGATTGGCata  <  2:67418/140‑1 (MQ=37)
                                                                                                                               |                                                                                                                                        
CTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGCCACCTGTTAGCATTCCTGTATACCTGAAACGACAATGTTTATCTACGAACTTTAAGAACACCCAAGATAAAAATTGTCAACTATATCATATATAACACATTACTAATTCGAGGCTATATGAACAGCATACTGATAATCACATCTCTCCTTATCATATTCAGCATTTTTAGTCATGCCCTAATAAAATTAGGGATTGGCATA  >  NC_000913/1632710‑1632973

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