Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,333,425 Δ1 bp coding (526/669 nt) basR ← DNA‑binding transcriptional dual regulator BasR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,333,4240T.96.9% 122.1 / ‑2.0 32coding (527/669 nt)basRDNA‑binding transcriptional dual regulator BasR
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (18/13);  total (19/13)
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.82e-01

GCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGC  >  NC_000913/4333345‑4333498
                                                                               |                                                                          
gCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGtg                                                                           >  2:890846/1‑80 (MQ=255)
 cccACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGtgt                                                                          <  2:890807/81‑3 (MQ=255)
 cccACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGtg                                                                           >  1:890844/1‑79 (MQ=255)
    aCTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGtgta                                                                         >  2:890845/1‑76 (MQ=255)
    aCTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTaga                                                                      <  2:890808/81‑1 (MQ=255)
       ttGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGaa                                                                   >  2:890863/1‑81 (MQ=255)
         gTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAAt                                                                  <  1:890809/80‑1 (MQ=255)
          tCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAAtt                                                                 <  2:890810/80‑1 (MQ=255)
             cgcAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTccc                                                             >  1:890853/1‑81 (MQ=255)
                aGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGAt                                                          >  1:890811/1‑81 (MQ=255)
                  aTTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATg                                                         >  1:890850/1‑80 (MQ=255)
                  aTTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGc                                                        <  2:890812/81‑1 (MQ=255)
                  aTTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGc                                                        <  2:890813/81‑1 (MQ=255)
                   ttGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCa                                                       >  1:890905/1‑81 (MQ=255)
                    tgtgGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCAc                                                      >  1:890856/1‑81 (MQ=255)
                      tgGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCAcc                                                     <  1:890814/80‑1 (MQ=255)
                      tgGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCg                                                    >  1:890878/1‑81 (MQ=255)
                               aCTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGActgcct                                            >  2:890865/1‑81 (MQ=255)
                                   ccAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCttt                                       >  1:890869/1‑81 (MQ=255)
                                   ccAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCttt                                       >  2:890855/1‑81 (MQ=255)
                                    cAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTg                                      <  1:890815/81‑1 (MQ=255)
                                      gggTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAg                                    >  2:890876/1‑81 (MQ=255)
                                         tgtTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCAt                                 >  1:890879/1‑81 (MQ=255)
                                           ttGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCAtt                                >  2:890851/1‑80 (MQ=255)
                                           ttGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTa                               <  2:890816/81‑1 (MQ=255)
                                               tCGAGGGTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCg                           >  2:890864/1‑81 (MQ=255)
                                                   gggTTCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGAc                       <  2:890817/81‑1 (MQ=255)
                                                       tCATTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGAcag                     >  1:890858/1‑79 (MQ=255)
                                                         aTTGTCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGAcagca                   >  2:890877/1‑79 (MQ=255)
                                                             tCCCAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCat             <  2:890818/81‑1 (MQ=255)
                                                               ccAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCatat           >  2:890870/1‑81 (MQ=255)
                                                                cAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATAtt          <  2:890820/81‑1 (MQ=255)
                                                                cAGTTATAGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATAtt          <  1:890819/81‑1 (MQ=255)
                                                                    tataGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCttt      >  2:890919/1‑81 (MQ=255)
                                                                       aGATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTggg   >  1:890872/1‑81 (MQ=255)
                                                                        gATGTCGTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGc  <  2:890821/81‑1 (MQ=255)
                                                                               |                                                                          
GCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGC  >  NC_000913/4333345‑4333498

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