Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,334,340 A→G 6.2% R145R (CGT→CGC)  PP_3806 ← isochorismatase family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,334,3400AG 6.2% 238.7 / 2.5 79R145R (CGT→CGC) PP_3806isochorismatase family hydrolase
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base G (3/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCGAAATGGCCAGTAACCGGCCGTTCAGGGCTGTTCACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGTTTCCATGCCGCAGACGATCAC  >  NC_002947/4334196‑4334486
                                                                                                                                                |                                                                                                                                                  
tCACCGAAATGGCCAGTAACCGGCCGTTCAGGGCTGTTCACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCa                                                                                                                                                >  4:179958/1‑149 (MQ=255)
  aCCGAAATGGCCAGTAACCGGCCGTTCAGGGCTGTTCACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCAt                                                                                                                                               <  5:134020/148‑1 (MQ=255)
   ccGAAATGGCCAGTAACCGGCCGTTCAGGGCTGTTCACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGTGCCGCCGGCGCGCag                                                                                                                                               <  8:58726/147‑2 (MQ=255)
   ccGAAATGGCCAGTAACCGGCCGTTCAGGGCTGTTCACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGTGCCGCCGGCGCGCag                                                                                                                                               >  7:58726/1‑146 (MQ=255)
     gAAATGGCCAGTAACCGGCCGTTCAGGGCTGTTCACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATgcgc                                                                                                                                           <  5:315104/149‑1 (MQ=255)
       aaTGGCCAGTAACCGGCCGTTCAGGGCTGTTCACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATgcgcgc                                                                                                                                         >  2:15178/1‑149 (MQ=255)
                              gggCTGTTCACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCg                                                                                                                  >  8:64490/1‑149 (MQ=255)
                              gggCTGTTCACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTACTg                                                                                                                      >  8:247670/1‑145 (MQ=255)
                              gggCTGTTCACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTACTg                                                                                                                      <  7:247670/145‑1 (MQ=255)
                                gCTGTTCACCGATCAGGTAGGATTTGCTGATGTGGCGCAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCggg                                                                                                                <  6:262551/149‑1 (MQ=255)
                                      cACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCg                                                                                                           <  6:31043/148‑1 (MQ=255)
                                      cACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCg                                                                                                           <  3:179958/148‑1 (MQ=255)
                                       aCCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGc                                                                                                         >  3:115748/1‑149 (MQ=255)
                                               ggtaggtTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAg                                                                                                 >  1:82567/1‑149 (MQ=255)
                                               ggtaggtTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAg                                                                                                 <  4:115748/149‑1 (MQ=255)
                                               ggtaggtTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAg                                                                                                 <  3:161865/149‑1 (MQ=255)
                                                 taggtTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGc                                                                                               <  6:208873/149‑1 (MQ=255)
                                                 taggtTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGc                                                                                               <  1:15178/149‑1 (MQ=255)
                                                 taggtTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACAATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAgg                                                                                                >  6:236800/1‑148 (MQ=255)
                                                     ttttGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATcc                                                                                           >  7:246181/1‑149 (MQ=255)
                                                     ttttGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATcc                                                                                           >  7:300410/1‑149 (MQ=255)
                                                             aTGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGAcc                                                                                     <  5:221341/147‑1 (MQ=255)
                                                               gtgGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACg                                                                                  >  5:24296/1‑148 (MQ=255)
                                                                tgGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGaa                                                                                <  2:82567/149‑1 (MQ=255)
                                                                  gCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAaca                                                                              >  5:72837/1‑149 (MQ=255)
                                                                    ggAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACAct                                                                            <  8:246181/149‑1 (MQ=255)
                                                                       aCTGCGGGTGCCCGGCACTGCCCATCAATTCGACCAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACActtgc                                                                         <  2:4693/149‑1 (MQ=255)
                                                                       aCTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAacac                                                                             >  8:212550/1‑145 (MQ=255)
                                                                       aCTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAacac                                                                             <  7:212550/145‑1 (MQ=255)
                                                                       aCTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACActtgc                                                                         <  5:58407/149‑1 (MQ=255)
                                                                              gTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGcc                                                                    >  1:151491/1‑147 (MQ=255)
                                                                                gCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAgc                                                                >  5:123837/1‑149 (MQ=255)
                                                                                  ccGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAgcgc                                                              >  7:32134/1‑149 (MQ=255)
                                                                                  ccGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAgcgc                                                              >  2:246998/1‑149 (MQ=255)
                                                                                  ccGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAgcgc                                                              >  5:186714/1‑149 (MQ=255)
                                                                                     gCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCt                                                                                                                       >  6:115771/1‑89 (MQ=255)
                                                                                     gCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCt                                                                                                                       <  5:115771/89‑1 (MQ=255)
                                                                                         tGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAgg                                                       >  8:125488/1‑149 (MQ=255)
                                                                                         tGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAgg                                                       >  8:304514/1‑149 (MQ=255)
                                                                                           cccATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGc                                                      <  7:125488/148‑1 (MQ=255)
                                                                                              aTCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAg                                                  >  8:175720/1‑149 (MQ=255)
                                                                                              aTCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAg                                                  >  3:19147/1‑149 (MQ=255)
                                                                                              aTCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAg                                                  >  4:293516/1‑149 (MQ=255)
                                                                                                  aTTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACg                                              >  3:199294/1‑149 (MQ=255)
                                                                                                   ttCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCtt                                                                                                                          <  6:195511/72‑1 (MQ=255)
                                                                                                   ttCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCtt                                                                                                                          >  5:195511/1‑72 (MQ=255)
                                                                                                   ttCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCg                                                    >  7:281807/1‑142 (MQ=255)
                                                                                                   ttCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCg                                                    <  8:281807/142‑1 (MQ=255)
                                                                                                       aacaacACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTg                                         <  6:72837/149‑1 (MQ=255)
                                                                                                       aacaacACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTg                                         <  3:293516/149‑1 (MQ=255)
                                                                                                       aacaacACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTg                                         >  2:177264/1‑149 (MQ=255)
                                                                                                        acaacaCCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTgc                                        <  1:177264/149‑1 (MQ=255)
                                                                                                         caacaCCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTgca                                       <  7:64490/149‑1 (MQ=255)
                                                                                                           acacCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTgcagc                                     <  4:19147/149‑1 (MQ=255)
                                                                                                            cacCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTgcagca                                    <  1:246998/149‑1 (MQ=255)
                                                                                                            cacCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTgcagca                                    >  2:2766/1‑149 (MQ=255)
                                                                                                               cATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGcc                                                             >  6:137205/1‑121 (MQ=255)
                                                                                                               cATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGcc                                                             <  5:137205/121‑1 (MQ=255)
                                                                                                                   tcGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCACGCATCCTCGACCACGAACACTTGCTTGCCCAgc                                                                >  8:14190/1‑114 (MQ=255)
                                                                                                                     gcgGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACAt                           <  6:24296/149‑1 (MQ=255)
                                                                                                                     gcgGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACAt                           <  6:123837/149‑1 (MQ=255)
                                                                                                                     gcgGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACAt                           >  4:18503/1‑149 (MQ=255)
                                                                                                                       gggTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCg                                                                                                    >  1:270484/1‑74 (MQ=255)
                                                                                                                       gggTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCg                                                                                                    <  2:270484/74‑1 (MQ=255)
                                                                                                                       gggTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATgg                         <  8:32134/149‑1 (MQ=255)
                                                                                                                            aCTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGttt                     <  5:236800/148‑1 (MQ=255)
                                                                                                                            aCTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGTTTc                    >  7:293145/1‑149 (MQ=255)
                                                                                                                              tACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGTTTCCa                  >  6:232827/1‑149 (MQ=255)
                                                                                                                                 atGCGTGCCGCCGGCGCGCATGCGCGCCAGGCATGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATTGGTTTCCATGc               >  6:263848/2‑149 (MQ=255)
                                                                                                                                 atGCGTGCCGCCGGCGCGCATGCGAGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCCGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGTTTCCATGc               >  1:270337/2‑149 (MQ=255)
                                                                                                                                   gcgcgcCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACActtg                                                                          >  4:8307/1‑88 (MQ=255)
                                                                                                                                   gcgcgcCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACActtg                                                                          <  3:8307/88‑1 (MQ=255)
                                                                                                                                   gcgcgcCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGTTTCCATg                <  1:48017/146‑1 (MQ=255)
                                                                                                                                   gcgcgcCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGTTTCCATg                >  2:48017/1‑146 (MQ=255)
                                                                                                                                    cgtgccgccGGCGCGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGTTTCCATGCCg             <  5:263848/145‑1 (MQ=255)
                                                                                                                                            cggcACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGTTTCCATGCCGCAGACGAt     >  2:193199/1‑148 (MQ=255)
                                                                                                                                            cggcACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGTTTCCATGCCGCAGACGAt     >  2:188557/1‑148 (MQ=255)
                                                                                                                                              gcACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGTTTCCATGCCGCAGACGATCAc  <  8:293145/149‑1 (MQ=255)
                                                                                                                                              gcACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGTTTCCATGCCGCAGACGATCAc  <  3:18503/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TCACCGAAATGGCCAGTAACCGGCCGTTCAGGGCTGTTCACCGATCAGGTAGGTTTTGCTGATGTGGCGGAACTGCGGGTGCCCGGCACTGCCCATCAATTCGAACAACACCATCTCGCGGGTCACTACCTGCGCGCCGGCGGCACGCATGCGCGCCAGGCCTGCGGCCTTGCTGGCCGGGGTGCGGCTGTCGCAGGCATCCTCGACCACGAACACTTGCTTGCCCAGCGCCAACAGGCCGAGCACGGTCTGCAGCACGCAAACATGGGTTTCCATGCCGCAGACGATCAC  >  NC_002947/4334196‑4334486

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