Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,559,548 T→C 60.7% I293T (ATT→ACT)  hemF → coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,559,5480TC60.7% 18.0 / 28.7 28I293T (ATT→ACT) hemFcoproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base T (5/6);  new base C (10/7);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTACGCTGGGAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTATTAAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCATCCGCTGATGCAGCGTCAGTGACG  >  W3110S.gb/2559483‑2559616
                                                                 |                                                                    
ggTACGCTGGGAATATGATTATCAGCCAAAAGATGGCAGCCAAGAAGCGGCGTTAAGTGAGTTTACTAAgg                                                                 >  1:886219/1‑71 (MQ=255)
         ggAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTACTAAGGTCAGGGa                                                          <  1:992557/69‑1 (MQ=255)
         ggAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTATTAAGGTCAGGGATTg                                                       <  1:2378544/71‑1 (MQ=255)
          gAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTATTAAGGTCAGGGATTg                                                       <  1:757295/71‑1 (MQ=255)
             tatGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTATTAAGGTCAgg                                                            <  1:1563669/63‑1 (MQ=255)
             tatGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTACTAAGGTCAGGGATTgg                                                      <  1:2764634/69‑1 (MQ=255)
               tGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTATTAAGGTCAGGGAtt                                                        >  1:679658/1‑65 (MQ=255)
               tGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTACTAAGGTCAGGGATTGGgtg                                                   >  1:727804/1‑70 (MQ=255)
               tGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTACTAAGGTCAGGGATTGGgtg                                                   >  1:111099/1‑70 (MQ=255)
                             aaGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTATTAAGGTCAGGGATTGGGTGTAACTCCCTCAcccc                                    <  1:2401909/71‑1 (MQ=255)
                              aGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTATTAAGGTCAGGGa                                                          >  1:2536020/1‑48 (MQ=255)
                                aTGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTACTAAGGTCAGGGATTGGGTGTaa                                                >  1:2791190/1‑56 (MQ=255)
                                aTGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTACTAAGGTCAGGGATTGGGTGTaa                                                >  1:1118333/1‑56 (MQ=255)
                                  ggCAGCCCAGAAGCGGCGTTAAGTGAGTTTATTAAGGTCAGGGATTGGGTGTAACTcc                                            >  1:464089/1‑58 (MQ=255)
                                       cccAGAAGCGGCGTTAAGTGAGTTTACTAAGGTCAGGGATTGGGTGTAACTCCCt                                          <  1:1275966/55‑1 (MQ=255)
                                       cccAGAAGCGGCGTTAAGTGAGTTTACTAAGGTCAGGGATTGGGTGTAACTCCCt                                          <  1:2759680/55‑1 (MQ=255)
                                        ccAGAAGCGGCGTTAAGTGAGTTTATTAAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCAt                         >  1:1081618/1‑71 (MQ=255)
                                         cAGAAGCGGCGTTAAGTGAGTTTACTAAGGTCAGGGATTGGGTTTAACTcc                                            <  1:1608479/51‑1 (MQ=255)
                                         cAGAAGCGGCGTTAAGTGAGTTTACTAAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCATc                        >  1:31551/1‑71 (MQ=255)
                                         cAGAAGCGGCGTTAAGTGAGTTTACTAAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCATc                        >  1:1048059/1‑71 (MQ=255)
                                              gcggcgTTAAGTGAGTTTACTAAGGTCAGGGATTGGGTGTACCTCCCTCACCCCCACTCCCGCATc                        >  1:2735917/1‑66 (MQ=255)
                                                      aaGTGAGTTTATTAAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCATCCGCTGATGCAGCg           >  1:2694806/1‑71 (MQ=255)
                                                        gTGAGTTTATTAAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCATCCGCTGATGCAGCGTc         <  1:286860/71‑1 (MQ=255)
                                                        gTGAGTTTATTAAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCATCCGCTGATGCAGCGTc         <  1:1864856/71‑1 (MQ=255)
                                                          gAGTTTACTAAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCATCCGCTGATGCAGCGTCAg       >  1:2191299/1‑71 (MQ=255)
                                                          gAGTTTACTAAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCATCCGCTGATGCAGCGTCAg       >  1:2183299/1‑71 (MQ=255)
                                                            gTTTACTAAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCa                          <  1:259615/50‑1 (MQ=255)
                                                            gTTTACTAAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCa                          <  1:1646054/50‑1 (MQ=255)
                                                               tattaAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCATCCGCTGATGCAGCGTCAGTGACg  >  1:790211/1‑71 (MQ=255)
                                                               tattaAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCATCCGCTGATGCAGCGTCAGTGACg  >  1:1376954/1‑71 (MQ=255)
                                                                 |                                                                    
GGTACGCTGGGAATATGATTATCAGCCAAAAGATGGCAGCCCAGAAGCGGCGTTAAGTGAGTTTATTAAGGTCAGGGATTGGGTGTAACTCCCTCACCCCCACTCCCGCATCCGCTGATGCAGCGTCAGTGACG  >  W3110S.gb/2559483‑2559616

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