Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,029,729 (C)6→5 coding (191/921 nt) yedA → amino acid exporter for phenylalanine, threonine

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,029,7240C.95.3% 165.3 / 0.8 43coding (186/921 nt)yedAamino acid exporter for phenylalanine, threonine
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (22/19);  total (23/20)
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.98e-01

GTATATCATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029586‑2029858
                                                                                                                                          |                                                                                                                                      
gTACATCATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGGTGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTc                                                                                                                                        <  2:279582/139‑1 (MQ=255)
           ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGTCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGt                                                                                                                            >  1:305204/1‑139 (MQ=255)
           ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGt                                                                                                                            >  1:1151438/1‑139 (MQ=255)
           ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGt                                                                                                                            >  1:240863/1‑139 (MQ=255)
           ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGt                                                                                                                            >  1:262987/1‑139 (MQ=255)
           ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGt                                                                                                                            >  2:283536/1‑139 (MQ=255)
            gggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTc                                                                                                                           <  2:287781/139‑1 (MQ=255)
              gCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCg                                                                                                                         <  2:288327/139‑1 (MQ=255)
              gCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCg                                                                                                                         <  2:240863/139‑1 (MQ=255)
                tCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCgct                                                                                                                       >  1:321561/1‑139 (MQ=255)
                        ttttGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATg                                                                                                               <  1:283536/139‑1 (MQ=255)
                              cATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAACCTCCCCCGCTACGTCCGCTGCTCAATGCCgcgc                                                                                                         >  2:710710/1‑139 (MQ=255)
                              cATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCgcgc                                                                                                         >  2:362288/1‑139 (MQ=255)
                                     aTTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGTCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTgg                                                                                                  <  2:305204/139‑1 (MQ=255)
                                            tGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTccccc                                                                                                                                    <  1:1058024/99‑1 (MQ=255)
                                            tGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTccccc                                                                                                                                    >  2:1058024/1‑99 (MQ=255)
                                                   cTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctgg                                                                                    >  1:994306/1‑139 (MQ=255)
                                                   cTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctgg                                                                                    <  1:1044544/139‑1 (MQ=255)
                                                   cTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctgg                                                                                    <  1:710710/139‑1 (MQ=255)
                                                                 tgatgGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGc                                                                      <  2:262987/139‑1 (MQ=255)
                                                                     ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATgg                                                                  >  1:951304/1‑139 (MQ=255)
                                                                      gcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctgg                                                                                    >  1:595939/1‑120 (MQ=255)
                                                                      gcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctgg                                                                                    <  2:595939/120‑1 (MQ=255)
                                                                      gcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGt                                                                 <  2:994306/139‑1 (MQ=255)
                                                                       cgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTg                                                                >  2:233037/1‑139 (MQ=255)
                                                                          gcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGAc                                                              >  2:466563/1‑138 (MQ=255)
                                                                          gcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGAc                                                              <  1:466563/138‑1 (MQ=255)
                                                                          gcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACg                                                             >  1:568081/1‑139 (MQ=255)
                                                                                         cAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCaa                                              <  2:568081/139‑1 (MQ=255)
                                                                                               gTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCaaa                                             <  2:711544/135‑1 (MQ=255)
                                                                                               gTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAAGGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCaaa                                             >  1:711544/1‑135 (MQ=255)
                                                                                                       ttGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCgg                                >  1:831640/1‑139 (MQ=255)
                                                                                                         gCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCcgc                                                                                                           <  1:482827/62‑1 (MQ=255)
                                                                                                         gCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCcgc                                                                                                           >  2:482827/1‑62 (MQ=255)
                                                                                                          cTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCAt                             >  1:162663/1‑139 (MQ=255)
                                                                                                          cTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCAt                             >  1:205484/1‑139 (MQ=255)
                                                                                                          cTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCAt                             >  2:867404/1‑139 (MQ=255)
                                                                                                                 ttttGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCgccgcc                      <  1:864446/139‑1 (MQ=255)
                                                                                                                    tGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTa                   <  1:867404/139‑1 (MQ=255)
                                                                                                                          tGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGtt             >  2:238848/1‑139 (MQ=255)
                                                                                                                           gcgcGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTg            <  2:831640/139‑1 (MQ=255)
                                                                                                                           gcgcGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTg            >  2:913277/1‑139 (MQ=255)
                                                                                                                            cgcgGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATcg                           <  2:395589/123‑1 (MQ=255)
                                                                                                                            cgcgGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATcg                           >  1:395589/1‑123 (MQ=255)
                                                                                                                                     aaaCTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  1:233037/139‑1 (MQ=255)
                                                                                                                                     aaaCTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  1:913277/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                      
GTATATCATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029586‑2029858

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