Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,194,577 2 bp→GC coding (97‑98/774 nt) yjaD → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,194,5770CG99.2% 364.9 / NA 133R33A (CGG→GCG) yjaDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (68/64);  total (69/64)
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
*NC_0009134,194,5780GC99.2% 339.5 / ‑4.9 131R33A (CGG→GCG) yjaDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (67/63);  minor base A (0/1);  total (67/64)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 134 total aligned reads displayed.

GATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAACGGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGG  >  NC_000913/4194505‑4194645
                                                                        ||                                                                   
gATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTCGCCGGCGGGTGAAATGCCACATGGCGAACcgg                                                                   >  2:3667092/1‑76 (MQ=255)
    aCGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCaaa                                                               >  1:3617223/1‑76 (MQ=255)
    aCGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCaaa                                                               >  1:5294342/1‑76 (MQ=255)
    aCGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCaaa                                                               >  1:5171685/1‑76 (MQ=255)
    aCGGCTGGTGGGGCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGGGGCaaa                                                               >  1:912322/1‑76 (MQ=255)
      ggCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAtt                                                             >  1:3663646/1‑76 (MQ=255)
      ggCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAtt                                                             >  2:3048480/1‑76 (MQ=255)
       gCTGGGGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAttt                                                            <  1:2051209/76‑1 (MQ=255)
        cTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTc                                                           <  2:4517282/76‑1 (MQ=255)
        cTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGAGGCAAATTTc                                                           <  1:5409648/76‑1 (MQ=255)
         tggtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCg                                                          >  1:1103242/1‑76 (MQ=255)
          ggtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGa                                                         <  1:1523969/76‑1 (MQ=255)
          ggtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGa                                                         >  2:2554474/1‑76 (MQ=255)
          ggtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGa                                                         >  2:4620677/1‑76 (MQ=255)
           gtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGAt                                                        >  1:3433925/1‑76 (MQ=255)
            ggggTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGTCGAAGCGGCAAATTTCTATc                                                       <  2:4916699/75‑1 (MQ=255)
            ggggTCGCCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAAGTGCCATATGGCGTAGCGGCAAAGTTCGATc                                                       <  2:4909684/75‑1 (MQ=255)
              ggTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCtt                                                     >  1:378065/1‑76 (MQ=255)
              ggTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCtt                                                     <  2:4759026/76‑1 (MQ=255)
               gtcgtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAATTCGATCTtg                                                    >  2:5106911/1‑76 (MQ=255)
                tcgtcgGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtgt                                                   <  2:3762541/76‑1 (MQ=255)
                 cgtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtgtg                                                  <  2:1589822/76‑1 (MQ=255)
                 cgtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATATtgtg                                                  >  1:4806036/1‑76 (MQ=255)
                  gtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTgg                                                 <  2:3617223/76‑1 (MQ=255)
                  gtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTgg                                                 >  1:5744129/1‑76 (MQ=255)
                  gtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATGTCGATCTTGTgg                                                 >  1:733683/1‑76 (MQ=255)
                   tcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTggg                                                <  1:576991/76‑1 (MQ=255)
                   tcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTggg                                                >  1:5490069/1‑76 (MQ=255)
                    cAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCaaa                                                               <  2:1359274/60‑1 (MQ=255)
                    cAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCaaa                                                               >  1:1359274/1‑60 (MQ=255)
                    cAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGt                                               >  1:233437/1‑76 (MQ=255)
                    cAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGt                                               <  1:1564999/76‑1 (MQ=255)
                     aGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTc                                              >  2:4422902/1‑76 (MQ=255)
                     aGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTc                                              >  1:2108862/1‑76 (MQ=255)
                     aGCCATGAACAAAAATGATGGATGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTc                                              <  1:2800169/76‑1 (MQ=255)
                       ccATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGAAAATTTCGATCTTGTGGGTCAg                                            >  2:3725268/1‑76 (MQ=255)
                        cATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgc                                           <  2:1853327/76‑1 (MQ=255)
                        cATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgc                                           <  1:4027514/76‑1 (MQ=255)
                        cATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgc                                           >  2:3582480/1‑76 (MQ=255)
                         aTGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcg                                          >  2:2065491/1‑76 (MQ=255)
                         aTGAACAAAAATTATGGTTGCCGAAGGGAGAAGTGCCATATGGCGAAGCGGCAACTTTAGATCTTGTCGGTCAgcg                                          >  2:5219316/1‑76 (MQ=255)
                          tGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgc                                         >  2:3148745/1‑76 (MQ=255)
                          tGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgc                                         <  1:3148745/76‑1 (MQ=255)
                           gAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgcg                                        >  1:240571/1‑76 (MQ=255)
                           gAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgcg                                        <  1:5377791/76‑1 (MQ=255)
                            aaCAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgcgc                                       >  2:2938012/1‑76 (MQ=255)
                            aaCAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgcgc                                       >  1:2814592/1‑76 (MQ=255)
                            aaCAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgcgc                                       <  1:5127293/76‑1 (MQ=255)
                            aaCAAAAATTATGGCTGCCACAGATAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGTGTCAgctcgc                                       <  1:3667092/76‑1 (MQ=255)
                               aaaaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACt                                    >  1:1740122/1‑76 (MQ=255)
                                aaaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTa                                   >  1:3792462/1‑76 (MQ=255)
                                 aaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTAc                                  >  1:3392197/1‑76 (MQ=255)
                                  aaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACa                                 >  2:4808000/1‑76 (MQ=255)
                                     tATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAAt                              <  2:5171684/76‑1 (MQ=255)
                                     tATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAAt                              >  1:3113291/1‑76 (MQ=255)
                                     tATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAAt                              >  1:271323/1‑76 (MQ=255)
                                      agggTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATc                             <  2:22208/74‑1 (MQ=255)
                                       tggtggCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCg                            <  2:1702866/76‑1 (MQ=255)
                                       tGGTCGCCGATGGGCGAATTGCCATATGGAGAAGCGGTAAATGTCGATCTTGGGGGTCATCGCTCACTATAAATCg                            <  1:5795368/76‑1 (MQ=255)
                                            gCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGaa                       >  1:4659631/1‑76 (MQ=255)
                                             ccGAAGGGAGAATTGCCATATNGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAAt                      >  1:1300101/1‑76 (MQ=255)
                                             ccGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAAt                      <  2:3716541/76‑1 (MQ=255)
                                              cGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATg                     >  2:1573813/1‑76 (MQ=255)
                                              cGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATg                     >  2:2373671/1‑76 (MQ=255)
                                               gAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATgg                    >  2:4532792/1‑76 (MQ=255)
                                               gAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATgg                    <  2:2814592/76‑1 (MQ=255)
                                               gAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATgg                    >  1:3432426/1‑76 (MQ=255)
                                               gAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATgg                    >  2:4987822/1‑76 (MQ=255)
                                               gAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATgg                    >  1:4584095/1‑76 (MQ=255)
                                                aaGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGc                   >  2:2210830/1‑76 (MQ=255)
                                                aaGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGc                   <  1:181142/76‑1 (MQ=255)
                                                aaGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGc                   >  2:3242533/1‑76 (MQ=255)
                                                 aGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCa                  <  2:5804046/76‑1 (MQ=255)
                                                 aGGGAGAATTGCCATATGGCGAAGCGGCAAATGTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCa                  >  1:2291288/1‑76 (MQ=255)
                                                  gggAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGGGCGCACTACAAATCGGCGAATGGCAg                 <  2:3323085/76‑1 (MQ=255)
                                                   ggAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgg                >  1:1446277/1‑76 (MQ=255)
                                                   ggAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgg                <  1:1768536/76‑1 (MQ=255)
                                                   ggAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgg                <  2:5664746/76‑1 (MQ=255)
                                                    gagaATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAggg               <  1:5211542/76‑1 (MQ=255)
                                                     agaATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgggg              <  1:4172924/76‑1 (MQ=255)
                                                      gaATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAggggg             >  2:5471227/1‑76 (MQ=255)
                                                      gaATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAggggg             <  1:2065491/76‑1 (MQ=255)
                                                       aaTTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGa            >  1:4610323/1‑76 (MQ=255)
                                                       aaTTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGa            <  1:161288/76‑1 (MQ=255)
                                                       aaTTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGa            <  2:1480151/76‑1 (MQ=255)
                                                        atggCCATATGGCGAAGCGGCAAATTTCGACCTTCGGGGCCAGCGCGCTCTACAAAGCGGCGCACGACAGGGGCaa           >  2:597039/4‑76 (MQ=255)
                                                        aTTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGaa           >  1:3276910/1‑76 (MQ=255)
                                                         ttGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAAc          >  1:4155790/1‑76 (MQ=255)
                                                         ttGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAAc          >  1:4913707/1‑76 (MQ=255)
                                                          tGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAAcc         <  1:1432154/76‑1 (MQ=255)
                                                           gCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCt        >  2:4111480/1‑76 (MQ=255)
                                                           gCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCt        >  1:4528261/1‑76 (MQ=255)
                                                            ccATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTg       <  2:480561/76‑1 (MQ=255)
                                                            ccATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTg       <  1:1413878/76‑1 (MQ=255)
                                                            ccATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTg       >  2:5445155/1‑76 (MQ=255)
                                                            ccATATGGCGAAGCGGAAAATTTCGATCTTGTGGGTCAGCGCGCACGACAAATCGGCGAATGGCAGCGGGAAACta       >  2:2984878/1‑75 (MQ=255)
                                                             cATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGt      >  2:3952482/1‑76 (MQ=255)
                                                               tatGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGttt    >  2:1510838/1‑76 (MQ=255)
                                                                atGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTg   <  2:1300101/76‑1 (MQ=255)
                                                                 tGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTgg  <  2:2878554/76‑1 (MQ=255)
                                                                        ||                                                                   
GATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAACGGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGG  >  NC_000913/4194505‑4194645

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