Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000964 4,155,395 (A)5→6 intergenic (‑186/+38) trnY‑Lys ← / ← purA tRNA‑Lys/adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009644,155,3951.A97.1% 117.4 / ‑2.5 35intergenic (‑186/+38)trnY‑Lys/purAtRNA‑Lys/adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (21/13);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

CGAACCTTCGACCCTCTGATTAAAAGTCAGATGCTCTACCAACTGAGCTAATGGCTCTTAAAGATGTCGGAATTTCTATAGAGCAAACGACTTAAATCTAGTTCTCATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAA‑CAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTAATATCCTCAGTCCAGCCCGGCATTTCTTCATATACCGGCTCAC  >  NC_000964/4155153‑4155629
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          
cGAACCTTCGACCCTCTGATTAAAAGTCAGATGCTCTACCAACTGAGCTAATGGCTCTTAAAGATGTCGGAATTTCTATAGAGCAAACGACTTAAATCTAGTTCTCATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGccc                                                                                                                                                                                                                                     <  2:55324/251‑1 (MQ=255)
   aCCTTCGGCCCGCTGATTAAAAGGCAGATGCTCTACCAACTGAGCTAATGGCTCTTAAAGATGTCGGAATTTCTATCGAGCAAACGACTTAAATCTAGTTCTCATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATATCTTTTGTGGTTTTTTTTCCTATGTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCtt                                                                                                                                                                                                                                  <  2:232177/251‑1 (MQ=255)
          aCCCTCTGATTAAAAGTCAGATGCTCTACCAACTGAGCTAATGGCTCTTAAAGATGTCGGAATTTCTATAGAGCAAACGACTTAAATCTAGTTCTCATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAgg                                                                                                                                                                                                                           >  1:302082/1‑251 (MQ=255)
                     aaaaGTCAGATGCTCTACCAACTGAGCTAATGGCTCTTAAAGATGTCGGAATTTCTATAGAGCAAACGACTTAAATCTAGTTCTCATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGAc                                                                                                                                                                                                                >  1:79179/1‑251 (MQ=255)
                         gTCAGATGCTCTACCAACTGAGCTAATGGCTCTTAAAGATGTCGGAATTTCTATAGAGCAAACGACTTAAATCTAGTTCTCATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACatat                                                                                                                                                                                                            <  1:328086/251‑1 (MQ=255)
                                                  atgcctctTAAAGATGTCGGAATTTCTATAGAGCAAACGACTTAAATCTAGTTCTCATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTacacac                                                                                                                                                                                   >  2:223915/5‑251 (MQ=255)
                                                    ggCTCTTAAAGATGTCGGAATTTCTATAGAGCAAACGACTTAAATCTAGTTCTCATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTg                                                                                                                                                                                 >  1:332208/1‑251 (MQ=255)
                                                                 gTCGGAATTTCTATAGAGCAAACGACTTAAATCTAGTTCTCATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACAtttg                                                                                                                                                                    <  1:169461/251‑1 (MQ=255)
                                                                        tttCTATAGAGCAAACGACTTAAATCTAGTTCTCATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTtgtgg                                                                                                                                                             >  1:26047/1‑250 (MQ=255)
                                                                                         aCTTAAATCTAGTTCTCATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGAc                                                                                                                                              >  1:203690/1‑249 (MQ=255)
                                                                                                         caacTTGTTCTCTAATGAAGCCATTTCCCTAACAACGGGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCgg                                                                                                                            <  1:69999/248‑1 (MQ=255)
                                                                                                           tCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGa                                                                                                                           >  2:300979/1‑250 (MQ=255)
                                                                                                              tgtacgttATTGAAGCAGTTACAGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCCATAACTTTTGTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGATTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTTAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAAt                                                                                                                         <  2:167476/242‑1 (MQ=255)
                                                                                                               gTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGc                                                                                                                       <  2:145799/250‑1 (MQ=255)
                                                                                                                        tGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTg                                                                                                             >  1:150076/1‑251 (MQ=255)
                                                                                                                                     gTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACCTTTGTTTGTGCGCGGTCTGGACCGACAGAGAAACTAGAAAGCGGAATGCCTGTCAGATGGGACACACGCTcc                                                                                                >  2:343151/1‑250 (MQ=255)
                                                                                                                                                   tataATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGgcgcgc                                                                                  >  1:74531/1‑251 (MQ=255)
                                                                                                                                                               aaCAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAg                                                                      <  1:5911/251‑1 (MQ=255)
                                                                                                                                                                      aGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAgctcgctc                                                               >  1:242735/1‑251 (MQ=255)
                                                                                                                                                                              cAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAATCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCtt                                                       >  2:138267/1‑251 (MQ=255)
                                                                                                                                                                                  aaCTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGCAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGAGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGc                                                    >  2:152288/1‑250 (MQ=255)
                                                                                                                                                                                   aCTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCAc                                                  <  1:74227/251‑1 (MQ=255)
                                                                                                                                                                                          tGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTAat                                           >  1:119545/1‑251 (MQ=255)
                                                                                                                                                                                                     cATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGAATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTAATATCCTCAGTcc                                >  2:78029/1‑251 (MQ=255)
                                                                                                                                                                                                     cATATTTTTTTCGCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTAATATCCTCAGTcc                                <  1:6439/251‑1 (MQ=255)
                                                                                                                                                                                                        atTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTAATATCCTCAGTCCAGc                             <  1:55272/251‑1 (MQ=255)
                                                                                                                                                                                                           tttttCTCAGACCAACTTATATGCGGCTTTCAAACACAAAA‑CAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAGGCGGAATGCCTGTCAGCTGAGACACACGATAAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTAATATCCTCAGTCCAGCCCg                          >  2:67841/1‑251 (MQ=255)
                                                                                                                                                                                                           tttttCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTAATATCCTCAGTCCAGCCCg                          >  1:178790/1‑251 (MQ=255)
                                                                                                                                                                                                              ttCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTATTACCCCCATTCCAGCCCGGct                       >  1:156302/1‑250 (MQ=255)
                                                                                                                                                                                                                      ccAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATg                                                                                        <  1:75574/178‑1 (MQ=255)
                                                                                                                                                                                                                      ccAACTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATg                                                                                        >  2:75574/1‑178 (MQ=255)
                                                                                                                                                                                                                         aCTTATATGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTAATATCCTCAGTCCAGCCCGGCATTTCTTCatat            >  2:346903/1‑251 (MQ=255)
                                                                                                                                                                                                                         aCTTATACGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTAATATCCTCAGTCCAGCCCGGCATTTCTTCatat            <  1:317149/251‑1 (MQ=255)
                                                                                                                                                                                                                               atGCGGCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTAATATCCTCAGTCCAGCCCGGCATTTCTTCATATACCGGc      >  1:33186/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ggCTTTCAAACA‑AAAAACAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTAATATCCTCAGTCCAGCCCGGCATTTCTTCATATACCGGCTcac  >  1:326828/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          
CGAACCTTCGACCCTCTGATTAAAAGTCAGATGCTCTACCAACTGAGCTAATGGCTCTTAAAGATGTCGGAATTTCTATAGAGCAAACGACTTAAATCTAGTTCTCATCTTGTACTCTATTGAAGCAGTTTCCGTAACAACGAGATTTATAATATCACGAACAAATAGAAGTGGCAATAACTTTTTTGGTTTTTTTTCATATTTTTTTCTCAGACCAACTTATATGCGGCTTTCAAACA‑AAAA‑CAAGCCCCTTAAATAGGGGCTTGCAGACATATTCTATTTAGTTCGCACGGTACACACTGCGAAGGACATTTGTTTGTGAGCGGTCTGGACCGACAGAGAAAATAGAAAGCGGAATGCCTGTCAGCTGAGACACACGCTCAAGATAATGGCGCGCATTTTCCGGAAGCTCGCTCAAGCTCTTCGCACCTGTAATATCCTCAGTCCAGCCCGGCATTTCTTCATATACCGGCTCAC  >  NC_000964/4155153‑4155629

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