Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,071,029 G→A 5.9% T246M (ACG→ATG)  rutD ← putative reactive intermediate detoxifying aminoacrylate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,071,0290GA 5.9% 215.7 / 2.5 68T246M (ACG→ATG) rutDputative reactive intermediate detoxifying aminoacrylate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (34/30);  new base A (2/2);  total (36/32)
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.55e-01

CCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGC  >  NC_000913/1070884‑1071173
                                                                                                                                                 |                                                                                                                                                
ccGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGttt                                                                                                                                               >  4:105039/1‑149 (MQ=255)
ccGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGttt                                                                                                                                               <  1:99588/149‑1 (MQ=255)
 cGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTc                                                                                                                                              <  3:105039/149‑1 (MQ=255)
 cGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTc                                                                                                                                              >  4:380482/1‑149 (MQ=255)
   tcAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCgg                                                                                                                                            <  7:145576/148‑1 (MQ=255)
      gCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGAt                                                                                                                                         >  1:286323/1‑149 (MQ=255)
             tgAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACATTTCGGGATCGGTCAc                                                                                                                                  >  1:22767/1‑149 (MQ=255)
             tgAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACATTTCGGGATCGGTCAc                                                                                                                                  >  1:29059/1‑149 (MQ=255)
                    gtgCATCGTTGAACAGGGTGCTAATCGCACCTGGGATAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGTTGATTTAACAGGCTGGCAAGCCAGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAg                                                                                                                           <  3:244310/147‑1 (MQ=255)
                    gggCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAg                                                                                                                           <  7:54870/149‑1 (MQ=255)
                    gggCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAg                                                                                                                           >  4:245067/1‑149 (MQ=255)
                      gCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGc                                                                                                                         <  3:179948/149‑1 (MQ=255)
                               aaCAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGc                                                                                                                 <  3:380482/148‑1 (MQ=255)
                               aaCAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGc                                                                                                                 <  2:286323/148‑1 (MQ=255)
                               aaCAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGc                                                                                                                 >  1:105199/1‑148 (MQ=255)
                                 cAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCAt                                                                                                              <  4:197154/149‑1 (MQ=255)
                                  aGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATa                                                                                                             <  8:176032/149‑1 (MQ=255)
                                     gTGCTAAGCGCACCTGGGCTAACGGGGTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAggg                                                                                                          <  3:244418/149‑1 (MQ=255)
                                            gcgcACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTAGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAg                                                                                                            <  1:369925/140‑1 (MQ=255)
                                            gcgcACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTAGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAg                                                                                                            >  2:369925/1‑140 (MQ=255)
                                                     ggCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCt                                                                                            <  3:245067/147‑1 (MQ=255)
                                                     ggCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCt                                                                                            <  2:105199/147‑1 (MQ=255)
                                                      gCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGc                                                                                         >  7:183970/1‑149 (MQ=255)
                                                        tAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCg                                                                                                                                             <  3:288744/95‑1 (MQ=255)
                                                        tAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCg                                                                                                                                             >  4:288744/1‑95 (MQ=255)
                                                        tAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCCCCATTTTCTGGCTa                                                                                       >  2:322930/1‑149 (MQ=255)
                                                         aaCGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTAt                                                                                      >  2:34620/1‑149 (MQ=255)
                                                          aCGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATc                                                                                     <  5:93115/149‑1 (MQ=255)
                                                           cGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCg                                                                                    <  6:268084/149‑1 (MQ=255)
                                                              cTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGACCGTTGAGTAACAGAGCATTAAACGTTTGGGTTTCGGTCACGTTGCAGGCGTTTCCGCCATAGGGCATAACCATTTTCTGGCTATCggg                                                                                  >  1:308301/1‑148 (MQ=255)
                                                               ttCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGAGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGAACGTTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCAttt                                                                                               >  3:316893/1‑134 (MQ=255)
                                                                tCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAGCGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAg                                                                               >  5:78415/1‑149 (MQ=255)
                                                                 cGTTCATAGCAATTCCTTACTGGGCGGCTTCTCGGTGATGTAACAGGCTTGCAATCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGGAGGCGTGTCCGCCATAGGGGATCACCATTTTCTGGCTATCGGGCa                                                                                <  4:37817/147‑1 (MQ=255)
                                                                   ttCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCAcc                                                                                                   <  6:331700/126‑1 (MQ=255)
                                                                   ttCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCAcc                                                                                                   >  5:331700/1‑126 (MQ=255)
                                                                     cATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTACACGTTTCGGGCTCGGGCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGggcg                                                                           <  5:276901/148‑1 (MQ=255)
                                                                       tAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCa                                                                        >  8:18697/1‑149 (MQ=255)
                                                                       tAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCa                                                                        >  5:346836/1‑149 (MQ=255)
                                                                        aGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGggcg                                                                           <  6:360640/145‑1 (MQ=255)
                                                                        aGCAAATCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGggcg                                                                           >  5:360640/1‑145 (MQ=255)
                                                                          cAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGTATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATg                                                                      >  5:149324/1‑148 (MQ=255)
                                                                           aaTTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATg                                                                      >  7:219727/1‑147 (MQ=255)
                                                                             ttCCTTTCAGGGCGGCTTCACGGTGATGTAACAGGCTGGCACGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACCGTGCAGGCGTGACCGCCATATGGCATAACCATTTTCTGGCTATCGCTCAGCGCGGCATGAAGt                                                                  >  5:9119/1‑149 (MQ=255)
                                                                                 ttACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGGAGTTCAc                                                              >  2:219104/1‑149 (MQ=255)
                                                                                      gggCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACATTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCTCTGGaa                                                         <  2:29059/149‑1 (MQ=255)
                                                                                      gggCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACATTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCGCTGGaa                                                         <  2:22767/149‑1 (MQ=255)
                                                                                       ggcggcTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGACCATTAAAAGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAAc                                                        <  6:295392/149‑1 (MQ=255)
                                                                                          ggcTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATg                                                     <  7:18697/149‑1 (MQ=255)
                                                                                             ttCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTc                                                                                                                    <  5:54589/83‑1 (MQ=255)
                                                                                             ttCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTc                                                                                                                    >  6:54589/1‑83 (MQ=255)
                                                                                                       tGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACcag                                        >  2:71394/1‑149 (MQ=255)
                                                                                                         tAACAGGCTGGCAAGCCCGTTGCGTAACAGAGCATTAACCGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACcagca                                      <  6:305682/149‑1 (MQ=255)
                                                                                                              ggCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATg                                                                      >  6:7549/1‑112 (MQ=255)
                                                                                                              ggCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATg                                                                      <  5:7549/112‑1 (MQ=255)
                                                                                                               gCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGatcat                                >  2:319747/1‑149 (MQ=255)
                                                                                                                 tGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGAtcatca                              >  2:282902/1‑149 (MQ=255)
                                                                                                                 tGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGAtcatca                              >  1:52541/1‑149 (MQ=255)
                                                                                                                       tccaggtGAGTAAACGTGGATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCAGCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGgcg                        <  2:308301/143‑1 (MQ=255)
                                                                                                                              gAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgatga                 >  8:393939/1‑149 (MQ=255)
                                                                                                                               aGTAACAGAGCAGTAAACGGTTCGTGTTCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAgatgat                <  8:183970/149‑1 (MQ=255)
                                                                                                                                 tAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGAttt              >  7:135074/1‑149 (MQ=255)
                                                                                                                                       agCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCg        >  4:331602/1‑149 (MQ=255)
                                                                                                                                         cATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGaa                                                                    <  6:131995/87‑1 (MQ=255)
                                                                                                                                         cATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGaa                                                                    >  5:131995/1‑87 (MQ=255)
                                                                                                                                           ttAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCggg      >  7:208412/1‑147 (MQ=255)
                                                                                                                                           ttAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCggg      >  8:256344/1‑147 (MQ=255)
                                                                                                                                           ttAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCa    >  8:350469/1‑149 (MQ=255)
                                                                                                                                              aaCGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCCCCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGc  <  1:322930/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCGCCAGCCGTTGTGAGTGCGGGCATCGGTGAACAGGGTGCTAAGCGCACCTGGGCTAACGGCTTCGTTCATAGCAATTCCTTACAGGGCGGCTTCACGGTGATGTAACAGGCTGGCAAGCCCGTTGAGTAACAGAGCATTAAACGTTTCGGGATCGGTCACGTTGCAGGCGTGTCCGCCATAGGGCATCACCATTTTCTGGCTATCGGGCAGGGCGGCATGAAGTTCACTGGAACATGCTGTTGGCACCAGCAGATCATCACTGGCGCAGATGATTTGCACCGGGCAGC  >  NC_000913/1070884‑1071173

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: