Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 181,805 T→C 13.2% intergenic (+19/‑87) 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,8050TC13.2% 102.6 / 6.5 44intergenic (+19/‑87)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (21/17);  new base C (2/4);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

AATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAT  >  NC_002947/181664‑181954
                                                                                                                                              |                                                                                                                                                     
aaTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGcatgc                                                                                                                                               <  2:25761/150‑4 (MQ=11)
          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                                                                                                          <  2:556417/148‑1 (MQ=40)
                                       ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACaga                                                                                                          >  1:556417/1‑148 (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                                                                                                                                       >  1:802764/1‑74 (MQ=35)
                                                                                    tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAg                                                                                                                                       <  2:802764/74‑1 (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                                                                                                                      >  1:57347/1‑60 (MQ=38)
                                                                                                                   tCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGt                                                                                                                      <  2:57347/60‑1 (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)
                                                                                                                                              |                                                                                                                                                     
AATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAT  >  NC_002947/181664‑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: