Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,604,910 C→T 8.0% A67A (GCC→GCT U712_08250 → Hypothetical Protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,604,9100CT 8.0% 192.9 / 5.1 62A67A (GCC→GCTU712_08250Hypothetical Protein
Reads supporting (aligned to +/- strand):  ref base C (25/32);  new base T (3/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

GTCAGCCCGGAGTCTGCGCCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAGGTCCCTCAGGAGTTGGTAAAGGAACGGT  >  CP006881/1604776‑1605048
                                                                                                                                      |                                                                                                                                          
gTCAGCCCGGAGTCTGCGCCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTc                                                                                                                                         <  6:131495/138‑1 (MQ=255)
gTCAGCCCGGAGTCTGCGCCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTc                                                                                                                                         <  6:299447/138‑1 (MQ=255)
    gCCCGGAGTCTGCGCCACTCAAACGGCTGATTCAGGCTGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGcc                                                                                                                                    <  8:342112/139‑1 (MQ=255)
       cGGAGTCTGAGCCAATCAAACGGATGATTCCGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCGga                                                                                                                                 <  5:334176/139‑1 (MQ=255)
         gAGTCTGCGCCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTgaga                                                                                                                               <  1:235017/139‑1 (MQ=255)
         gAGTCTGCGCCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTgaga                                                                                                                               <  3:207700/139‑1 (MQ=255)
         gAGTCTGCGCCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTgaga                                                                                                                               <  2:336644/139‑1 (MQ=255)
             cTGCGCCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGt                                                                                                                           <  5:375558/139‑1 (MQ=255)
               gcgcCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTg                                                                                                                         >  4:16147/1‑139 (MQ=255)
               gcgcCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTg                                                                                                                         >  3:298202/1‑139 (MQ=255)
                cgcCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGt                                                                                                                           <  3:7131/136‑1 (MQ=255)
                cgcCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGc                                                                                                                        >  8:24015/1‑139 (MQ=255)
                 gcCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTCTCTTATCTGCCGTCCAGCCTGAGACAGTTGca                                                                                                                       <  7:90110/139‑1 (MQ=255)
                   cAATCAAACGGATGATTCAGGATGCAAGAGACAGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGcaca                                                                                                                     <  2:50944/139‑1 (MQ=255)
                    aaTCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTc                                                                                                                                         >  5:333143/1‑118 (MQ=255)
                    aaTCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTc                                                                                                                                         <  6:333143/118‑1 (MQ=255)
                      tCAAACGGATGATTCAGGATGCAAGAGCCCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAg                                                                                                                  >  1:29566/1‑139 (MQ=255)
                                       gATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAaagaag                                                                                                 <  3:16147/139‑1 (MQ=255)
                                        aTGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAaagaag                                                                                                 >  7:256820/1‑138 (MQ=255)
                                              gagaCCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTAt                                                                                          <  7:24015/139‑1 (MQ=255)
                                                   cgctGAATTCAAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATg                                                                                     <  6:41695/139‑1 (MQ=255)
                                                   cgcgGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGAt                                                                                      >  4:109847/1‑138 (MQ=255)
                                                        aaTGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAggg                                                                                <  7:28044/139‑1 (MQ=255)
                                                         aTGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAggg                                                                                <  1:43241/138‑1 (MQ=255)
                                                          tGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCa                                                                              <  6:137703/139‑1 (MQ=255)
                                                          tGCTAATTGACGCTACATACGGAAGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCa                                                                              <  3:23254/139‑1 (MQ=255)
                                                                    cGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTg                                                                    >  6:37061/1‑139 (MQ=255)
                                                                      cTACATACGGACGAAGAACCCGAGCAGTTGTCGTCATGGATAGTGATCACATTATCTAATCTGCCGGCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGcc                                                                  >  1:284621/1‑139 (MQ=255)
                                                                              ggACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATga                                                            >  3:264382/1‑137 (MQ=255)
                                                                              ggACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATga                                                            >  1:178527/1‑137 (MQ=255)
                                                                                ccgcagaaCCCGTGCAGTTGTCGTCATGGAGAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAATAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATgaaaga                                                        <  4:12071/135‑1 (MQ=255)
                                                                                 cGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGcc                                                                                                                                    <  1:309769/62‑1 (MQ=255)
                                                                                 cGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGcc                                                                                                                                    >  2:309769/1‑62 (MQ=255)
                                                                                 cGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGcaca                                                                                                                     <  1:152634/77‑1 (MQ=255)
                                                                                 cGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGcaca                                                                                                                     >  2:152634/1‑77 (MQ=255)
                                                                                 cGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATgaaaga                                                        >  8:39075/1‑138 (MQ=255)
                                                                                   aagaaCCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATgaaaga                                                        >  6:161411/1‑136 (MQ=255)
                                                                                   aagaaCCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGaaagaaag                                                     <  5:37061/139‑1 (MQ=255)
                                                                                    agaaCCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAaagaaaga                                                    <  5:112659/139‑1 (MQ=255)
                                                                                     gaaCCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTAt                                                                                          <  5:141567/100‑1 (MQ=255)
                                                                                     gaaCCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTAt                                                                                          >  6:141567/1‑100 (MQ=255)
                                                                                            tGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGttat                                             >  2:28628/1‑138 (MQ=255)
                                                                                            tGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGttat                                             >  7:353459/1‑138 (MQ=255)
                                                                                                gTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATc                                        <  8:353459/139‑1 (MQ=255)
                                                                                                 ttGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTga                                                                                                                                 >  3:47489/1‑49 (MQ=255)
                                                                                                 ttGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTga                                                                                                                                 <  4:47489/49‑1 (MQ=255)
                                                                                                   gtcgtcATGGATAGTGATCACATTATCTTATCTGCTGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATCAAAGAAAGAGGGTTATTAATCGtt                                     >  1:100883/1‑139 (MQ=255)
                                                                                                         aTGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCa                               <  1:28628/139‑1 (MQ=255)
                                                                                                                 tGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAaagaaaga                                                    <  4:113760/110‑1 (MQ=255)
                                                                                                                 tGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAaagaaaga                                                    >  3:113760/1‑110 (MQ=255)
                                                                                                                 tGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAaagaaaga                                                    <  2:341181/110‑1 (MQ=255)
                                                                                                                 tGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAaagaaaga                                                    >  1:341181/1‑110 (MQ=255)
                                                                                                                 tGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAGGTCCCTc                       >  6:387697/1‑139 (MQ=255)
                                                                                                                      acaTTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGGTATTAATCGTTCTCTCAGGTCCCTCAGGa                   <  8:256820/138‑1 (MQ=255)
                                                                                                                           atcttatctGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAGGTCCCTCAGGAGTTgg              >  3:338134/1‑138 (MQ=255)
                                                                                                                           atcttatctGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAGGTCCCTCAGGAGTTgg              >  7:334391/1‑138 (MQ=255)
                                                                                                                                atctGCTGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAgg                             >  7:256181/1‑118 (MQ=255)
                                                                                                                                atctGCTGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAgg                             >  1:178381/1‑118 (MQ=255)
                                                                                                                                atctGCTGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAgg                             <  8:256181/118‑1 (MQ=255)
                                                                                                                                atctGCTGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAgg                             <  2:178381/118‑1 (MQ=255)
                                                                                                                                atctGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAGGTCCCTCAGGAGTTGGTaa           <  4:264382/136‑1 (MQ=255)
                                                                                                                                  ctGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAGGTCCCTCAGGAGTTGGTAAAGGaa      >  1:240717/1‑139 (MQ=255)
                                                                                                                                      cGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAACGAGGGTTATTAATCGTTCTCTCAgg                             <  4:359725/112‑1 (MQ=255)
                                                                                                                                      cGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAgg                             >  3:359725/1‑112 (MQ=255)
                                                                                                                                      cGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAGGTCCCTCAGGAGTTGGTAAAGGAACGGt  >  4:247917/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                          
GTCAGCCCGGAGTCTGCGCCAATCAAACGGATGATTCAGGATGCAAGAGACCGCGGAATGCTAATTGACGCTACATACGGACGAAGAACCCGTGCAGTTGTCGTCATGGATAGTGATCACATTATCTTATCTGCCGTCCAGCCTGAGACAGTTGCACACAGACTTTCTGTTAAAGAAGAAATTATGGATGAAGGGCAGGGGTAATTGCCGCATGAAAGAAAGAGGGTTATTAATCGTTCTCTCAGGTCCCTCAGGAGTTGGTAAAGGAACGGT  >  CP006881/1604776‑1605048

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: