Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,347,615 G→A 7.4% T433M (ACG→ATG)  rnb ← RNase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,347,6150GA 7.4% 163.8 / 2.8 54T433M (ACG→ATG) rnbRNase II
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base A (2/2);  total (26/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 = 8.99e-01

CCGGGTTCAGTGCTAATTTCAGCAAATGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATG  >  NC_000913/1347468‑1347764
                                                                                                                                                   |                                                                                                                                                     
ccGGGTTCAGTGCTAATTTCAGCAAATGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGTGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGtt                                                                                                                                                     <  2:237887/150‑1 (MQ=255)
      tCAGTGCTAATTTCAGCAAATGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGc                                                                                                                                               <  1:10100/150‑1 (MQ=255)
      tCAGTGCTAATTTCAGCAAATGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGc                                                                                                                                               >  2:1642366/1‑150 (MQ=255)
       cAGTGCTAATTTCAGCAAATGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCa                                                                                                                                              <  2:1532694/150‑1 (MQ=255)
              aaTTTCAGCAAATGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGc                                                                                                                                       <  1:1786382/150‑1 (MQ=255)
              aaTTTCAGCAAATGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGc                                                                                                                                       <  1:852835/150‑1 (MQ=255)
                   cAGCAAATGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCCCGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAgcg                                                                                                                                  >  2:1321180/1‑150 (MQ=255)
                   cAGCAAATGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAgcg                                                                                                                                  >  2:1536375/1‑150 (MQ=255)
                   cAGCAAATGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAgcg                                                                                                                                  <  1:1151209/150‑1 (MQ=255)
                     gCAAATGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGcgcg                                                                                                                                <  1:743736/150‑1 (MQ=255)
                          tGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCACTTTGCAGAAACCGTACAGCGTGAGCACTTCTTCCGCATCGACATGCAGACCGTGCGATTTCAGCAACGCTGCCAGCTCGTCGGc                                                                                                                           >  2:1843665/1‑150 (MQ=255)
                              tGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTc                                                                                                                       <  1:183871/150‑1 (MQ=255)
                                 aaGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGcc                                                                                                                    <  2:1249296/150‑1 (MQ=255)
                                   gCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCgg                                                                                                                  <  2:1012330/150‑1 (MQ=255)
                                     gacgaATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGAt                                                                                                                <  2:733819/150‑1 (MQ=255)
                                     gacgaAGGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGAt                                                                                                                <  1:1755168/150‑1 (MQ=255)
                                      acgaATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATc                                                                                                               <  2:750255/150‑1 (MQ=255)
                                      acgaATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATc                                                                                                               <  2:743212/150‑1 (MQ=255)
                                      acgaATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGGCCAGCGGGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGGCGGCATTCGCCGGATc                                                                                                               <  1:536247/150‑1 (MQ=255)
                                           tGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGc                                                                                                          >  1:585035/1‑150 (MQ=255)
                                                cTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCata                                                                                                     >  1:1622734/1‑150 (MQ=255)
                                                  gTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATg                                                                                                   <  1:1753988/150‑1 (MQ=255)
                                                   tagAGGAAACCAGATGGTTGCGAGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGAACTTCTTCGGCATAGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGAGCGTCGGCATTCGCCGGATCAAAGCCCATATGc                                                                                                  >  2:685450/3‑150 (MQ=255)
                                                    cGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCa                                                                                                 <  2:194845/150‑1 (MQ=255)
                                                       ggACACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCGGGATCAAAGCCCATATGCACGt                                                                                              <  2:1132881/150‑1 (MQ=255)
                                                           aCCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTtata                                                                                          <  2:131646/150‑1 (MQ=255)
                                                             cAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGa                                                                                        >  1:423588/1‑150 (MQ=255)
                                                                  ggttgCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCAGATGCACGTTATAGATGCCa                                                                                   >  1:1959391/1‑150 (MQ=255)
                                                                            ccAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCt                                                                         >  2:118470/1‑150 (MQ=255)
                                                                            ccAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCt                                                                         >  2:125552/1‑150 (MQ=255)
                                                                             cAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCATTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCtt                                                                        >  2:1124795/1‑150 (MQ=255)
                                                                             cAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCATTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCtt                                                                        >  2:1895556/1‑150 (MQ=255)
                                                                                 tCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATcg                                                                    >  1:531374/1‑150 (MQ=255)
                                                                                  cACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATcgc                                                                   >  1:1808561/1‑150 (MQ=255)
                                                                                         cAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTAc                                                            >  1:1969574/1‑150 (MQ=255)
                                                                                              tGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCggg                                                       <  1:1321180/150‑1 (MQ=255)
                                                                                                cAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGcc                                                     >  2:1757297/1‑150 (MQ=255)
                                                                                                cAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGcc                                                     <  2:585035/150‑1 (MQ=255)
                                                                                                     aCCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGcaca                                                <  1:1263996/150‑1 (MQ=255)
                                                                                                       cGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGAGCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACaaa                                              <  2:563492/150‑1 (MQ=255)
                                                                                                        gTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGctaacat                                             >  1:1608014/1‑144 (MQ=255)
                                                                                                            aGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTa                                         <  2:531374/150‑1 (MQ=255)
                                                                                                             gCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAg                                        >  2:1428420/1‑150 (MQ=255)
                                                                                                                      cttcttCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAatca                               >  2:1206749/1‑150 (MQ=255)
                                                                                                                      cttcttCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACCTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAatca                               >  2:655173/1‑150 (MQ=255)
                                                                                                                        tcttcGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAatcatc                             >  1:77963/1‑150 (MQ=255)
                                                                                                                                 tCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTcg                    <  1:102101/150‑1 (MQ=255)
                                                                                                                                  cGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTcga                   >  1:1632954/1‑150 (MQ=255)
                                                                                                                                   gACATGCAGACCGTGCATTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTcgac                  <  1:1895556/150‑1 (MQ=255)
                                                                                                                                   gACATGCAGACCGTGCATTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTcgac                  <  1:1124795/150‑1 (MQ=255)
                                                                                                                                       tGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATa              >  2:1670086/1‑150 (MQ=255)
                                                                                                                                       tGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATa              >  2:478688/1‑150 (MQ=255)
                                                                                                                                          aGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACgg           >  1:1253087/1‑150 (MQ=255)
                                                                                                                                             ccGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTg        <  2:521501/150‑1 (MQ=255)
                                                                                                                                                   gTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCTGCAATCATCGCTTCTTCGACGATACGGTTGGCAATg  <  2:423588/150‑1 (MQ=255)
                                                                                                                                                   gTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATg  >  1:1919362/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGGGTTCAGTGCTAATTTCAGCAAATGACTGGAAGCGACGAATGCGGCTGTCGAGGAAACCAGTTGGTTGCGCGTCCAGTTCACGACGCAGTTTGCAGAAACCGTCCAGCGTGAGCACTTCTTCGGCATCGACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATG  >  NC_000913/1347468‑1347764

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