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.94.5% 207.8 / 4.7 55coding (186/921 nt)yedAamino acid exporter for phenylalanine, threonine
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (26/26);  total (27/28)
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

GTATATCATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029586‑2029858
                                                                                                                                          |                                                                                                                                      
gTATATCATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTc                                                                                                                                        >  1:1090297/1‑139 (MQ=255)
  atatCATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTccc                                                                                                                                      >  2:1514328/1‑139 (MQ=255)
           ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGt                                                                                                                            >  2:1421481/1‑139 (MQ=255)
           ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGt                                                                                                                            >  1:675468/1‑139 (MQ=255)
           ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGACTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGTCACCGGGCCGACCAAAACTCCCCCGCTACgg                                                                                                                            >  2:1481762/1‑138 (MQ=255)
            gggCCCAACTTATTTTTTCCTTCGGATTGGGGTGGAAAGCTGGCCCCCGTTATTGATGGCGGGCGTTCGACTCCTGGCAGCCGGTATTTTATTGCTGGAATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTc                                                                                                                           <  1:1481762/139‑1 (MQ=255)
             ggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTcc                                                                                                                          >  2:1518444/1‑139 (MQ=255)
              tcTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCg                                                                                                                         <  1:680643/138‑1 (MQ=255)
              gCTCAACCTATTTTGTCATTCGGATTGGCTTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTCTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCg                                                                                                                         <  2:1090297/139‑1 (MQ=255)
                             tCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCcgcg                                                                                                          >  2:775281/1‑139 (MQ=255)
                              cATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcc                                                                                                                                       >  2:613378/1‑110 (MQ=255)
                              cATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcc                                                                                                                                       <  1:613378/110‑1 (MQ=255)
                                 tCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGa                                                                                                      <  1:775281/139‑1 (MQ=255)
                                  cGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGAt                                                                                                     >  2:1292831/1‑139 (MQ=255)
                                    gATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTg                                                                                                   <  2:675468/139‑1 (MQ=255)
                                    gATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTg                                                                                                   <  1:1518444/139‑1 (MQ=255)
                                       tGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCg                                                                                                                         >  2:933237/1‑114 (MQ=255)
                                       tGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTCTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCg                                                                                                                         <  1:933237/114‑1 (MQ=255)
                                                   cTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCTGCGCTGATTGGCCTGTTATTgctgg                                                                                    >  2:1157590/1‑139 (MQ=255)
                                                   cTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctgg                                                                                    >  1:558006/1‑139 (MQ=255)
                                                   cTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctgg                                                                                    >  1:1216623/1‑139 (MQ=255)
                                                      gCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctggctg                                                                                 <  1:1292831/139‑1 (MQ=255)
                                                            gTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTa                                                                           <  2:1216623/139‑1 (MQ=255)
                                                                 tgatgGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGc                                                                      <  2:1606206/139‑1 (MQ=255)
                                                                     ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcccccg                                                                                                                                   <  2:717503/75‑2 (MQ=255)
                                                                     ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcccccg                                                                                                                                   >  1:717503/1‑74 (MQ=255)
                                                                     ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATgg                                                                  >  1:1459154/1‑139 (MQ=255)
                                                                     ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATgg                                                                  >  1:1381301/1‑139 (MQ=255)
                                                                     ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATgg                                                                  >  2:1016708/1‑139 (MQ=255)
                                                                             ttCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTccc                                                                                                                                      <  1:163931/64‑1 (MQ=255)
                                                                             ttCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTccc                                                                                                                                      >  2:163931/1‑64 (MQ=255)
                                                                                gATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCAGCTACGTCCGCTGCTCAATGCTGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGcc                                                       <  1:1157590/139‑1 (MQ=255)
                                                                                     cTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACa                                                  <  1:812524/139‑1 (MQ=255)
                                                                                     cTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACa                                                  <  2:1459154/139‑1 (MQ=255)
                                                                                          aGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCaaa                                             <  2:558006/139‑1 (MQ=255)
                                                                                            ccGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTgg                                                                                                  >  1:1532284/1‑84 (MQ=255)
                                                                                            ccGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTgg                                                                                                  <  2:1532284/84‑1 (MQ=255)
                                                                                            ccGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAAt                                           >  2:427405/1‑139 (MQ=255)
                                                                                                      aTTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGt                                         <  2:399436/131‑1 (MQ=255)
                                                                                                      aTTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGt                                         >  1:399436/1‑131 (MQ=255)
                                                                                                          cTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGc                                                                                                 <  2:524686/71‑1 (MQ=255)
                                                                                                          cTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGc                                                                                                 >  1:524686/1‑71 (MQ=255)
                                                                                                          cTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGt                                                                                             >  1:1438086/1‑75 (MQ=255)
                                                                                                          cTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGt                                                                                             <  2:1438086/75‑1 (MQ=255)
                                                                                                          cTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGtt                                        >  2:1252884/1‑128 (MQ=255)
                                                                                                          cTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGtt                                        <  1:1252884/128‑1 (MQ=255)
                                                                                                                       tACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgc                                                                                       >  2:1332503/1‑68 (MQ=255)
                                                                                                                       tACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgc                                                                                       <  1:1332503/68‑1 (MQ=255)
                                                                                                                         cTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGt              <  1:1016708/139‑1 (MQ=255)
                                                                                                                           gcgcGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCcgc                                                                                                           <  2:691007/44‑1 (MQ=255)
                                                                                                                           gcgcGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCcgc                                                                                                           >  1:691007/1‑44 (MQ=255)
                                                                                                                           gcgcGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGccgccg                     <  2:1544542/130‑1 (MQ=255)
                                                                                                                           gcgcGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGccgccg                     >  1:1544542/1‑130 (MQ=255)
                                                                                                                           gcgcGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTg            >  1:1110461/1‑139 (MQ=255)
                                                                                                                            cgcgGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgct                                                                                      <  1:807991/64‑1 (MQ=255)
                                                                                                                            cgcgGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgct                                                                                      >  2:807991/1‑64 (MQ=255)
                                                                                                                            cgcgGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGt                                                                                <  1:1310775/70‑1 (MQ=255)
                                                                                                                            cgcgGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGt                                                                                >  2:1310775/1‑70 (MQ=255)
                                                                                                                            cgcgGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGc           >  2:345318/1‑139 (MQ=255)
                                                                                                                            cgcgGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGc           >  2:1068368/1‑139 (MQ=255)
                                                                                                                                   acaAACCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTg    <  1:1139576/139‑1 (MQ=255)
                                                                                                                                     aaaCTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  2:1110461/139‑1 (MQ=255)
                                                                                                                                     aaaCTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  1:1421481/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                      
GTATATCATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029586‑2029858

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