Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,297,191 A→C 9.5% S58A (TCG→GCG)  TpbB ← diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,297,1910AC 9.5% 245.4 / 12.0 81S58A (TCG→GCG) TpbBdiguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (33/40);  new base C (4/4);  total (37/44)
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

TCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCAC  >  NC_002947/5297043‑5297338
                                                                                                                                                    |                                                                                                                                                   
tCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGa                                                                                                                                                     >  8:147500/1‑149 (MQ=255)
        ccAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCg                                                                                                                                              >  4:228216/1‑148 (MQ=255)
        ccAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGa                                                                                                                                             <  2:29586/149‑1 (MQ=255)
            tGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGCCCGCTCGATCaa                                                                                                                                         <  2:215983/149‑1 (MQ=255)
              cgcCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGCCCGCTCGATCAAc                                                                                                                                        <  5:727/148‑1 (MQ=255)
              cgcCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGCCCGCTCGATCAAc                                                                                                                                        <  2:239141/148‑1 (MQ=255)
              cgcCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAAc                                                                                                                                        <  5:69790/148‑1 (MQ=255)
                ccAGGTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGc                                                                                                                                     <  7:147500/149‑1 (MQ=255)
                 cAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGACGCTTCTACGGTGTAGTTGATCGAGCGGTCGATCAACTg                                                                                                                                      >  3:223459/1‑147 (MQ=255)
                       cTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTa                                                                                                                               <  4:177914/148‑1 (MQ=255)
                       cTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTa                                                                                                                               <  3:138465/148‑1 (MQ=255)
                       cTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTAt                                                                                                                              <  4:71866/149‑1 (MQ=255)
                        tCGCCCTCGTTGTTGAAGCCCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCTAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATg                                                                                                                             >  4:29397/1‑149 (MQ=255)
                        tCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATg                                                                                                                             >  2:20329/1‑149 (MQ=255)
                            ccTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTg                                                                                                                         >  5:281422/1‑149 (MQ=255)
                            ccTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTg                                                                                                                         <  4:262984/149‑1 (MQ=255)
                            ccTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTg                                                                                                                         <  4:267338/149‑1 (MQ=255)
                            ccTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTg                                                                                                                         >  2:215168/1‑149 (MQ=255)
                                 ttgttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACTGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATa                                                                                                                    <  4:103014/149‑1 (MQ=255)
                                    ttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGac                                                                                                                 <  5:279154/149‑1 (MQ=255)
                                      gAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                                                                                        >  5:290987/1‑140 (MQ=255)
                                      gAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                                                                                        <  6:290987/140‑1 (MQ=255)
                                      gAAGACCTTGGCCTCGGCCACTTCTTCGGTAATGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACAACACAGCCGCTTCTACGGTGTAGTTGATCGAGAGGGCGATCAACTGCATGTTATGGTTGGCAAAGacac                                                                                                               >  6:319015/1‑149 (MQ=255)
                                       aaGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACg                                                                                                              >  4:288838/1‑149 (MQ=255)
                                        aGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGc                                                                                                             <  4:223459/149‑1 (MQ=255)
                                         gACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGTTTATGGTTGGCATAGACACg                                                                                                              >  8:257029/1‑147 (MQ=255)
                                            cTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAggg                                                                                                         >  4:282013/1‑149 (MQ=255)
                                                        cACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCca                                                                                             <  2:107988/149‑1 (MQ=255)
                                                           ttcttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcag                                                                                            >  4:306428/1‑147 (MQ=255)
                                                           ttcttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcag                                                                                            >  6:142347/1‑147 (MQ=255)
                                                             cttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcagcag                                                                                         <  3:29397/148‑1 (MQ=255)
                                                                cGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAgg                                                                                        <  7:266279/146‑1 (MQ=255)
                                                                 ggTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGt                                                                                      >  8:179408/1‑147 (MQ=255)
                                                                 ggTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGa                                                                                    <  3:306428/149‑1 (MQ=255)
                                                                  gTACTGGCAATCAACTGCAGCGATTAGTTGGCGGCTGCGCTGTAATCGAACACTACAGCCGCTTCTACTGGTTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTTTCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTTAg                                                                                   <  5:319015/149‑1 (MQ=255)
                                                                   tACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                  <  1:215168/149‑1 (MQ=255)
                                                                   tACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                  <  5:142347/149‑1 (MQ=255)
                                                                        gCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGgag                                                                               <  6:242566/147‑1 (MQ=255)
                                                                            tCACCGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAccc                                                                         <  1:140552/149‑1 (MQ=255)
                                                                                 gCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAccc                                                                         >  7:251772/1‑144 (MQ=255)
                                                                                 gCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAccc                                                                         <  8:251772/144‑1 (MQ=255)
                                                                                 gCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAg                                                                    <  5:294984/149‑1 (MQ=255)
                                                                                  ccAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAgg                                                                   <  8:259890/149‑1 (MQ=255)
                                                                                   cAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGc                                                                  <  2:206275/149‑1 (MQ=255)
                                                                                   cAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGc                                                                  <  3:228216/149‑1 (MQ=255)
                                                                                        aTTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACg                                                             >  1:285075/1‑149 (MQ=255)
                                                                                         ttCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAg                                                                    >  4:302051/1‑141 (MQ=255)
                                                                                         ttCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAg                                                                    <  3:302051/141‑1 (MQ=255)
                                                                                            gTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGgcaa                                                         <  2:285075/149‑1 (MQ=255)
                                                                                             ttGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                  >  7:114882/1‑123 (MQ=255)
                                                                                             ttGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                  <  8:114882/123‑1 (MQ=255)
                                                                                             ttGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                  <  8:110601/123‑1 (MQ=255)
                                                                                             ttGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                                                                                  >  7:110601/1‑123 (MQ=255)
                                                                                                      gcgcTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTg                                                                                                                         <  1:143514/75‑1 (MQ=255)
                                                                                                      gcgcTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTg                                                                                                                         >  2:143514/1‑75 (MQ=255)
                                                                                                       cgcTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGc                                                                                                        >  8:219128/1‑91 (MQ=255)
                                                                                                       cgcTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGc                                                                                                        <  7:219128/91‑1 (MQ=255)
                                                                                                       cgcTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAgg                                                                                        <  7:193981/107‑1 (MQ=255)
                                                                                                       cgcTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAgg                                                                                        >  8:193981/1‑107 (MQ=255)
                                                                                                              aTCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGt                                       >  6:120429/1‑149 (MQ=255)
                                                                                                              aTCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGt                                       >  7:276923/1‑149 (MQ=255)
                                                                                                                 gAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAgg                                                                   >  8:55953/1‑118 (MQ=255)
                                                                                                                  aaCACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGc                                    >  7:106752/1‑148 (MQ=255)
                                                                                                                      ccacAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACgg                               >  2:285082/1‑149 (MQ=255)
                                                                                                                       cacaGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGc                              >  6:65031/1‑149 (MQ=255)
                                                                                                                       cacaGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGc                              >  1:168644/1‑149 (MQ=255)
                                                                                                                       cacaGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGCTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGc                              >  7:97732/1‑149 (MQ=255)
                                                                                                                        acaGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGcc                             >  4:271882/1‑149 (MQ=255)
                                                                                                                          aGCCGCTTCTACGTTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCa                           <  6:281422/149‑1 (MQ=255)
                                                                                                                                  cTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGa                    <  1:285082/148‑1 (MQ=255)
                                                                                                                                       ctgcaGTTGATCGCCCGCTCGATCAACTACAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAg              >  2:99547/5‑149 (MQ=255)
                                                                                                                                          tAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGt           >  3:146550/1‑149 (MQ=255)
                                                                                                                                          tAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGt           <  4:146550/149‑1 (MQ=255)
                                                                                                                                          tAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGt           <  5:120429/149‑1 (MQ=255)
                                                                                                                                           aGTTGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTc          >  4:128043/1‑149 (MQ=255)
                                                                                                                                            gatGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCg         >  8:135551/3‑149 (MQ=255)
                                                                                                                                             ttGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGc                                                  <  3:13475/107‑1 (MQ=255)
                                                                                                                                             ttGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGc                                                  >  4:13475/1‑107 (MQ=255)
                                                                                                                                             ttGATCGCCCGCTCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCgg        >  2:266299/1‑149 (MQ=255)
                                                                                                                                               gATCGAGCGGGCGATCAACTGCAGTTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCggg       <  7:257029/148‑1 (MQ=255)
                                                                                                                                                aTCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGgc      <  2:168644/148‑1 (MQ=255)
                                                                                                                                                   gAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCAc  <  8:97732/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAGCAAGGCGACGCTGAGGTGGCCACGGCCCAGCACCGAGCGCAGGGTCGGGCGCAC  >  NC_002947/5297043‑5297338

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