Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,796,337 C→T *150* (TAG→TAA ygaU ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,796,3370CT98.2% 168.5 / ‑4.3 54*150* (TAG→TAAygaUuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (29/24);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

GGAGGATAAACACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTC  >  NC_000913/2796201‑2796466
                                                                                                                                        |                                                                                                                                 
ggaggaTAAACACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAattaatta                                                                                                                                 <  2:1185548/139‑1 (MQ=255)
   ggaTAAACACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTActc                                                                                                                              <  2:176921/139‑1 (MQ=255)
        aaCACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCg                                                                                                                         <  2:1440250/139‑1 (MQ=255)
                         aCCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGc                                                                                                        >  1:1055383/1‑139 (MQ=255)
                         aCCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGc                                                                                                        <  2:901005/139‑1 (MQ=255)
                         aCCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGc                                                                                                        >  2:961364/1‑139 (MQ=255)
                         aCCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCGTCCGGAATACGCAACACTTGc                                                                                                        >  1:901005/1‑139 (MQ=255)
                           ccTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGccc                                                                                                      >  1:716609/1‑139 (MQ=255)
                             tGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGa                                                                                                                       >  1:1158781/1‑120 (MQ=255)
                             tGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGa                                                                                                                       <  2:1158781/120‑1 (MQ=255)
                             tGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCg                                                                                                    >  2:464352/1‑139 (MQ=255)
                                cTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGc                                                                                                                 >  1:1354491/1‑123 (MQ=255)
                                cTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGc                                                                                                                 <  2:1354491/123‑1 (MQ=255)
                                     gcggTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGa                                                                                                                       >  2:370466/1‑112 (MQ=255)
                                     gcggTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGa                                                                                                                       <  1:370466/112‑1 (MQ=255)
                                        gTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTa                                                                                         >  1:382697/1‑139 (MQ=255)
                                         tACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTAt                                                                                        <  2:382697/139‑1 (MQ=255)
                                         tACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTACCCCGGATAAATTTTAt                                                                                        <  1:961364/139‑1 (MQ=255)
                                              gcgcGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAattaatta                                                                                                                                 >  2:1118804/1‑93 (MQ=255)
                                              gcgcGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAattaatta                                                                                                                                 <  1:1118804/93‑1 (MQ=255)
                                                   agaTACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGa                                                                                                                       <  1:1081468/98‑1 (MQ=255)
                                                   agaTACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGa                                                                                                                       >  2:1081468/1‑98 (MQ=255)
                                                        ccGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCa                                                                         >  2:1158136/1‑139 (MQ=255)
                                                          gcgcgcgcCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATc                                                                       <  2:716609/139‑1 (MQ=255)
                                                            gcgcgcCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCa                                                                                                                >  1:1133892/1‑96 (MQ=255)
                                                            gcgcgcCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCa                                                                                                                <  2:1133892/96‑1 (MQ=255)
                                                                            aTGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAattaatta                                                                                                                                 >  2:698636/1‑63 (MQ=255)
                                                                                 ccAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCa                                                                         <  2:985131/114‑1 (MQ=255)
                                                                                 ccAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCa                                                                         >  1:985131/1‑114 (MQ=255)
                                                                                   aGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTAt                                              >  2:361510/1‑139 (MQ=255)
                                                                                               cGCCTATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCg                                                                                     <  1:818057/88‑1 (MQ=255)
                                                                                               cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCggg                                                                                   >  2:1437342/1‑90 (MQ=255)
                                                                                               cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCggg                                                                                   <  1:1437342/90‑1 (MQ=255)
                                                                                               cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCg                                                                                     >  2:818057/1‑88 (MQ=255)
                                                                                               cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCAt                                                                        >  2:565189/1‑101 (MQ=255)
                                                                                               cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCAt                                                                        <  1:565189/101‑1 (MQ=255)
                                                                                               cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAg                                  >  1:81654/1‑139 (MQ=255)
                                                                                                          aTTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGt                               >  2:29502/1‑131 (MQ=255)
                                                                                                          aTTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGt                               <  1:29502/131‑1 (MQ=255)
                                                                                                              tCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGc                                                                                                                 <  1:109118/45‑1 (MQ=255)
                                                                                                              tCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGc                                                                                                                 >  2:109118/1‑45 (MQ=255)
                                                                                                                  gcagGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGttt               >  1:326483/1‑139 (MQ=255)
                                                                                                                   cgggCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTg              >  1:726744/3‑139 (MQ=255)
                                                                                                                   cagGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGc                                                                          >  2:150130/1‑79 (MQ=255)
                                                                                                                   cagGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGc                                                                          <  1:150130/79‑1 (MQ=255)
                                                                                                                    agGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGcc                                                                                                       >  1:12557/1‑49 (MQ=255)
                                                                                                                    agGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGcc                                                                                                       <  2:12557/49‑1 (MQ=255)
                                                                                                                      gCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTgg             <  1:456993/137‑1 (MQ=255)
                                                                                                                      gCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTgg             >  2:456993/1‑137 (MQ=255)
                                                                                                                        gATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCg                                                                                     >  1:1309820/1‑63 (MQ=255)
                                                                                                                        gATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCg                                                                                     <  2:1309820/63‑1 (MQ=255)
                                                                                                                        gATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAAt         >  2:311392/1‑139 (MQ=255)
                                                                                                                        gATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAAt         >  1:441175/1‑139 (MQ=255)
                                                                                                                          tggctggcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATgg       <  2:1042856/139‑1 (MQ=255)
                                                                                                                            gctggcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCa     >  1:451856/1‑139 (MQ=255)
                                                                                                                               ggcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTc  <  1:311392/139‑1 (MQ=255)
                                                                                                                                gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTc                                                         <  2:1020160/83‑1 (MQ=255)
                                                                                                                                gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTc                                                         >  1:1020160/1‑83 (MQ=255)
                                                                                                                                        |                                                                                                                                 
GGAGGATAAACACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTC  >  NC_000913/2796201‑2796466

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