Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,136,489 T→A 5.8% F233Y (TTC→TAC)  PP_3640 → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,136,4890TA 5.8% 214.9 / 2.2 69F233Y (TTC→TAC) PP_3640AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (28/37);  new base A (2/2);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

GCAGGGCAATCGCTGGGCGGGCTGGTCGCCCTCGCCCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAGGCGCGCAGGTGGGGTATCAGTCGGTGGCCGCATTCA  >  NC_002947/4136344‑4136635
                                                                                                                                                 |                                                                                                                                                  
gCAGGGCAATCGCTGGGCGGGCTGGTCGCCCTCGCCCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCaa                                                                                                                                                 <  2:154548/149‑1 (MQ=255)
gCAGGGCAATCGCGGGGCGGGCTGGTCGCCCTCGCGCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCaa                                                                                                                                                 <  4:49897/149‑1 (MQ=255)
     gCAATCGCTGGGCGGGCTGGTCGCCCTCGCCCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGtt                                                                                                                                            <  3:311573/149‑1 (MQ=255)
       aaTCGCTGGGCGGGCTGGTCGCCCTCGCCCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGttt                                                                                                                                           >  7:7051/1‑148 (MQ=255)
         tCGCTGGGCGGGCTGGTCGCCCTCGCCCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCaa                                                                                                                                                 <  5:139581/140‑1 (MQ=255)
         tCGCTGGGCGGGCTGGTCGCCCTCGCCCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCaa                                                                                                                                                 >  6:139581/1‑140 (MQ=255)
                gcgggcTGGTCGCCCTCGCCCGCCAGCCGGGATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTg                                                                                                                                 <  1:234466/149‑1 (MQ=255)
                   ggcTGGTCGCCCTCGCCCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGcc                                                                                                                              <  8:7051/149‑1 (MQ=255)
                     cTGGTCGCCCTCGCCCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCgg                                                                                                                            <  7:38880/149‑1 (MQ=255)
                      tGGTCGCCCTCGCCCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCggg                                                                                                                           >  3:87775/1‑149 (MQ=255)
                       ggTCGCCCTCGCCCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGc                                                                                                                          <  2:149735/149‑1 (MQ=255)
                        gTCGCCCTCGCCCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGc                                                                                                                          <  6:69837/148‑1 (MQ=255)
                           gccctcgcccGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGt                                                                                                                       >  6:304994/1‑148 (MQ=255)
                             cctcgcccGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCgg                                                                                                                            <  6:312900/141‑1 (MQ=255)
                             cctcgcccGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCgg                                                                                                                            >  5:312900/1‑141 (MQ=255)
                             cctcgcccGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCAc                                                                                                                    <  7:201710/149‑1 (MQ=255)
                             cctcgcccGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCAc                                                                                                                    >  4:29045/1‑149 (MQ=255)
                             cctcgcccGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCAc                                                                                                                    >  3:196852/1‑149 (MQ=255)
                                cgcccgccAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCgg                                                                                                                 >  8:8000/1‑149 (MQ=255)
                                 ggccgccAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGTGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGCCACCgg                                                                                                                 <  8:138912/146‑1 (MQ=255)
                                       cAGCCGGCATTTGGCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGt                                                                                                          <  3:32499/149‑1 (MQ=255)
                                       cAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGt                                                                                                          <  2:120501/149‑1 (MQ=255)
                                          ccGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACTATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCaa                                                                                                                                                 <  8:304977/107‑1 (MQ=255)
                                          ccGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACAATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCaa                                                                                                                                                 >  7:304977/1‑107 (MQ=255)
                                               aTTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTgg                                                                                                  >  1:155803/1‑149 (MQ=255)
                                                     cGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTgc                                                                                            <  8:291563/149‑1 (MQ=255)
                                                       gATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTg                                                                                             >  5:142561/1‑146 (MQ=255)
                                                       gATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTg                                                                                             <  6:142561/146‑1 (MQ=255)
                                                           gCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCTTAAGTGCCAGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATgc                                                                                      >  1:272798/1‑149 (MQ=255)
                                                                aGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCAt                                                                                 >  2:1993/1‑149 (MQ=255)
                                                                 gCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATg                                                                                >  7:217050/1‑149 (MQ=255)
                                                                   aaTTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTCCTGGCGCTGCGCATGCGCCATgc                                                                               >  4:285589/1‑148 (MQ=255)
                                                                    aTTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTgcgc                                                                                          >  6:36841/1‑136 (MQ=255)
                                                                    aTTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTgcgc                                                                                          <  5:36841/136‑1 (MQ=255)
                                                                     ttGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATgccag                                                                            >  4:148498/1‑149 (MQ=255)
                                                                       gcgcGAGCCTCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATgccagc                                                                           <  8:276544/145‑1 (MQ=255)
                                                                       gATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATgccagc                                                                           <  2:241433/148‑1 (MQ=255)
                                                                        aTCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCATGGCGCATGCGCAGCgccagcag                                                                         <  3:29045/149‑3 (MQ=255)
                                                                            aGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATgccag                                                                            >  1:234068/1‑142 (MQ=255)
                                                                            aGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATgccag                                                                            <  2:234068/142‑1 (MQ=255)
                                                                              cATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTtgct                                                                   <  3:128591/149‑1 (MQ=255)
                                                                                 gcGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATgc                                                                               <  8:234406/133‑1 (MQ=255)
                                                                                 ccGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATgc                                                                               >  7:234406/1‑134 (MQ=255)
                                                                                  cGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAg                                                               >  5:111662/1‑149 (MQ=255)
                                                                                    gggCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGc                                                             <  6:204689/149‑1 (MQ=255)
                                                                                      gCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCATGGTCGGCTTTTTTCAAGCGTTACAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCgg                                                           >  1:27012/1‑149 (MQ=255)
                                                                                      gCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTACAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTtgctgc                                                                 <  2:305027/143‑1 (MQ=255)
                                                                                      gCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTACAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTtgctgc                                                                 >  1:309705/1‑143 (MQ=255)
                                                                                      gCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTACAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTtgctgc                                                                 >  1:305027/1‑143 (MQ=255)
                                                                                      gCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTACAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTtgctgc                                                                 <  2:309705/143‑1 (MQ=255)
                                                                                       cAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCggg                                                          <  2:155803/149‑1 (MQ=255)
                                                                                        aaaCCTGGACCATGGAGAACATGGCCGGTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCgggg                                                         <  4:134028/149‑1 (MQ=255)
                                                                                          aCCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTTTCCCGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCa                                                       <  4:237248/149‑1 (MQ=255)
                                                                                              ggACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATAcc                                                   >  5:207907/1‑149 (MQ=255)
                                                                                                 ccATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGtt                                                <  1:197652/149‑1 (MQ=255)
                                                                                                   aTGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGa                                              <  7:120523/149‑1 (MQ=255)
                                                                                                       agaACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGAGCAg                                          >  8:59307/1‑149 (MQ=255)
                                                                                                           cATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGtt                                                <  6:183506/139‑1 (MQ=255)
                                                                                                           cATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGtt                                                >  5:183506/1‑139 (MQ=255)
                                                                                                            aTGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAgg                                     >  7:270319/1‑149 (MQ=255)
                                                                                                            aTGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAgg                                     <  7:179532/149‑1 (MQ=255)
                                                                                                             tGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCGAGTTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAGgc                                    <  8:217050/149‑1 (MQ=255)
                                                                                                                ccGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAGgcgcg                                 >  3:303609/1‑149 (MQ=255)
                                                                                                                  gCTTGCACCGGGTTGTCACGGTCGGCTTTTTTGAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAGGCGCGCa                               <  4:196852/149‑1 (MQ=255)
                                                                                                                    ttGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGGTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAGGCGCGCAgg                             >  8:212935/1‑149 (MQ=255)
                                                                                                                          cGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAGGCGCGCAGGTgggg                        >  6:251060/1‑148 (MQ=255)
                                                                                                                                   aCGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAGGCGCGCAGGTGGGGTATCAGTCg               <  1:1993/148‑1 (MQ=255)
                                                                                                                                   aCGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAGGCGCGCAGGTGGGGTATCAGTCg               <  5:304994/148‑1 (MQ=255)
                                                                                                                                    cGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAGGCGCGCAGGTGGGGTATCAGTCGGt             <  4:303609/149‑1 (MQ=255)
                                                                                                                                            ttttttCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATgcca                                                                             >  7:295673/1‑77 (MQ=255)
                                                                                                                                            ttttttCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATgcca                                                                             <  8:295673/77‑1 (MQ=255)
                                                                                                                                               tttCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGtt                                                <  3:160519/103‑1 (MQ=255)
                                                                                                                                               tttCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAGGCGCGCAGGTGGGGTATCAGTCGGTGGCCGCATTCa  >  1:23109/1‑149 (MQ=255)
                                                                                                                                               tttCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCAAGCCAGTTGCTGCAGGCGGGGCATACCGtt                                                >  4:160519/1‑103 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCAGGGCAATCGCTGGGCGGGCTGGTCGCCCTCGCCCGCCAGCCGGCATTTGCCGGATTGCTTGAGCAATTGATCGAGCATCCGGGGCAAACCTGGACCATGGAGAACATGGCCGCTTGCACCGGGTTGTCACGGTCGGCTTTTTTCAAGCGTTTCAGCGAGCTGGCCGGGCAGTCACCGGGACAAGTGCTGCTGGCGCTGCGCATGCGCCATGCCAGCCAGTTGCTGCAGGCGGGGCATACCGTTGAGCAGGTAGGCGCGCAGGTGGGGTATCAGTCGGTGGCCGCATTCA  >  NC_002947/4136344‑4136635

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: