Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,966,751 C→T A112V (GCG→GTG)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,7510CT98.1% 186.3 / ‑3.9 53A112V (GCG→GTG) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (26/26);  total (27/26)
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.16e-01

CCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCG  >  NC_000913/3966608‑3966897
                                                                                                                                               |                                                                                                                                                  
cctacattCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGctgctg                                                                                                                                             >  2:453592/7‑151 (MQ=255)
      ttCCGCAACCAGCCCCTCCCCCGCCGGTCCTAAGACCATCTACGTTCCCCCGACACAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGa                                                                                                                                            <  2:830536/146‑1 (MQ=255)
        ccgaaaacacccccTCCCCCGCCGGTCCTGATGACCTCTCCTTTCCCCCTAGCAAAACCCGCCGTTTCACCCTCCGCACTGGTAATACCATCTCGGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGt                                                                                                                                        <  2:213424/136‑1 (MQ=255)
          gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACg                                                                                                                                   >  2:629081/1‑151 (MQ=255)
          gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACg                                                                                                                                   >  2:289588/1‑151 (MQ=255)
            aaccAGCTCCTCCCTCGCCGTTCTGGATGCCTTTTCGTTTTCCCATGGCAAAATCCGCCGTTCCACCCTCCGCACTGGTGATCCCATCTCGGATAGGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGaaa                                                                                                                                          <  2:602579/139‑1 (MQ=255)
            aGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGa                                                                                                                                  >  2:653450/1‑150 (MQ=255)
              aCAGCTCCTCCCTCGCCGGTCCTGATGACATCTCCTTTCCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGaa                                                                                                                                 <  1:370397/149‑1 (MQ=255)
               cAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGtt                                                                                                                              >  2:330966/1‑151 (MQ=255)
                   tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACt                                                                                                                          >  2:695522/1‑151 (MQ=255)
                    ccacccccGCCGGCCCTTATGACATCTCCGTTTCCCCAACCCAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATCCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTa                                                                                                                             <  2:861097/140‑1 (MQ=255)
                        ccTCGCGGGCCCTGATAACATCTCCTTTCCCCCTACCCAAACCCGCCGTTTCAACCCCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCg                                                                                                                       <  1:230294/149‑1 (MQ=255)
                                     gatgaCATCTCCGTTTCCCCTAGCCAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAc                                                                                                         <  1:220770/150‑1 (MQ=255)
                                     gatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACg                                                                                                        >  2:369850/1‑151 (MQ=255)
                                           aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                                                                                                   >  2:79949/1‑150 (MQ=255)
                                           aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                                                                                                   >  1:87781/1‑150 (MQ=255)
                                               aCGTTCCCCCAACCCAACCCCGCCTTTTCACCCTCCGCCCTGGTGATCCCATCTCTGTTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaaca                                                                                               <  1:105331/150‑1 (MQ=255)
                                                cGTTCCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaaca                                                                                               <  1:364318/149‑1 (MQ=255)
                                                  ttcccccTAGCCAAATCCGCCGTTTCACCCCCCGCACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaaa                                                                                            <  1:3604/147‑1 (MQ=255)
                                                    tCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc                                                                                         >  1:568441/1‑151 (MQ=255)
                                                     ccccTACCCAAACCCGCCTTTTCACCCGCCGCACTGTTGATACCATCTCGGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt                                                                                           <  2:863730/148‑1 (MQ=255)
                                                      cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGCTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCct                                                                                        >  2:75091/1‑150 (MQ=255)
                                                              aattccacagTTACAACCACCGAGCAGGTGAGACCATGTTTGGAAGGATTGACCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAGCCTCTTTg                                                                                   <  1:592357/138‑1 (MQ=255)
                                                              aaaTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga                                                                                >  1:806537/1‑150 (MQ=255)
                                                                acccgccgTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGaa                                                                               <  2:222730/147‑1 (MQ=255)
                                                                       gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAAcccc                                                                      >  1:888453/1‑151 (MQ=255)
                                                                                 ccGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACg                                                              <  1:729402/149‑1 (MQ=255)
                                                                                    cACTGGTGATCCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAAc                                                         <  1:653597/151‑1 (MQ=255)
                                                                                       tGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAActct                                                      >  1:584524/1‑151 (MQ=255)
                                                                                           gATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTc                                                  >  2:745408/1‑151 (MQ=255)
                                                                                             tACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTg                                                >  1:797540/1‑151 (MQ=255)
                                                                                              aCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGc                                               >  2:801522/1‑151 (MQ=255)
                                                                                               ccATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCg                                              <  1:140051/151‑1 (MQ=255)
                                                                                                  tctcTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTAt                                           >  1:797201/1‑151 (MQ=255)
                                                                                                      tgttAAGTTCCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATgg                                         <  1:815551/146‑1 (MQ=255)
                                                                                                       ggTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTa                                            <  2:631107/145‑1 (MQ=255)
                                                                                                          aaGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGt                                    <  1:377644/150‑1 (MQ=255)
                                                                                                           aGATTCGACCGCCGAAAGAAGGTGAACGCTCTTTCGTGCTGCTGAAAGTTACCGAAGTTAACTTCGACAAACCTGCAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGCACGTgg                                  >  2:626550/1‑151 (MQ=255)
                                                                                                            gATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTgg                                  >  2:688851/1‑150 (MQ=255)
                                                                                                             aTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGt                                 <  1:244580/150‑1 (MQ=255)
                                                                                                              ttCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTa                                <  1:73472/150‑1 (MQ=255)
                                                                                                               tCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTaa                               <  1:458405/150‑1 (MQ=255)
                                                                                                               tCGCCCGCCGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAAc                              >  2:805142/1‑151 (MQ=255)
                                                                                                                cgcccgccGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAAc                              <  1:803662/150‑1 (MQ=255)
                                                                                                                cgcccgccGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAAGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACg                             >  2:127911/1‑151 (MQ=255)
                                                                                                                  cccgccGAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAAc                              <  2:583067/148‑1 (MQ=255)
                                                                                                                   ccgccgAAAGAAGGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTCAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGtt                          >  2:450736/1‑151 (MQ=255)
                                                                                                                             aagGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGa                   >  2:546962/1‑148 (MQ=255)
                                                                                                                              agGTGAACGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGa                   >  2:643344/1‑147 (MQ=255)
                                                                                                                                     cGCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTg        >  1:505659/1‑151 (MQ=255)
                                                                                                                                      gCTATTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGc       <  2:346994/151‑1 (MQ=255)
                                                                                                                                         aTTTTGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTcg    >  1:94934/1‑151 (MQ=255)
                                                                                                                                            ttGTGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTcgcg  <  1:490690/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCG  >  NC_000913/3966608‑3966897

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