Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,108,438 C→A noncoding (56/76 nt) pheV → tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,108,4380CA96.7% 171.2 / NA 60noncoding (56/76 nt)pheVtRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (14/44);  total (14/46)
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.99e-01

GTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGA  >  NC_000913/3108366‑3108508
                                                                        |                                                                      
gTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGAt                                                                     <  2:3045629/76‑1 (MQ=0)
                  cccGGATAGCTCAGTCGGTAGAACAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGcaccac                                                   <  1:1268418/76‑1 (MQ=11)
                  cccGGATAGCGCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGcaccac                                                   <  2:386778/76‑1 (MQ=255)
                   ccGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACt                                                  >  1:1113908/1‑76 (MQ=12)
                   ccGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACt                                                  <  1:1360850/76‑1 (MQ=12)
                    cGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGCGCCCGGGCACCACTa                                                 >  2:263273/1‑76 (MQ=11)
                    cGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTa                                                 <  1:1860813/76‑1 (MQ=19)
                     ggATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGCGTCCGGGCACCACTaa                                                >  1:4479966/1‑76 (MQ=255)
                     ggATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTaa                                                >  1:539428/1‑76 (MQ=255)
                     ggATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTaa                                                >  1:2136275/1‑76 (MQ=255)
                        tAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTc                                             >  1:4315634/1‑76 (MQ=255)
                          gCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCtt                                           >  2:3397780/1‑76 (MQ=255)
                          gCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCtt                                           <  1:3901221/76‑1 (MQ=255)
                            tCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTaa                                         >  1:1000647/1‑76 (MQ=255)
                              aGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGa                                       >  2:2074975/1‑76 (MQ=255)
                               gTCGGTAGAGGAGGGGATTGAAAATCCCCGTGTCCGTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGaa                                      <  2:1218626/76‑1 (MQ=17)
                               gTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGaa                                      <  2:1188927/76‑1 (MQ=255)
                                gctgTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAc                                     <  1:361696/73‑1 (MQ=255)
                                 cGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAcc                                    <  2:2925691/76‑1 (MQ=255)
                                  ggTAGAGGAGGGTGTTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAGccc                                   <  1:1729022/76‑1 (MQ=17)
                                  ggTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAccc                                   <  2:4356205/76‑1 (MQ=255)
                                      gagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGccc                               <  1:2563704/76‑1 (MQ=255)
                                      gagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGccc                               <  1:4386996/76‑1 (MQ=255)
                                      gagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGccc                               <  1:3647005/76‑1 (MQ=255)
                                       agCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCca                              <  1:1531161/76‑1 (MQ=255)
                                       agCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCca                              <  1:4115134/76‑1 (MQ=255)
                                        gCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGCGTCCGGGCACCACTAATTCTTAAGAACCCTCCcac                             <  2:2868615/76‑1 (MQ=255)
                                         cAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcaca                            <  1:2738520/76‑1 (MQ=255)
                                           ggggATTGAATATCCCCGTGTCCTTGTTTAGATTCCTAGTCCGGGCCCCACTAATTCTTAAGAACCCTCCCTCaat                          <  2:3487303/76‑2 (MQ=18)
                                           ggggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAg                          <  2:3689579/76‑1 (MQ=255)
                                            gggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAgg                         <  1:919846/76‑1 (MQ=255)
                                             ggATTGACAACCTCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGACCCCGCCCCCAAGGc                        <  2:2878791/76‑1 (MQ=18)
                                             ggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCGACAAGGc                        <  1:2661288/76‑1 (MQ=21)
                                             ggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGc                        >  2:2430240/1‑76 (MQ=255)
                                                ttGAGAATCCCCGTGTCCTTGGTTAGATGCCGAGTCCGGGCACCACGAATTCTTAAGACCCCGCCCACAAGGCggg                     <  1:4772355/76‑1 (MQ=21)
                                                ttGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCggg                     <  2:3757777/76‑1 (MQ=255)
                                                ttGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCggg                     <  1:4766160/76‑1 (MQ=255)
                                                 tGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                    <  2:639992/76‑1 (MQ=255)
                                                 tGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                    <  1:60856/76‑1 (MQ=255)
                                                  gAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtt                   <  2:2830036/76‑1 (MQ=255)
                                                  gAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtt                   <  1:4700063/76‑1 (MQ=255)
                                                  gAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtt                   <  2:1616921/76‑1 (MQ=255)
                                                  gAAAAACCCCGTGGCCTTGGTTAGACTCAGAGACCGCGCACCACTACTTCTTAAGAACCCGCCCACAAGTCgggct                   >  2:3514860/1‑74 (MQ=255)
                                                    aaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtttt                 <  2:4175576/76‑1 (MQ=255)
                                                     aaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttttt                <  2:406388/76‑1 (MQ=255)
                                                     aaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttttt                <  1:4797046/76‑1 (MQ=255)
                                                      aTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTg               <  1:1450356/76‑1 (MQ=255)
                                                      aTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTg               <  1:4131183/76‑1 (MQ=255)
                                                      aTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCGTAAGAACCCGCCCACAAGGCGGGTTTTTg               <  1:4132715/76‑1 (MQ=25)
                                                       tCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGc              >  2:2566019/1‑76 (MQ=255)
                                                            gtgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTg         <  2:1426047/76‑1 (MQ=255)
                                                             tgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTgg        <  2:1088214/76‑1 (MQ=255)
                                                              gtCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGCACCCGCCCACAAGGCGGGTTTTTGCTTTTGGa       <  1:3570467/76‑1 (MQ=255)
                                                              gtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGa       <  2:4294992/76‑1 (MQ=255)
                                                                ccTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc     <  1:4249693/76‑1 (MQ=255)
                                                                ccTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc     >  2:3251138/1‑76 (MQ=255)
                                                                 cTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt    >  1:3678427/1‑76 (MQ=255)
                                                                 cTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt    <  2:4308006/76‑1 (MQ=255)
                                                                   tGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGa  <  1:960208/76‑1 (MQ=255)
                                                                   tGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGa  <  1:3765642/76‑1 (MQ=255)
                                                                        |                                                                      
GTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGA  >  NC_000913/3108366‑3108508

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