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.2% 161.7 / ‑5.7 54noncoding (56/76 nt)pheVtRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base A (22/29);  minor base T (1/1);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGAC  >  NC_000913/3108383‑3108509
                                                       |                                                                       
gcgcGGATAGCTCATTCGGTATGGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGCGTCCGGGctcca                                                     <  1:2831616/73‑1 (MQ=1)
 cccGGGTAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGcaccac                                                    <  2:231875/76‑1 (MQ=255)
    ggATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTTGATTCCGAGTCCGGGCACCCCTaa                                                 >  2:2926506/1‑76 (MQ=255)
       tAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTc                                              <  2:955415/76‑1 (MQ=255)
         gCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCtt                                            <  2:2616633/76‑1 (MQ=255)
          cTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTa                                           >  2:4780571/1‑76 (MQ=255)
           tCAGTCGGTAGAGCAGGGTATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTaa                                          <  1:2646897/76‑1 (MQ=255)
            cAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAg                                         >  1:1478297/1‑76 (MQ=255)
            cAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAg                                         >  2:2968955/1‑76 (MQ=255)
             aGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGa                                        <  1:858731/76‑1 (MQ=255)
              gTCGGTAGAGCAGGGGATTGACAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGaa                                       >  2:3488217/1‑76 (MQ=17)
              gTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGaa                                       >  1:4178026/1‑76 (MQ=255)
              gTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGaa                                       >  2:2528182/1‑76 (MQ=255)
                cGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAcc                                     >  2:567151/1‑76 (MQ=255)
                 ggcagacCAGGCGACTGAAAATCCCCGTGTTCGTGCCGAGATTCCAAGGCCGCACCCCACTAAAACTTAAGAAccc                                    >  2:3007932/4‑76 (MQ=255)
                 ggTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAccc                                    >  2:3001740/1‑76 (MQ=255)
                  gTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCTGGGCACCACTAATTCTTAAGAACCCg                                   <  1:464535/76‑1 (MQ=18)
                    agagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGcc                                 <  1:599580/76‑1 (MQ=255)
                    agagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGcc                                 <  2:2318019/76‑1 (MQ=255)
                     gagCGGGGGATTGAATATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCCCTAATTCTTCAGACCCCGcgc                                <  2:1105205/76‑3 (MQ=14)
                       gCAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcac                              >  2:1745688/1‑76 (MQ=255)
                        cAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcaca                             <  1:2062004/76‑1 (MQ=255)
                        cAGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcaca                             <  2:4285717/76‑1 (MQ=255)
                         cggggATTTAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                            <  1:5031383/75‑1 (MQ=21)
                         aGGGGATTGAAAGTCCGCGTGTCCTTGGTTTGATTCCGAGTCCGGGCATCCTTATTTCTTAAGAACCCGCCCACaa                            <  2:1857044/76‑1 (MQ=17)
                         aGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                            >  1:2804620/1‑76 (MQ=255)
                         aGGGGATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACaa                            <  2:1376312/76‑1 (MQ=255)
                            ggATTGAGAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGc                         <  2:2310330/76‑1 (MQ=255)
                            ggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGc                         >  2:1909714/1‑76 (MQ=255)
                            ggATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGc                         <  1:3961234/76‑1 (MQ=255)
                             gATTGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCg                        >  2:4014141/1‑76 (MQ=255)
                               ttGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCggg                      <  1:4398483/76‑1 (MQ=255)
                                tGAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                     <  1:202323/76‑1 (MQ=255)
                                 gAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtt                    <  2:127422/76‑1 (MQ=255)
                                 gAAAATCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtt                    >  2:4171592/1‑76 (MQ=255)
                                   aaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGtttt                  >  1:2084154/1‑76 (MQ=255)
                                   aaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                     >  2:2031693/1‑73 (MQ=255)
                                   aaaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                     <  1:2031693/73‑1 (MQ=255)
                                    aaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttttt                 <  2:4178025/76‑1 (MQ=255)
                                    aaTCCCCGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGctttt                 >  2:1162391/1‑76 (MQ=255)
                                        cccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGcc                                 <  1:3895050/56‑1 (MQ=255)
                                        cccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGcc                                 >  2:3895049/1‑56 (MQ=255)
                                        cccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtt             >  2:3333070/1‑76 (MQ=255)
                                        cccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCtt             <  2:4513679/76‑1 (MQ=255)
                                         ccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCttt            <  2:1553232/76‑1 (MQ=255)
                                         ccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCttt            <  1:1396564/76‑1 (MQ=255)
                                         ccGTGTCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGCTTTTGCttt            >  1:163650/1‑76 (MQ=255)
                                          cGTGTCCTTGGTTAGATTCCGAGTCGGGGCACCACTAATTCTTAAGACCCCGCCCACAAGGCGGGTTTTTGCtttt           <  1:1215642/76‑1 (MQ=38)
                                            tgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTgg         >  2:392884/1‑76 (MQ=255)
                                            tgtCCTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTgg         >  2:2406007/1‑76 (MQ=255)
                                                cTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGGACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt     <  2:2887972/76‑1 (MQ=38)
                                                cTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt     <  2:3440484/76‑1 (MQ=255)
                                                cTTGGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt     <  1:148529/76‑1 (MQ=255)
                                                   ggTTAGATTCCGAGTCCGGGCACCACTAAGTCTGAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGAc  <  2:2093440/76‑1 (MQ=39)
                                                       |                                                                       
GCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGAC  >  NC_000913/3108383‑3108509

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