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% 108.9 / 3.7 44intergenic (+26/‑80)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (20/19);  new base G (3/2);  total (23/21)
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.92e-01

CGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAT  >  NC_002947/181674‑181954
                                                                                                                                           |                                                                                                                                              
cgacCATAGAGCATTGGCACCACCTGTTACTATGCCGAACTCAGGAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATACTCAAGATTCATTACGTAAAAACCCCATCTCCACATCCAGATaagc                                                                                                                                     <  2:520749/150‑4 (MQ=11)
cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAggg                                                                                                                                     <  2:5536/150‑1 (MQ=11)
   ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGtt                                                                                                                                  >  2:725320/1‑150 (MQ=32)
    cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGttt                                                                                                                                 >  1:830621/1‑150 (MQ=32)
        gagCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTc                                                                                                                             <  2:830621/150‑1 (MQ=32)
             ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                        <  2:610807/150‑1 (MQ=32)
             ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATCCAGGTAGGGGTTTTGTCTTTaa                                                                                                                        >  2:101492/1‑150 (MQ=17)
                   ccaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                        <  2:800106/144‑1 (MQ=32)
                   ccaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                                        >  1:800106/1‑144 (MQ=32)
                     accTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAAc                                                                                                                >  2:328605/1‑150 (MQ=35)
                       cTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTa                                                                                                              >  2:602490/1‑150 (MQ=38)
                             ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACaga                                                                                                          >  1:556417/1‑148 (MQ=40)
                             ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACaga                                                                                                          <  2:556417/148‑1 (MQ=40)
                              catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGAt                                                                                                       <  1:220597/150‑1 (MQ=255)
                                            tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTcc                                                                                         <  1:725320/150‑1 (MQ=255)
                                              gtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTa                                                                                       <  1:328605/150‑1 (MQ=255)
                                               tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAg                                                                                      >  1:686774/1‑150 (MQ=255)
                                                       tGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGc                                                                              <  2:310549/150‑1 (MQ=255)
                                                          aTCGCCGATGGTAGTGTGGGGTTTCCGCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACGCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAg                                                                           >  2:530835/1‑150 (MQ=255)
                                                            cGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTTGGTCATCGTCAAGATTCAGTACGGAGTACCCCTATGTGCGCATGTAGGGAGGGGTGTTGTATTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGca                                                                         >  2:470778/1‑150 (MQ=18)
                                                             gctgATGGTAGAGTCGGGGTTCCCCATGTGAGAGGAGGTCATCGTCAAGATTCATTTCCAAAAACCCCTATCTGCACATGCAGGGAGGGGTTTTGTCTTTAAGTAGTAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGGCAggaa                                                                        >  2:613011/4‑147 (MQ=18)
                                                                gATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCCAGCTTCTATTTCAAAAACCCCTACCGGGAGATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGGAGGCACGTCCTAGGGACGGGCCAGcacaca                                                                     <  1:101492/150‑1 (MQ=21)
                                                                         gtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCtt                                                            <  1:602490/150‑1 (MQ=255)
                                                                          tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAg                                                                                                                                       <  2:802764/74‑1 (MQ=35)
                                                                          tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAg                                                                                                                                       >  1:802764/1‑74 (MQ=35)
                                                                                   cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACGTTCAGGTAGGGTTTTTGTCTTTAAGTATGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGACAGCGCACAGAATTTCTTGACGACCATa                                                  >  2:390132/1‑150 (MQ=25)
                                                                                    ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCt                                                             <  2:499199/138‑1 (MQ=255)
                                                                                    ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCt                                                             >  1:499199/1‑138 (MQ=255)
                                                                                    ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATag                                                 >  2:677925/1‑150 (MQ=255)
                                                                                          gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCa                                                                                                  <  1:640328/95‑1 (MQ=255)
                                                                                          gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCa                                                                                                  >  2:640328/1‑95 (MQ=255)
                                                                                          gagaCTAGGACATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTGGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCAAACAGCATTTCTTCACGACCATAGAGCAtt                                           <  2:504435/150‑1 (MQ=25)
                                                                                              gTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGaa                                       <  1:677925/150‑1 (MQ=255)
                                                                                                  gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAaccac                                   >  1:58657/1‑150 (MQ=255)
                                                                                                         tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGt                                                                                                                      <  2:57347/60‑1 (MQ=38)
                                                                                                         tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGt                                                                                                                      >  1:57347/1‑60 (MQ=38)
                                                                                                         tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcc                            >  2:320582/1‑150 (MQ=255)
                                                                                                               ttCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatccc                      >  1:346044/1‑150 (MQ=255)
                                                                                                                          aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGTGCATTGGAaccac                                   >  2:637300/1‑127 (MQ=255)
                                                                                                                          aaaaCCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGTGCATTGGAaccac                                   <  1:637300/127‑1 (MQ=255)
                                                                                                                              ccccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGa        >  2:841976/1‑150 (MQ=37)
                                                                                                                               cccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                   >  1:97849/1‑138 (MQ=37)
                                                                                                                               cccTACCTGCATGTGCAGATAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa                   <  2:97849/138‑1 (MQ=37)
                                                                                                                                    tCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt  <  1:320582/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
CGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAT  >  NC_002947/181674‑181954

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