Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,702,790 C→T 69.4% intergenic (‑71/+201) rnc ← / ← lepB RNase III/leader peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,702,7900CT69.4% 62.2 / 40.8 49intergenic (‑71/+201)rnc/lepBRNase III/leader peptidase
Reads supporting (aligned to +/- strand):  ref base C (11/4);  new base T (18/16);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

AATAAATGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCG  >  W3110S.gb/2702732‑2702848
                                                          |                                                          
aataaatGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGc                                                    >  1:1585635/1‑67 (MQ=255)
aataaatGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGc                                                    >  1:1201839/1‑67 (MQ=255)
  taaatGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGca                                                <  1:3468625/69‑1 (MQ=255)
  taaatGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGca                                                <  1:2486544/69‑1 (MQ=255)
  taaatGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGca                                                <  1:1629711/69‑1 (MQ=255)
  taaatGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGc                                                 <  1:2608985/68‑1 (MQ=255)
          tCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGcagtttcag                                         >  1:780446/1‑68 (MQ=255)
          tCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGcagtttcag                                         >  1:1352516/1‑68 (MQ=255)
           caacaaTTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGc                                                    <  1:1897398/56‑1 (MQ=255)
              caaTTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAAc                                                            >  1:3546874/1‑45 (MQ=255)
              caaTTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGcagtttcagtt                                       >  1:433356/1‑66 (MQ=255)
                  tCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGc                                                 >  1:930674/1‑52 (MQ=255)
                  tCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCgt                                 >  1:139341/1‑67 (MQ=255)
                  tCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCgg                                 >  1:1488123/1‑68 (MQ=255)
                  tCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCgg                                 >  1:652898/1‑68 (MQ=255)
                  tCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCgg                                 >  1:1786766/1‑68 (MQ=255)
                   cAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGta                               <  1:2660670/69‑1 (MQ=255)
                   cAGCACACGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGta                               <  1:242368/69‑1 (MQ=255)
                   cAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGta                               <  1:2469854/69‑1 (MQ=255)
                      cacaCGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTatat                            <  1:2982984/69‑1 (MQ=255)
                      cacaCGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTatat                            <  1:826168/69‑1 (MQ=255)
                       acacGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGca                                                <  1:1358125/48‑1 (MQ=255)
                       acacGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGTATATg                           <  1:3618897/69‑1 (MQ=255)
                          cGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGtata                             <  1:3661046/64‑1 (MQ=255)
                          cGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGtata                             >  1:510939/1‑64 (MQ=255)
                          cGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGTATATgggg                        >  1:417149/1‑69 (MQ=255)
                          cGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGTATATgggg                        >  1:536074/1‑69 (MQ=255)
                          cGAAACAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGTATATgggg                        >  1:870830/1‑69 (MQ=255)
                            aaaCAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGtat                              >  1:2272218/1‑61 (MQ=255)
                            aaaCAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGtat                              >  1:1268253/1‑61 (MQ=255)
                            aaaCAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGtat                              >  1:754321/1‑61 (MQ=255)
                            aaaCAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGtat                              >  1:2118265/1‑61 (MQ=255)
                            aaaCAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGtat                              >  1:265899/1‑61 (MQ=255)
                            aaaCAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGtat                              >  1:2648249/1‑61 (MQ=255)
                             aaCAGACCTGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGAcc                     <  1:3522772/69‑1 (MQ=255)
                                     tGATATACATGTTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGt             <  1:1971497/69‑1 (MQ=255)
                                     tGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGt             <  1:2849099/69‑1 (MQ=255)
                                     tGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGt             <  1:1628565/69‑1 (MQ=255)
                                     tGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGt             <  1:2172924/69‑1 (MQ=255)
                                     tGATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGt             <  1:2376004/69‑1 (MQ=255)
                                      gATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGta                               >  1:2950467/1‑50 (MQ=255)
                                      gATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGta                               >  1:2702071/1‑50 (MQ=255)
                                      gATATACATGGTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGta                               >  1:2544267/1‑50 (MQ=255)
                                         ataCATGGTTCTGCTAATTGCTTCGCTGcagtttcagttt                                      >  1:556739/1‑40 (MQ=255)
                                         ataCATGGTTCTGCTAATTGCTTCGCTGcagtttcagttt                                      >  1:1843150/1‑40 (MQ=255)
                                         ataCATGGTTCTGCTAATTGCTTCGCTGcagtttcagttt                                      >  1:1942816/1‑40 (MQ=255)
                                            cATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCgtgt      >  1:1057757/1‑69 (MQ=255)
                                                gTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGta                               >  1:3209444/1‑40 (MQ=255)
                                                gTTCTGCTAATTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCg  <  1:2046332/69‑1 (MQ=255)
                                                  tCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGta                               >  1:1655876/1‑38 (MQ=255)
                                                          |                                                          
AATAAATGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCG  >  W3110S.gb/2702732‑2702848

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