Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,931,799 T→C N71N (AAT→AAC rbsR → ribose operon transcriptional repressor RbsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,931,7990TC98.2% 186.2 / ‑3.8 55N71N (AAT→AACrbsRribose operon transcriptional repressor RbsR
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (23/31);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

GCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAA  >  NZ_CP009273/3931663‑3931926
                                                                                                                                        |                                                                                                                               
gCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAAccc                                                                                                                               >  1:1247894/1‑139 (MQ=255)
    cGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCttt                                                                                                                            <  1:109529/138‑1 (MQ=255)
    cGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCttt                                                                                                                            <  2:1247894/138‑1 (MQ=255)
    cGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAAATCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCttt                                                                                                                            <  1:719364/138‑1 (MQ=255)
        aGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTAtt                                                                                                                       <  1:473257/139‑1 (MQ=255)
        aGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTAtt                                                                                                                       <  2:537173/139‑1 (MQ=255)
        aGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTAtt                                                                                                                       <  1:147386/139‑1 (MQ=255)
          tGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTAtt                                                                                                                       >  1:521378/1‑137 (MQ=255)
                 aTTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTgg                                                                                                              >  1:277098/1‑139 (MQ=255)
                 aTTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTgg                                                                                                              >  2:1260802/1‑139 (MQ=255)
                    aCCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGc                                                                                                           >  1:975746/1‑139 (MQ=255)
                     ccGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCg                                                                                                          <  1:662787/139‑1 (MQ=255)
                          aaaGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCg                                                                                                     <  1:956620/139‑1 (MQ=255)
                                gAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCttt                                                                                                                            >  1:190335/1‑110 (MQ=255)
                                gAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCttt                                                                                                                            <  2:190335/110‑1 (MQ=255)
                                gAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGaa                                                                                                <  2:1265338/138‑1 (MQ=255)
                                        gATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAgctgc                                                                                       >  2:848669/1‑139 (MQ=255)
                                              aGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGaa                                                                                 <  1:1064394/139‑1 (MQ=255)
                                                aaCTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcg                                                                               <  1:993791/139‑1 (MQ=255)
                                                aaCTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcg                                                                               <  2:491068/139‑1 (MQ=255)
                                                aaCTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcg                                                                               <  1:1016557/139‑1 (MQ=255)
                                                aaCTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcg                                                                               <  2:983309/139‑1 (MQ=255)
                                                aaCTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcg                                                                               <  1:848669/139‑1 (MQ=255)
                                                aaCTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcg                                                                               <  1:17418/139‑1 (MQ=255)
                                                aaCTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcg                                                                               <  1:346384/139‑1 (MQ=255)
                                                           gcgcCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGtc                                                                    <  1:12386/139‑1 (MQ=255)
                                                           gcgcCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGtc                                                                    <  1:322503/139‑1 (MQ=255)
                                                            cgcCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGtct                                                                   >  1:297402/1‑139 (MQ=255)
                                                               cATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGt                                                                >  1:598899/1‑139 (MQ=255)
                                                                aTCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTGTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCCTCGAACGCGGTTATAGTCTCGTc                                                               <  2:910825/139‑1 (MQ=255)
                                                                         ggCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCaa                                                       >  1:854360/1‑138 (MQ=255)
                                                                         ggCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCaa                                                       >  1:244991/1‑138 (MQ=255)
                                                                             cgTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATAcc                                                   >  1:773329/1‑138 (MQ=255)
                                                                               tAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGaa                                                                                                <  2:600797/91‑1 (MQ=255)
                                                                               tAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGaa                                                                                                >  1:600797/1‑91 (MQ=255)
                                                                               tAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTta                                                                         >  1:1036866/1‑114 (MQ=255)
                                                                               tAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTta                                                                         <  2:1036866/114‑1 (MQ=255)
                                                                                                aaaCACATACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGAt                               <  1:145295/139‑1 (MQ=255)
                                                                                                      aTACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcg                                                                               >  1:143933/1‑85 (MQ=255)
                                                                                                      aTACCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcg                                                                               <  2:143933/85‑1 (MQ=255)
                                                                                                        aCCATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGc                        >  1:1201656/1‑138 (MQ=255)
                                                                                                          cATTGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAAt                     <  1:785846/139‑1 (MQ=255)
                                                                                                            ttGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcg                                                                               <  2:491321/79‑1 (MQ=255)
                                                                                                            ttGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcg                                                                               >  1:491321/1‑79 (MQ=255)
                                                                                                            ttGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCt                   >  1:1107737/1‑139 (MQ=255)
                                                                                                            ttGGCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCt                   >  2:888880/1‑139 (MQ=255)
                                                                                                              ggCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTg                                                                                                        >  1:839106/1‑52 (MQ=255)
                                                                                                              ggCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTg                                                                                                        <  2:839106/52‑1 (MQ=255)
                                                                                                              ggCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTg                                                                                                        <  2:876805/52‑1 (MQ=255)
                                                                                                              ggCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTg                                                                                                        >  1:876805/1‑52 (MQ=255)
                                                                                                               gCATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGa                <  1:918396/139‑1 (MQ=255)
                                                                                                                cATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCg                         <  1:973323/129‑1 (MQ=255)
                                                                                                                cATGTTGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCg                         >  2:973323/1‑129 (MQ=255)
                                                                                                                    ttGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGc           >  2:1007643/1‑139 (MQ=255)
                                                                                                                    ttGATCACTGCCAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGc           >  2:951779/1‑139 (MQ=255)
                                                                                                                              ccAGTACCAACCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCaaa  >  2:299015/1‑138 (MQ=255)
                                                                                                                                        |                                                                                                                               
GCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAA  >  NZ_CP009273/3931663‑3931926

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