Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,797,369 G→T 10.2% I192I (ATC→ATA)  pheT ← phenylalanine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,797,3690GT10.2% 132.8 / 3.9 53I192I (ATC→ATA) pheTphenylalanine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base T (2/4);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAAT  >  NC_000913/1797221‑1797518
                                                                                                                                                    |                                                                                                                                                     
cGCAACGACGCAGTTTTTCTTTCATCCACAGCTGAGTTGGCGCGTTAACGTTACTGCCTTTTACCCCACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGa                                                                                                                                                      >  1:637810/1‑150 (MQ=255)
cGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGa                                                                                                                                                      >  1:426357/1‑150 (MQ=255)
  cAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGAtt                                                                                                                                                    <  1:342213/150‑1 (MQ=255)
     cgacgCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTcc                                                                                                                                                 <  1:1281894/150‑1 (MQ=255)
       acgCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCgg                                                                                                                                               <  2:428878/150‑1 (MQ=255)
       acgCAGTTTTCCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCgg                                                                                                                                               <  1:487518/150‑1 (MQ=255)
        cgCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGt                                                                                                                                              >  1:334245/1‑150 (MQ=255)
          cAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTg                                                                                                                                            >  2:20975/1‑150 (MQ=255)
          cAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGACTTCCGGTTg                                                                                                                                            <  2:267905/150‑1 (MQ=255)
           aGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAAATCGTTCCGGTTgg                                                                                                                                           <  2:3632/150‑2 (MQ=255)
           aGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCAGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAAATCGTTCCGGTTgg                                                                                                                                           <  1:824453/150‑2 (MQ=255)
            gTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCTCACCAACCGGAACGATTTCCGGTTGaa                                                                                                                                          >  2:115534/1‑150 (MQ=255)
             tttttCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAAATCGTTCCGGTTggtg                                                                                                                                         <  1:434151/150‑4 (MQ=255)
              ttttCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCAATGGTCGCACCAACCGGAACGATTTCCGGTTGAAcc                                                                                                                                        <  2:1256467/150‑1 (MQ=255)
                 tctttcATCCACAGCGGAGTTGGCGCGTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACTATTTCCGGTTGAACCAGc                                                                                                                                     <  2:28593/150‑1 (MQ=255)
                    ttcATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGc                                                                                                                                  <  1:1248591/150‑1 (MQ=255)
                         tcacaGAGGAGTTGGGGTTTTAATGTTTATTTCTTTTACCACAGGGACTAGATAACGCGTGCTGGGTTCGGGCTCTTCGTCTGTAATCTGGAGCGTGTCTTCGATGGTCGCACCATCGGTAAAGTTTTCAGGTTGAACCAGCGGCAGCTg                                                                                                                             <  2:966880/149‑1 (MQ=255)
                                  aGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCTGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCAc                                                                                                                    <  2:864567/150‑1 (MQ=255)
                                          cTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGt                                                                                                            >  2:570317/1‑150 (MQ=255)
                                          cTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGt                                                                                                            >  1:1227175/1‑150 (MQ=255)
                                            ttAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCa                                                                                                          <  1:1175154/150‑1 (MQ=255)
                                             tAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCAc                                                                                                         >  2:812795/1‑150 (MQ=255)
                                                  ttAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCAGGCGCTTATACTGTAATCGGCAGCGTGTCGTCGATGGTGTCACAAACAGGAACTATTTCAGTTTGAACCATCGGCAGCTGTTTCAGCACGGCAACGTAACGCGCa                                                                                                    >  2:858943/1‑150 (MQ=255)
                                                       gCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACAcc                                                                                               <  2:730119/150‑1 (MQ=255)
                                                               ccACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGc                                                                                       >  2:1258892/1‑150 (MQ=255)
                                                                  cacGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTa                                                                                    <  2:375159/150‑1 (MQ=255)
                                                                  cacGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTa                                                                                    <  1:812795/150‑1 (MQ=255)
                                                                       ccAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGGCGTCGATGGGCGCACCACCCGGACCGATTTCCGGTTCAACCACCGCCAGAGGGTGCAGCACGCCAACGACCCGCCCAACACCACGGATGCGTacgagg                                                                               >  2:738992/1‑145 (MQ=255)
                                                                           gATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCgg                                                                           <  2:637810/150‑1 (MQ=255)
                                                                                       aGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCg                                                               <  1:8491/150‑1 (MQ=255)
                                                                                       aGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCg                                                               <  2:348978/150‑1 (MQ=255)
                                                                                       aGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCg                                                               <  2:1285074/150‑1 (MQ=255)
                                                                                       aGGCTTCCGGCGCTTCGACTGTAAGCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCg                                                               <  1:989044/150‑1 (MQ=255)
                                                                                        ggCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGt                                                              >  2:557903/1‑150 (MQ=255)
                                                                                               ggCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGcgg                                                       >  1:397432/1‑148 (MQ=255)
                                                                                                gcgcTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGa                                                      >  2:1145703/1‑150 (MQ=255)
                                                                                                gcgcTTCGACTGTAATCGGCAGCGTGTCGTCCATGCTCGCACAAACCGGAACGATTTCCGGTTGAACCAGCGCAAGCTGATTCAGCACGGCAACGTCACCCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGa                                                      >  1:463285/1‑150 (MQ=255)
                                                                                                     tCGACTGTAATCGGCAGCGTGTCGGCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCg                                                 <  2:1062505/150‑1 (MQ=255)
                                                                                                        aCTGTAACCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCACCACGGCAACGTCACGCGCAACACCAAAGATGCCTAAGAAGTCGGCACGGTTTGCCGTCACGCTGATGTCGATg                                              >  1:140421/1‑150 (MQ=255)
                                                                                                              aTCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttg                                        >  2:421229/1‑150 (MQ=255)
                                                                                                              aTCGGCAGCGTGGCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttg                                        <  1:77927/150‑1 (MQ=255)
                                                                                                                cGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGGGTTGtc                                      <  2:698332/150‑1 (MQ=255)
                                                                                                                 ggCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGtca                                     <  1:90209/150‑1 (MQ=255)
                                                                                                                           cgtcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACCCCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCa                           >  2:504020/1‑150 (MQ=255)
                                                                                                                            gtcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAg                          >  1:104988/1‑150 (MQ=255)
                                                                                                                             gcgATGGTCGCACCAAGTTTAAAGATTTCCGTTTGAACCAGCGGCTGCTGGTTCAGCACGGCATCGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCAGTGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAgg                         <  2:766794/149‑1 (MQ=255)
                                                                                                                                aTGGTCGCAACACCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATTATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTAt                      <  1:117542/150‑1 (MQ=255)
                                                                                                                                  ggTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCGTAAGCAGTCGGCCCGGGTGGGCGTCCCGCTGATTTCGATGGTGTTGTCATCAAGTGTCAGGTATTc                    >  1:353833/1‑150 (MQ=255)
                                                                                                                                   gTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCa                   <  2:713290/150‑1 (MQ=255)
                                                                                                                                      gCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCACCACCAATGCTGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACgg                >  1:941423/1‑150 (MQ=255)
                                                                                                                                       gttcAACCGGAAATCGTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGa               >  2:1064949/4‑150 (MQ=255)
                                                                                                                                         ccAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCCCGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGata             >  1:1283514/1‑150 (MQ=255)
                                                                                                                                         ccAACCGGAACGATTTCCGGTTGAACCAGCGGAAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCAAGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGata             >  1:557089/1‑150 (MQ=255)
                                                                                                                                                    gATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAAt  >  1:32130/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCAACGACGCAGTTTTTCTTTCATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAAT  >  NC_000913/1797221‑1797518

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