Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,812 A→G 10.9% intergenic (+26/‑80) PP_r01 → / → PP_5SB 5S ribosomal RNA/5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947181,8120AG10.9% 138.6 / 7.1 55intergenic (+26/‑80)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base G (4/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCG  >  NC_002947/181668‑181960
                                                                                                                                                |                                                                                                                                                      
tCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGTTGGTAGTGTGGGGTTTCCGCATGTGGGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGc                                                                                                                                                     <  4:212599/148‑1 (MQ=12)
  ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAgg                                                                                                                                                  >  4:299939/1‑149 (MQ=32)
       gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAggg                                                                                                                                            <  4:283735/149‑1 (MQ=32)
       gacCATAGAGCATTGGAAACACCTGATCCAATCCCGAACTCAGTAGTGAAACGATGCATAGCCGATGGTAGTGTGGGGTTTCACAATGTTAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAggg                                                                                                                                            >  2:154165/1‑149 (MQ=11)
         ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGt                                                                                                                                          <  6:258828/149‑1 (MQ=32)
         ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGt                                                                                                                                          >  4:64722/1‑149 (MQ=32)
               agCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTc                                                                                                                                    >  2:138150/1‑149 (MQ=32)
                 cATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCtt                                                                                                                                  <  3:64722/149‑1 (MQ=32)
                  aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCttt                                                                                                                                 <  6:288437/149‑1 (MQ=32)
                         ccaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTACCCATGTGAGAGTAGGTCATCGTAAAGATTCATTTCGCAAAACCCCTATCTGCACATGAAGGTAGGGGTTTTGTCTTTAAGTAg                                                                                                                           >  8:271538/1‑148 (MQ=12)
                            ccTGATCCCATCCCGAAATCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACACCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAAc                                                                                                                       >  6:240991/1‑149 (MQ=25)
                                atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTAc                                                                                                                    <  3:299939/148‑1 (MQ=38)
                                  cccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGt                                                                                                                                                 >  6:138482/1‑118 (MQ=32)
                                  cccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGt                                                                                                                                                 <  5:138482/118‑1 (MQ=32)
                                   ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACaga                                                                                                                 <  4:162667/148‑1 (MQ=40)
                                         cGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGc                                                                                                         >  7:208079/1‑149 (MQ=255)
                                            aatCATTAGTGAAACGATTCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACACCTATCTGCACATGCAGGTAGTGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCa                                                                                                        <  7:271538/145‑1 (MQ=255)
                                                agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAg                                                                                                                                              <  2:252958/106‑1 (MQ=11)
                                                agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAg                                                                                                                                              >  1:252958/1‑106 (MQ=11)
                                                 gtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACag                                                                                                                  <  8:265089/133‑1 (MQ=40)
                                                  tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTc                                                                                                >  4:182320/1‑149 (MQ=255)
                                                       aaaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAgg                                                                                           >  3:49672/1‑149 (MQ=255)
                                                                aTCGCCGATGGTAGTTGGGGGTTTCCCCATGTGAGAGTACGTCAGCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCACGCACGTCCTAGGGACGGGCCa                                                                                  <  8:287046/149‑1 (MQ=18)
                                                                         ggTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAg                                                                          >  6:20710/1‑148 (MQ=255)
                                                                         ggTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAg                                                                          >  2:91383/1‑148 (MQ=255)
                                                                           tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCAc                                                                                                   <  3:154849/121‑1 (MQ=255)
                                                                           tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCAc                                                                                                   >  4:154849/1‑121 (MQ=255)
                                                                           tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAt                                                                       <  1:138150/149‑1 (MQ=255)
                                                                              tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGATATTCTCAGGCACGTCCTAGGGACTGGCCAGCAAACAGAAtta                                                                     >  3:4319/1‑147 (MQ=255)
                                                                                  gggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTg                                                                 <  4:49672/148‑1 (MQ=255)
                                                                                         cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                          >  8:129026/1‑148 (MQ=255)
                                                                                                gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCCCGTccctagg                                                                                          >  4:182689/1‑104 (MQ=255)
                                                                                                gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCCCGTccctagg                                                                                          <  3:182689/109‑6 (MQ=255)
                                                                                                 agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCTTCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAt                                                  <  1:91383/149‑1 (MQ=255)
                                                                                                 agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt                                                 <  1:154165/149‑1 (MQ=255)
                                                                                                     tAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGaa                                             >  6:333983/1‑149 (MQ=255)
                                                                                                         tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCGCACAGAATTTCTTGACGACCATaga                                                      <  2:40093/136‑1 (MQ=255)
                                                                                                         tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCGCACAGAATTTCTTGACGACCATaga                                                      >  1:40093/1‑136 (MQ=255)
                                                                                                           aTCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAaccacc                                        <  3:182320/148‑1 (MQ=255)
                                                                                                                 aaGATTGATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcc                                  <  5:240991/148‑1 (MQ=255)
                                                                                                                 aaGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTga                                                                <  2:120881/118‑1 (MQ=255)
                                                                                                                 aaGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTga                                                                >  1:120881/1‑118 (MQ=255)
                                                                                                                    aTTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatc                              >  2:49864/1‑149 (MQ=255)
                                                                                                                       ttaaagacAAAACCCCTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatccc                            >  2:213133/8‑148 (MQ=37)
                                                                                                                         tttCGCAAAACCCCTATCTGCACATGCAGGTAGGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTgacga                                                             >  1:46020/1‑114 (MQ=255)
                                                                                                                         tttCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTgacga                                                             <  2:46020/113‑1 (MQ=255)
                                                                                                                         tttCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                         >  4:321046/1‑149 (MQ=255)
                                                                                                                         tttCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                         >  7:61917/1‑149 (MQ=255)
                                                                                                                              cAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag                    <  8:61917/149‑1 (MQ=255)
                                                                                                                               aaaaCCCCTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag                    >  4:99470/1‑148 (MQ=255)
                                                                                                                                aaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag                    >  1:266289/1‑147 (MQ=255)
                                                                                                                                   ccccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat                                    >  4:227254/1‑128 (MQ=255)
                                                                                                                                   ccccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat                                    <  3:227254/128‑1 (MQ=255)
                                                                                                                                         tCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAAc           <  5:20710/147‑1 (MQ=255)
                                                                                                                                         tCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAAc           <  3:321046/147‑1 (MQ=255)
                                                                                                                                             cacaTGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATg       <  2:266289/147‑1 (MQ=255)
                                                                                                                                             cacaTGCAGGTAGGGGTTTTGTATTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCAAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt        <  4:4319/146‑1 (MQ=255)
                                                                                                                                                aTGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCg  >  2:301889/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
TCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCG  >  NC_002947/181668‑181960

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