Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,766,021	G→A	W1280* (TGG→TAG)	rhsA →	rhs element protein RhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,766,021	0	G	A	98.0%	169.0 / NA	50	W1280* (TGG→TAG)	rhsA	rhs element protein RhsA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/23);  minor base T (1/0);  total (27/23)
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

ATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTG  >  NC_000913/3765892‑3766146                                                                                                                                  |                                                                                                                              aTCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAAtagtg                                                                                                                            <  1:505054/133‑1 (MQ=255) aTCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCACCAAGCAGAAGATGATGTTGCAAAAGGATTGACAAAAtagtg                                                                                                                            >  2:141959/1‑133 (MQ=255)  tCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAAtagtgg                                                                                                                           <  2:396401/133‑1 (MQ=255)     aGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAAc                                                                                                                        >  2:487525/1‑133 (MQ=255)     aGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAAc                                                                                                                        >  2:396098/1‑133 (MQ=255)     aGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAAc                                                                                                                        <  1:265076/133‑1 (MQ=255)                 aGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATTATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATca                                                                                                            >  2:280537/1‑133 (MQ=255)                   aTCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATa                                                                                                                    <  2:89971/123‑1 (MQ=255)                   aTCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATa                                                                                                                    >  1:89971/1‑123 (MQ=255)                     ccTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACaaa                                                                                                        <  1:93454/133‑1 (MQ=255)                     ccTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACaaa                                                                                                        >  2:93454/1‑133 (MQ=255)                      ctctGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAg                                                                                                       >  1:41474/1‑133 (MQ=255)                      ctctGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTCACAAAATTGTGGAACGCTACAGCATCACAAAg                                                                                                       >  2:185579/1‑133 (MQ=255)                       tctGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAGAAGCATTGACAAAACATTGGAACGAAACAGCATCACAAAGa                                                                                                      >  2:226782/1‑133 (MQ=255)                       tctGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGa                                                                                                      >  1:477794/1‑133 (MQ=255)                       tctGGGGTTAGAAGTGTTTCCTGGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGa                                                                                                      <  1:168319/133‑1 (MQ=255)                           gggTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAatat                                                                                                  >  2:59683/1‑133 (MQ=255)                               tAGAAGTTTTTCCTATACGATTCCCCTTGCAAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGa                                                                                              <  1:167091/133‑1 (MQ=255)                               tAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATTATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGa                                                                                              <  1:280537/133‑1 (MQ=255)                               tAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGa                                                                                              <  2:386090/133‑1 (MQ=255)                               tAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGa                                                                                              <  2:427036/133‑1 (MQ=255)                               tAGAAGTTTTCCCTATACCATTCCTCTTCCCATTTCCATGGACCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGa                                                                                              <  2:132015/133‑1 (MQ=255)                               tAGAAGTTCTTCCTACACCATTACCCTTGCTAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGa                                                                                              <  1:141959/133‑1 (MQ=255)                                aGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTTAc                                                                                             >  1:402203/1‑133 (MQ=255)                                    gTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTTActctc                                                                                         <  2:402203/133‑1 (MQ=255)                                    gTTTTTCCGAGACCATTCCCCTGGCCAATTGCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGActctc                                                                                         <  1:175486/133‑1 (MQ=255)                                              gACCATTCCCCTTGCCAATTCCCTGGCCCATAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGaa                                                                               <  1:185579/133‑1 (MQ=255)                                              gACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGaa                                                                               <  2:27043/133‑1 (MQ=255)                                                   ttCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAACGAATATTTGACCCTATAATATTGAATAATc                                                                          >  1:175565/1‑133 (MQ=255)                                                         ttGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGAc                                                                    <  2:477794/133‑1 (MQ=255)                                                         ttGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGAc                                                                    <  2:477903/133‑1 (MQ=255)                                                          tGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACt                                                                   >  1:91775/1‑133 (MQ=255)                                                          tGCCAATACCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACt                                                                   <  2:452770/133‑1 (MQ=255)                                                           gCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTa                                                                  >  1:405178/1‑133 (MQ=255)                                                               aTTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCAt                                                              >  2:364974/1‑133 (MQ=255)                                                                    aTGGCTCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAATGTATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGgg                                                         >  1:320121‑M2/1‑131 (MQ=255)                                                                    aTGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGgg                                                         >  2:41474‑M2/1‑131 (MQ=255)                                                                         ccAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggt                                                    >  1:396098‑M2/1‑126 (MQ=255)                                                                         ccAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggt                                                    >  1:487525‑M2/1‑126 (MQ=255)                                                                          cAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtg                                                   <  1:261581‑M2/133‑9 (MQ=255)                                                                            aaaGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcg                                                 <  1:260304‑M2/133‑11 (MQ=255)                                                                               gCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaac                                              >  1:394555‑M2/1‑120 (MQ=255)                                                                                   tGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagt                                          <  2:10955‑M2/133‑18 (MQ=255)                                                                                     cacaGCAGGAAGCAGGTGGTAACGCTGCAATAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtcc                                        >  1:226782‑M2/1‑114 (MQ=255)                                                                                          cagcaAGCAGATGATAATGCTGCAAAAGCATTTACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgat                                   <  2:356197‑M2/133‑25 (MQ=255)                                                                                          cagcaAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgat                                   >  1:84109‑M2/1‑109 (MQ=255)                                                                                          cagcaAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgat                                   <  2:84109‑M2/133‑25 (MQ=255)                                                                                          cagcaAGCAGATGATAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgat                                   >  1:10955‑M2/1‑109 (MQ=255)                                                                                                   gatgatAATGCTGCAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatca                          >  2:261581‑M2/1‑100 (MQ=255)                                                                                                              tgcAAAAGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatc               <  2:54388‑M2/133‑45 (MQ=255)                                                                                                                 aaaaGCATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgg            >  1:437525‑M2/1‑86 (MQ=255)                                                                                                                      cATTGACAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagcca       >  1:451286‑M2/1‑81 (MQ=255)                                                                                                                           aCAAAATAGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaa  >  2:260304‑M2/1‑76 (MQ=255)                                                                                                                                  |                                                                                                                              ATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGTGGAAATGTTGCAGACACAGGGATAACTGATCGTG  >  NC_000913/3765892‑3766146

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