Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,250,567 G→A 5.5% P40S (CCG→TCG)  cidB ← holin CidA anti‑holin regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,250,5670GA 5.5% 220.7 / 2.0 72P40S (CCG→TCG) cidBholin CidA anti‑holin regulator
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base A (2/2);  total (34/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 = 7.59e-01

GCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCACGACGCAGCCTCCAGCGGCACACCATCT  >  NC_002947/5250423‑5250711
                                                                                                                                                |                                                                                                                                                
gcgcCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCAc                                                                                                                                               >  5:365221/1‑148 (MQ=255)
 cgcCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCa                                                                                                                                             <  3:181738/149‑1 (MQ=255)
                  gctgcCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGACACCAGGATCAACGGCGAcag                                                                                                                             >  2:25418/1‑148 (MQ=255)
                    tgcCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGAca                                                                                                                              <  8:345097/145‑1 (MQ=255)
                    tgcCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGAca                                                                                                                              >  7:345097/1‑145 (MQ=255)
                        aGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGTGGCACGGCCACAGCCAGCAGCAACACCGACa                                                                                                                                                >  8:173050/1‑123 (MQ=255)
                        aGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGACa                                                                                                                                                <  7:173050/123‑1 (MQ=255)
                         gATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAgc                                                                                                                     >  4:96928/1‑149 (MQ=255)
                           tCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCA‑CAGGATCAACGGCGACAGCAGGTAGCGCt                                                                                                                  <  4:80120/149‑1 (MQ=255)
                             ggCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGTCATAGGACGTATGCAGCGGCATGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCt                                                                                                                  <  5:272516/148‑1 (MQ=255)
                               cACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGc                                                                                                               <  6:114602/149‑1 (MQ=255)
                                  cgcgAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTAc                                                                                                            >  5:260849/1‑149 (MQ=255)
                                   gcgAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCT‑GCTACGc                                                                                                          >  5:220279/1‑149 (MQ=255)
                                     gAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGcc                                                                                                         >  8:323313/1‑149 (MQ=255)
                                         gCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATa                                                                                                     >  2:240805/1‑149 (MQ=255)
                                           gACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATAc                                                                                                    >  7:81675/1‑148 (MQ=255)
                                            aCAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATAcaa                                                                                                  <  2:381240/149‑1 (MQ=255)
                                             cAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATAcaac                                                                                                 >  6:197890/1‑149 (MQ=255)
                                              aGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTgg                                                                                                                >  4:351547/1‑133 (MQ=255)
                                              aGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTgg                                                                                                                <  3:351547/133‑1 (MQ=255)
                                               gTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATAcaacca                                                                                               <  8:81675/149‑1 (MQ=255)
                                                      ggCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTg                                                                                        >  5:253742/1‑149 (MQ=255)
                                                       gCCGAGCACGCCCATCAACCTGTGGGTGTTGCGCGCATACTCTGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGTCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGc                                                                                       <  1:157117/149‑1 (MQ=255)
                                                       gCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGc                                                                                       <  4:211045/149‑1 (MQ=255)
                                                          gAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCa                                                                                    <  6:365221/149‑1 (MQ=255)
                                                            gCACGCCCATCCACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAgg                                                                                  <  6:386966/149‑1 (MQ=255)
                                                                gCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAg                                                                              >  1:183123/1‑149 (MQ=255)
                                                                  acATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGc                                                                             >  2:213342/2‑148 (MQ=255)
                                                                   cATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGc                                                                             >  3:378494/1‑147 (MQ=255)
                                                                   cATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCa                                                                           >  4:218389/1‑149 (MQ=255)
                                                                     tCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGc                                                                         <  7:323313/149‑1 (MQ=255)
                                                                       aaCCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGCCAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGc                                                                         >  5:99002/1‑147 (MQ=255)
                                                                       aaCCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGc                                                                         >  8:336495/1‑147 (MQ=255)
                                                                        aCCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAg                                                                      <  4:68247/149‑1 (MQ=255)
                                                                         ccAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGc                                                                         >  1:183886/1‑145 (MQ=255)
                                                                         ccAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGc                                                                         <  2:183886/145‑1 (MQ=255)
                                                                          cAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCATGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGt                                                                    >  7:209497/1‑149 (MQ=255)
                                                                           aGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTc                                                                   <  7:336495/149‑1 (MQ=255)
                                                                            gTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTc                                                                   >  5:24288/1‑148 (MQ=255)
                                                                            gTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTc                                                                   >  5:29803/1‑148 (MQ=255)
                                                                            gTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCa                                                                  <  3:96928/149‑1 (MQ=255)
                                                                               ggTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGccc                                                                                     >  8:97981/1‑127 (MQ=255)
                                                                               ggTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGccc                                                                                     <  7:97981/127‑1 (MQ=255)
                                                                                gtgtTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGTCACAGCCAGCAGCAACACAGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACCGCTGCACAGCTAAGCCAGCAAGGTCAAGGc                                                              >  8:300062/1‑149 (MQ=255)
                                                                                   ttGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCaa                                                           >  4:351934/1‑149 (MQ=255)
                                                                                   ttGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCaa                                                           <  3:218389/149‑1 (MQ=255)
                                                                                     gcgcgcATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAAc                                                          <  1:213342/148‑1 (MQ=255)
                                                                                           aTACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACa                                                   <  2:162635/149‑1 (MQ=255)
                                                                                            tACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAg                                                  >  3:145620/1‑149 (MQ=255)
                                                                                            tACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCT‑GCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGc                                                 <  6:220279/149‑1 (MQ=255)
                                                                                             aCTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGCCAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGCACAGc                                                 >  2:29641/1‑149 (MQ=255)
                                                                                               tCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAg                                                                      >  6:340467/1‑126 (MQ=255)
                                                                                               tCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAg                                                                      <  5:340467/126‑1 (MQ=255)
                                                                                                    aTAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATgg                                          <  6:253742/149‑1 (MQ=255)
                                                                                                    aTAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATgg                                          <  5:197890/149‑1 (MQ=255)
                                                                                                     tAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATggg                                         <  6:260849/149‑1 (MQ=255)
                                                                                                          cGTATGCAGCGGCACGGTCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCaa                                    <  7:300062/149‑1 (MQ=255)
                                                                                                          cGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCaa                                    <  4:378494/149‑1 (MQ=255)
                                                                                                           gTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGACACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAACAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAg                                   <  1:25418/149‑1 (MQ=255)
                                                                                                                cAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCa                              <  6:24288/149‑1 (MQ=255)
                                                                                                                cAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCa                              <  6:29803/149‑1 (MQ=255)
                                                                                                                 aGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCac                             <  3:351934/149‑1 (MQ=255)
                                                                                                                 aGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCac                             >  6:142014/1‑149 (MQ=255)
                                                                                                                  gCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCac                             >  1:130173/1‑148 (MQ=255)
                                                                                                                      cACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCACGACGc                        <  1:157783/149‑1 (MQ=255)
                                                                                                                       aCGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCacgacg                         >  3:218143/1‑147 (MQ=255)
                                                                                                                            cacaGCCAGCGGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAg                                   >  1:83162/1‑132 (MQ=255)
                                                                                                                            cacaGCCAGCGGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAg                                   <  2:83162/132‑1 (MQ=255)
                                                                                                                            cacaGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCACGACGCAGCCTc                  <  1:240805/149‑1 (MQ=255)
                                                                                                                               agccagcAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCACGACGCAGCCTCCAg               <  8:209497/149‑1 (MQ=255)
                                                                                                                                            acCGGCACCAGGATCAACGGCGAGAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCACGACGCAGCCTCCAGCGGCACACCAtct  <  2:134402/149‑1 (MQ=255)
                                                                                                                                            acCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCACGACGCAGCCTCCAGCGGCACACCAtct  >  2:342694/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGGTAGCGCTGGCTACGCCGATACAACCAACGGCTGCCCAGGTAAGCCAGCAAGGTCAAGGCCAACCAGAACAGCGGCATGGGTTCAAGGCTCACGACGCAGCCTCCAGCGGCACACCATCT  >  NC_002947/5250423‑5250711

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: ‑