Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,966,530 2 bp→TT coding (114‑115/1260 nt) rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,5300CT87.5% 61.3 / ‑2.2 24I38I (ATC→ATTrhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (12/9);  total (14/10)
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_0009133,966,5310CT91.7% 56.5 / ‑3.2 24L39L (CTG→TTG) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (12/10);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACT  >  NC_000913/3966300‑3966777
                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                                      
gtgACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAaga                                                                                                                                                                                                                                      <  1:203271/250‑1 (MQ=255)
              aaaaGACATGAGTCGATGACCGTAAACAGACATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAACCCTGGCTCGTCTGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGAACGCAAAGAGGGGCGAAGAAATCt                                                                                                                                                                                                                       >  2:222268/1‑251 (MQ=255)
                    cATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTg                                                                                                                                                                                                                 <  1:31923/251‑1 (MQ=255)
                          tcatgacacgaaaaaatgCATTGATGACCCTGCCATACCATGCACAACATTAAGTTCGAGATTGACCCCAAGTTGAATCACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCTGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGTGCTGGAAAACCTTGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCTGGAATCTGCACGCAAAGAGTGGCGAAGATACCTTTGGTgctgtc                                                                                                                                                                                                            <  2:43907/233‑6 (MQ=255)
                                                 atgatCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATAACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAACTTATTTTTGCCATCCTGAACAAGCACGCACATGGTGGCGAGGATATCTTTGGTGGTGTCGTaaaggggagatttgaggatggt                                                                                                                                                                                    >  2:5241/1‑230 (MQ=255)
                                                   gatCCGGCCATACCATTCACAACATGAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCCCTATGAATCTTACCGAATTAAAGAATACGCCGGTCTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATACTGCAGGATGGAtt                                                                                                                                                                                  <  2:154277/251‑1 (MQ=255)
                                                                            tAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCTCAGACAg                                                                                                                                                          <  1:465/250‑1 (MQ=255)
                                                                                            cccAAGTTTAAGAACTCAAACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAACCCTGGCTCGTATGCGTAAGCAGGACATTTTTTTTGCCATTTTGAAGCAGCACGCAAAGCATGGCGAAGATATCTTTGGTGTTGGCGTACTGGAGATATTGCAGGATGGGTTTGGTTTCATCCGTTCAGCTAACAGCTCCTACCTCGCCGGTc                                                                                                                                         >  2:195827/1‑251 (MQ=255)
                                                                                               aaGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTGTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTCAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGAAAAGAGTGGCGAAGAAATCTTTGGTGATGACGTAGTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTACGAAGACAGCTCCTACCACGCCGGTCCtg                                                                                                                                      >  1:69463/1‑251 (MQ=255)
                                                                                                  tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatg                                                                                                                                   >  1:60937/1‑251 (MQ=255)
                                                                                                     aaGAACTCACTCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGAc                                                                                                                                 <  1:146722/250‑1 (MQ=255)
                                                                                                                  acTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTTAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCTGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTcccc                                                                                                                   <  2:48264/251‑1 (MQ=255)
                                                                                                                                cGAATTAAAGAATACGCCGGTTTCTGAGCTGATCCCTCTCGGCGAAATTATGGGGCTGGAAACCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCTATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTTGTGATGGCTTACTGGAGATATTGCAGAATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgcc                                                                                                     <  1:35461/251‑1 (MQ=255)
                                                                                                                                 gAATTAAAGAATTCGCCGTTTTCTGAGCTGATCACGCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATTCGTAAGCAGGACATTATTTTTGCCTTTTTGAAGCAGCACGCAAAGAGTGGCGAATATATCTTGGGTGATGGCTTCCTGGAAATATTGCAGGAGGGATTTGGTTTCCTCCTTTCCGCCGACAGCTCCTCCCTCGCCGTTCCTGATGACATCTACGTTGCCCCTAGCCAAACccgccg                                                                                                    <  2:264058/251‑1 (MQ=255)
                                                                                                                                              cGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATAACTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGGTTTGGTTTCCTCCGTTACGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCAACGTTTTCCCTAGCCAACTCCCCCGTTTCAACCCCCGc                                                                                       >  2:222791/1‑251 (MQ=255)
                                                                                                                                                          aGCTGATCAATCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTCTGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGAGGCGAAGATATCTTTGGTGGTGGAGTACTGGAGATATTGCAGGATGGATTTGGTTTACTCCGTTCCGCAGCCAGCTCCTACCTCGCCGGTCCTGGTGCCCTCTAAGTTTCCCCTAGCCAAATCCGCAGTTTCAACCTCCGCACTGGTGATAcc                                                                           >  2:291656/1‑251 (MQ=255)
                                                                                                                                                                           gAAAATATGGGGCTGGAAAACCTGGATCGTATGCGTAAGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcg                                                          >  1:106287/1‑251 (MQ=255)
                                                                                                                                                                                              accaggCTCATATGCGTAAGAAGGACATTATTTTTGACATTTTGAAGCAGCACACAAAGAGTGGCGAAGATATCTTTGGTGATGGCGAACTGGACATATAGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTCCGTTTACCCTAGCCAAATCCGCAGTTTCAACCTCCGAACTGGTGATACCATCTCTGGAAAGCTTAGCCCGCCGCAAGAAgggga                                        >  2:265940/5‑247 (MQ=255)
                                                                                                                                                                                                              tAAGCAGGAAATTATTTTTGCCATTTTGAAGCAGCACGCCAAGCGTGGCGAAGATTTTTTTGGTGATGGCTTACTGGAGATTTTGCAGGATGGGTTTGGTTTCCTCCTTTCCGCAGACAGCTCCTACCTCGCCGGACCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCATCCTCCGCCCTGGTGATACCATCTCTGGACAGATTCACCCGCCGCAAGAAGGTGAAAGCTATTTTGCgctgt                       >  2:289645/1‑250 (MQ=255)
                                                                                                                                                                                                               aaGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGCAAGAAGGTGAACGCTATTTTGCgctgct                      >  1:219056/1‑251 (MQ=255)
                                                                                                                                                                                                                agccggACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGCCAGCTCCTACCTCGCCGGTCCTGTTGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGa                                                              >  2:260283/5‑210 (MQ=255)
                                                                                                                                                                                                                aGCAGGACATTATTTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGa                                                              <  1:260283/210‑1 (MQ=255)
                                                                                                                                                                                                                           attTTTGCCATTTTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCAGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg           >  1:170775/1‑250 (MQ=255)
                                                                                                                                                                                                                                   cATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGGGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGGTAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt  >  2:148623/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                                                      
GTGACCAGGCAGCGAAAAGACATGAGTCGATGACCGTAAACAGGCATGGATGATCCTGCCATACCATTCACAACATTAAGTTCGAGATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACT  >  NC_000913/3966300‑3966777

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