Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,466,100 A→T 7.7% noncoding (103/119 nt) PP_mr39 ← P15 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,466,1000AT 7.7% 158.2 / 2.9 52noncoding (103/119 nt)PP_mr39P15 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base T (2/2);  total (27/25)
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.97e-01

GATCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGG  >  NC_002947/3465951‑3466247
                                                                                                                                                     |                                                                                                                                                   
gATCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGaa                                                                                                                                                     <  1:76409/150‑1 (MQ=255)
     gAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGTGGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCGGa                                                                                                                                                >  1:580403/1‑150 (MQ=255)
       ggcggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGAcc                                                                                                                                              <  1:716323/150‑1 (MQ=255)
         cggcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCgg                                                                                                                                            <  1:561408/150‑1 (MQ=255)
           gcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCgggg                                                                                                                                          <  1:197369/150‑1 (MQ=255)
           gcAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCgggg                                                                                                                                          <  1:518662/150‑1 (MQ=255)
             aGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGCTTCCCGACCGGGGGt                                                                                                                                        <  1:960407/150‑1 (MQ=255)
                 tGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGa                                                                                                                                    <  2:830725/150‑1 (MQ=255)
                     ccGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTAt                                                                                                                                <  2:793091/150‑1 (MQ=255)
                         aCGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATtctc                                                                                                                            <  2:14419/150‑1 (MQ=255)
                             gATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGTTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGt                                                                                                                        >  1:369089/1‑150 (MQ=255)
                                  gCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcggc                                                                                                                   >  1:314030/1‑150 (MQ=255)
                                     aaaTCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCggcggc                                                                                                                >  1:220178/1‑150 (MQ=255)
                                         cGCCAGGCGTAGGTGTTTGCCTACTAATAAGTTGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCt                                                                                                            <  2:580403/150‑1 (MQ=255)
                                          gCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTc                                                                                                           >  1:596608/1‑150 (MQ=255)
                                           ccAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCg                                                                                                          >  1:36276/1‑150 (MQ=255)
                                             aGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGaa                                                                                                        <  2:596608/150‑1 (MQ=255)
                                                 gTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTg                                                                                                    <  2:36276/150‑1 (MQ=255)
                                                 gTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTg                                                                                                    >  1:189843/1‑150 (MQ=255)
                                                  tAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTgg                                                                                                   >  1:169639/1‑150 (MQ=255)
                                                    ggTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGAc                                                                                                 <  2:369089/150‑1 (MQ=255)
                                                       gtTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGcc                                                                                              >  2:1095694/1‑150 (MQ=255)
                                                             gtagtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCgggg                                                                                                                                          <  1:575933/100‑1 (MQ=255)
                                                             gtagtaATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCgggg                                                                                                                                          >  2:575933/1‑100 (MQ=255)
                                                                                  cATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTaa                                                                   >  1:457001/1‑150 (MQ=255)
                                                                                   aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTaaa                                                                  <  2:1025688/150‑1 (MQ=255)
                                                                                   aTGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTaaa                                                                  >  2:462324/1‑150 (MQ=255)
                                                                                       aTCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCaa                                                              <  2:189843/150‑1 (MQ=255)
                                                                                              cTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAAc                                                                 <  1:514763/140‑1 (MQ=255)
                                                                                              cTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAAc                                                                 >  2:514763/1‑140 (MQ=255)
                                                                                               tCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCaaa                                                      >  1:757759/1‑150 (MQ=255)
                                                                                                             tgtTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAg                                        >  1:881553/1‑150 (MQ=255)
                                                                                                              gtTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAgg                                       >  2:399421/1‑150 (MQ=255)
                                                                                                               ttGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcggc                                                                                                                   >  1:235044/1‑73 (MQ=255)
                                                                                                               ttGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTcggc                                                                                                                   <  2:235044/73‑1 (MQ=255)
                                                                                                                     aaaGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACg                                <  1:1041425/150‑1 (MQ=255)
                                                                                                                         gCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGc                            >  1:656711/1‑150 (MQ=255)
                                                                                                                           cTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCg                          >  2:1136115/1‑150 (MQ=255)
                                                                                                                           cTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCg                          <  1:399421/150‑1 (MQ=255)
                                                                                                                            tGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTTAGGATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGgc                                                                                <  2:296397/95‑1 (MQ=255)
                                                                                                                            tGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGgc                                                                                >  1:296397/1‑95 (MQ=255)
                                                                                                                              aatactcaaacCCCCGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCa                       >  2:966269/8‑150 (MQ=38)
                                                                                                                                tactcaaacCCCCGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATa                                                <  2:1088832/118‑1 (MQ=38)
                                                                                                                                tactcaaacCCCCGGTCGGCTTCCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATa                                                >  1:1088832/6‑123 (MQ=38)
                                                                                                                                 aacaaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATc                                                                                   >  1:698615/1‑87 (MQ=255)
                                                                                                                                 aacaaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATc                                                                                   <  2:698615/87‑1 (MQ=255)
                                                                                                                                    aaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGcc                 >  2:749427/1‑150 (MQ=255)
                                                                                                                                    aaaacCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGcc                 >  1:400708/1‑150 (MQ=255)
                                                                                                                                           ccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGc                                    >  1:939605/1‑124 (MQ=255)
                                                                                                                                           ccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGc                                    <  2:939605/124‑1 (MQ=255)
                                                                                                                                           ccGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAg          >  1:232353/1‑150 (MQ=255)
                                                                                                                                            cGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAgg         <  2:656711/150‑1 (MQ=255)
                                                                                                                                                   gAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAgg  <  2:457001/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GATCGGAGGCGGCAGGCTGGCCCGAACGGGATTGGCTAAATCGCCAGGCGTAGGTGTTTGCGTAGTAATAAGTGGCGTAGTTCATGAATCTGTCCTCAAAGTGAATCGGTGTTGGTCAAAGGCCTGAAAAACAAAACCCCCGGTCGGGAAGCCGACCGGGGGTTTGAGTATTCTCATGTTCGGCGGCCCCTCGAAGGTGGACGCCGGGTGGGTATCGGCGCCTAGTGGCTAAACCAATACCCAAAGTAATAGCTGGCAGGTGCCACGCAGCCGGCAACGGCCAGCACAGGGCGCAGG  >  NC_002947/3465951‑3466247

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