Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,606,094 (G)8→7 noncoding (72/87 nt) leuV ← tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,0870G.96.5% 110.2 / ‑2.0 29noncoding (79/87 nt)leuVtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (12/16);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

TTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGA  >  NC_000913/4605957‑4606221
                                                                                                                                  |                                                                                                                                      
ttGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAAAGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggac                                                                                                                                >  2:633880/1‑137 (MQ=255)
       aTTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAAcc                                                                                                                        <  2:479536/139‑1 (MQ=255)
          aCATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAAccccc                                                                                                                     <  1:541140/139‑1 (MQ=255)
             ttttGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                  <  1:756000/139‑1 (MQ=255)
             ttttGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                  <  1:730158/139‑1 (MQ=255)
               ttGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTc                                                                                                                <  1:633880/139‑1 (MQ=255)
                tGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTc                                                                                                                >  1:16994/1‑139 (MQ=255)
                   tACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                            <  2:686532/139‑1 (MQ=255)
                            gCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTa                                                                                                   >  1:234615/1‑139 (MQ=255)
                                ggCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAacac                                                                                               >  1:569529/1‑139 (MQ=255)
                                 gCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACAcc                                                                                              >  1:1033713/1‑139 (MQ=255)
                                       aTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                                  <  1:611912/129‑1 (MQ=255)
                                       aTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                                  >  2:611912/1‑129 (MQ=255)
                                                              aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                  <  1:673402/90‑1 (MQ=255)
                                                              aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                  >  2:673402/1‑90 (MQ=255)
                                                              aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                                 >  1:805980/1‑139 (MQ=255)
                                                                      tGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGca                                                         <  1:419586/139‑1 (MQ=255)
                                                                                 tGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTa                                              <  2:569529/139‑1 (MQ=255)
                                                                                            ggAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggactt                                                                                                                              >  1:464933/1‑45 (MQ=255)
                                                                                            ggAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggactt                                                                                                                              <  2:464933/49‑5 (MQ=255)
                                                                                                   aCAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCt                            >  1:197140/1‑139 (MQ=255)
                                                                                                      aGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGt                         <  2:805980/139‑1 (MQ=255)
                                                                                                            gAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                   <  2:164699/91‑1 (MQ=255)
                                                                                                            gAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                   >  1:164699/1‑91 (MQ=255)
                                                                                                                       aaTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                                 <  2:1069345/66‑1 (MQ=255)
                                                                                                                       aaTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                                 >  1:1069345/1‑66 (MQ=255)
                                                                                                                        aTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTTTACCAATTCCGCCACCTTCGCACAGt                                                     >  2:37058/1‑93 (MQ=255)
                                                                                                                        aTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTTTACCAATTCCGCCACCTTCGCACAGt                                                     <  1:37058/93‑1 (MQ=255)
                                                                                                                            gTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTg   >  1:616496/1‑139 (MQ=255)
                                                                                                                             tGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGAc      <  2:613789/135‑1 (MQ=255)
                                                                                                                             tGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGAc      >  1:613789/1‑135 (MQ=255)
                                                                                                                             tGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGa  <  2:245341/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                      
TTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGA  >  NC_000913/4605957‑4606221

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