Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,408,005 T→C 8.3% F105L (TTT→CTT)  ydaL → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,408,0050TC 8.3% 122.2 / 4.0 48F105L (TTT→CTT) ydaLconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (13/31);  new base C (2/2);  total (15/33)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATA  >  W3110S.gb/1407940‑1408074
                                                                 |                                                                     
cGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTAtt                                                                  >  1:1454287/1‑71 (MQ=255)
 gCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTc                                                                 <  1:2935673/71‑1 (MQ=255)
  cAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTc                                                                 >  1:1053553/1‑70 (MQ=255)
   aaCAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTcaa                                                               <  1:2944254/71‑1 (MQ=255)
   aaCAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTcaa                                                               <  1:4290942/71‑1 (MQ=255)
       ggCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAg                                                           >  1:348113/1‑71 (MQ=255)
                gAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGcc                                                         >  1:4001936/1‑64 (MQ=255)
                   tCTTTTGCTCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTCTTATTcaa                                                               >  1:4133566/1‑55 (MQ=255)
                   tCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAg                                                           <  1:3295257/59‑1 (MQ=255)
                   tCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTCTTATTcaa                                                               >  1:3330721/1‑55 (MQ=255)
                    cTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGGTCAGTTTTATTcaacaa                                                            <  1:4085376/57‑1 (MQ=255)
                     ttttGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGt                                             >  1:1739162/1‑71 (MQ=255)
                      tttGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGtt                                            >  1:1819581/1‑71 (MQ=255)
                        tGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTg                                          >  1:614119/1‑71 (MQ=255)
                        tGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTg                                          >  1:1447789/1‑71 (MQ=255)
                          cgcCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGc                                         <  1:86715/70‑1 (MQ=255)
                            ccAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTa                                      <  1:463659/71‑1 (MQ=255)
                             cAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTaa                                     >  1:1798448/1‑71 (MQ=255)
                              aGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAAc                                    <  1:4129496/71‑1 (MQ=255)
                               gCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACg                                   <  1:682860/71‑1 (MQ=255)
                                   gTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGa                                                 <  1:3616289/53‑1 (MQ=255)
                                    tGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTcaac                                                              <  1:172836/39‑1 (MQ=255)
                                    tGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTcaac                                                              <  1:1331632/39‑1 (MQ=255)
                                      gAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTcaa                                                               <  1:2718810/36‑1 (MQ=255)
                                      gAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATgg                                              <  1:68127/53‑1 (MQ=255)
                                        agagTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGatta                          <  1:960043/71‑1 (MQ=255)
                                        agagTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGatta                          <  1:1342313/71‑1 (MQ=255)
                                        agagTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGatta                          <  1:1461787/71‑1 (MQ=255)
                                         gagTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGattat                         <  1:463844/71‑1 (MQ=255)
                                         gagTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGattat                         <  1:1874373/71‑1 (MQ=255)
                                           gTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGattatt                        >  1:2460664/1‑70 (MQ=255)
                                            tGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATgg                                              <  1:433625/47‑1 (MQ=255)
                                              ccGCAAAATGGTGTTCAGTTTTATTAAACAAGCCCTGGCGGATgg                                              >  1:3240487/1‑45 (MQ=255)
                                                  aaaaTGGTGTTCAGTTTTATTCAACAAGCCATGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTaa                <  1:1250551/71‑1 (MQ=255)
                                                   aaaTGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACg                                   >  1:3587815/1‑51 (MQ=255)
                                                    aaTGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGt                                       <  1:2574830/46‑1 (MQ=255)
                                                     aTGGTGTTCAGTTTTAATCAACAAGCCCTGGCGGATgg                                              >  1:3795558/1‑38 (MQ=255)
                                                           ttCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCg                                        <  1:3989242/38‑1 (MQ=255)
                                                           ttCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGat                            <  1:1824505/50‑1 (MQ=255)
                                                           ttCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATTGTAAAGGGCggg        <  1:4439761/70‑1 (MQ=255)
                                                           ttCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTaa                <  1:4358003/62‑1 (MQ=255)
                                                           ttCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCggg        <  1:4541267/70‑1 (MQ=255)
                                                           ttCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCggg        <  1:670525/70‑1 (MQ=255)
                                                           ttCAGTCTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCAt                     <  1:801430/57‑1 (MQ=255)
                                                           ttCAGTCTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCAt                     <  1:4227287/57‑1 (MQ=255)
                                                               gTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGat                            <  1:2957803/46‑1 (MQ=255)
                                                               gTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGat                            <  1:1398585/46‑1 (MQ=255)
                                                               gTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGat                            <  1:924824/46‑1 (MQ=255)
                                                                ttttATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGt                  <  1:1384718/55‑1 (MQ=255)
                                                                ttttATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTa                 <  1:687813/56‑1 (MQ=255)
                                                                ttttATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATa  <  1:1518974/71‑1 (MQ=255)
                                                                ttttATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATa  <  1:3285293/71‑1 (MQ=255)
                                                                 |                                                                     
CGCAACAGGCGAGCCTGAATCTTTTGCGCCAGCCGGTGGAAGAGTGCCGCAAAATGGTGTTCAGTTTTATTCAACAAGCCCTGGCGGATGGTTTGCGTAACGTGCTGATTATTCATGGTAAAGGGCGGGATGATA  >  W3110S.gb/1407940‑1408074

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