Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,452,530 G→C 6.2% V1247L (GTA→CTA)  ytfN → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,452,5300GC 6.2% 197.1 / 2.6 64V1247L (GTA→CTA) ytfNconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/26);  new base C (2/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

TATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTT  >  AP009048/4452466‑4452598
                                                                |                                                                    
tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGa                                                                  <  1:228375/69‑1 (MQ=255)
tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGAc                                                                 >  1:381775/1‑70 (MQ=255)
 ataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTa                                                                    >  1:2900083/1‑66 (MQ=255)
 ataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGa                                                                  >  1:364462/1‑68 (MQ=255)
  taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGAcc                                                                <  1:3073247/69‑1 (MQ=255)
  taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGAc                                                                 <  1:3288614/68‑1 (MQ=255)
  taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCa                                                               >  1:2787027/1‑70 (MQ=255)
  taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAg                                                              >  1:2823125/1‑71 (MQ=255)
  taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAg                                                              >  1:4225690/1‑71 (MQ=255)
  taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAg                                                              >  1:4257586/1‑71 (MQ=255)
   aGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAg                                                                   <  1:1962120/65‑1 (MQ=255)
     cAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                           <  1:4270169/71‑1 (MQ=38)
      aaCACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAg                                                                   <  1:3872444/62‑1 (MQ=255)
       acacTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACt                                                         <  1:505772/71‑1 (MQ=255)
        cacTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTg                                                        >  1:3993306/1‑71 (MQ=255)
            cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAtt                                                    >  1:1263025/1‑71 (MQ=255)
            cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAtt                                                    >  1:3574659/1‑71 (MQ=255)
             acgttacgttATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttt                                                   >  1:3944379/1‑71 (MQ=255)
             acgttacgttATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAtt                                                    <  1:1888161/70‑1 (MQ=255)
             acgttacgttATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAt                                                     <  1:4616915/69‑1 (MQ=255)
              cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTa                                                                    >  1:7862/1‑53 (MQ=255)
               gttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGc                                                 >  1:3166020/1‑71 (MQ=255)
               gttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGc                                                 <  1:4764269/71‑1 (MQ=255)
                ttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGct                                                >  1:3028840/1‑71 (MQ=255)
                 tacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGct                                                >  1:831210/1‑70 (MQ=255)
                    gttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTAt                                            >  1:2184392/1‑71 (MQ=255)
                     ttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATc                                           <  1:1359347/71‑1 (MQ=255)
                       aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttt                                                   <  1:719870/61‑1 (MQ=255)
                       aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAg                                         >  1:3889355/1‑71 (MQ=255)
                        tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTa                                             <  1:4493873/66‑1 (MQ=255)
                        tCGCCTGATGCCTAAGCTATATCTGGAAGACGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGt                                        >  1:2799527/1‑71 (MQ=255)
                         cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTa                                             <  1:2339994/65‑1 (MQ=255)
                          gCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGtt                                       <  1:428932/70‑1 (MQ=255)
                           ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTc                                      >  1:3971944/1‑70 (MQ=255)
                            cTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttt                                                   >  1:67652/1‑56 (MQ=255)
                             tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAg                                   >  1:5013915/1‑71 (MQ=255)
                             tGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTAGACCAg                                                              <  1:2111638/44‑1 (MQ=38)
                                 gCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACt                                                         <  1:2016711/45‑1 (MQ=255)
                                   cTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTa                              >  1:1774430/1‑70 (MQ=255)
                                      aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGtt                                       >  1:3343744/1‑58 (MQ=255)
                                      aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCa                           <  1:4991903/70‑1 (MQ=255)
                                      aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCa                           <  1:1860162/70‑1 (MQ=255)
                                       gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCaa                          >  1:3087268/1‑70 (MQ=255)
                                       gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAAt                         <  1:2825162/71‑1 (MQ=255)
                                        cTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                        <  1:3102649/71‑1 (MQ=255)
                                        cTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                        <  1:1739029/71‑1 (MQ=255)
                                             tCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAg                             >  1:2793517/1‑61 (MQ=255)
                                                ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGAttt                                                   >  1:2466951/1‑36 (MQ=255)
                                                 gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAata                  <  1:3083517/68‑1 (MQ=255)
                                                 gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAata                  <  1:321085/68‑1 (MQ=255)
                                                 gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                <  1:2310630/70‑1 (MQ=255)
                                                    gCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGa                                    <  1:212205/47‑1 (MQ=255)
                                                    gCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAg                                   <  1:4549517/48‑1 (MQ=255)
                                                    gCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGc                            <  1:2527121/55‑1 (MQ=255)
                                                      cGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa          >  1:1155360/1‑71 (MQ=255)
                                                         gtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                >  1:2827487/1‑62 (MQ=255)
                                                         gtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACgg       >  1:4185317/1‑71 (MQ=255)
                                                         gtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACgg       >  1:1422961/1‑71 (MQ=255)
                                                          tCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACgg       >  1:811605/1‑70 (MQ=255)
                                                          tCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACgg       >  1:4117328/1‑70 (MQ=255)
                                                          tCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACgg       >  1:2954071/1‑70 (MQ=255)
                                                            tGGTCTACACCAGACACTGGATTTGCTCTATCAGTTCGAGTTTTAGCaa                          >  1:104698/1‑49 (MQ=37)
                                                             ggTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGt   >  1:4310403/1‑71 (MQ=255)
                                                             ggTCTACACCAGACACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACg        >  1:1200107/1‑66 (MQ=255)
                                                              gtgtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCa     <  1:1779184/68‑1 (MQ=255)
                                                              gtgtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGtt  <  1:5164024/71‑1 (MQ=255)
                                                              gtgtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGtt  >  1:3893399/1‑71 (MQ=255)
                                                              gtgtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGtt  >  1:4263916/1‑71 (MQ=255)
                                                               tgtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCaa                          <  1:2174361/46‑1 (MQ=255)
                                                                gtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGtt  <  1:5179845/69‑1 (MQ=255)
                                                                |                                                                    
TATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTT  >  AP009048/4452466‑4452598

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: