Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,234,432 A→G 6.6% R58R (CGT→CGC)  tyrB ← aromatic amino acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,234,4320AG 6.6% 175.2 / 2.0 60R58R (CGT→CGC) tyrBaromatic amino acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (25/31);  new base G (2/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

CGCCGGTACCACCGACGGCCTGCACGGTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGAATGCTTCGTTGAGGCCCAGAATAGGGTCGCGGGGT  >  NC_002947/2234284‑2234576
                                                                                                                                                    |                                                                                                                                                
cgcCGGTACCACCGACGGCCTGCACGGTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCa                                                                                                                                                  >  2:179805/1‑149 (MQ=255)
      tACCACCGACGGCCTGCACGGTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTgg                                                                                                                                            <  3:203434/149‑1 (MQ=255)
        ccaccGACGGCCTGCACGGTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCGCGGTGGGt                                                                                                                                          >  7:296453/1‑149 (MQ=255)
        ccaccGACGGCCTGCACGGTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCGCGCTGGGt                                                                                                                                          >  7:300649/1‑149 (MQ=255)
         caccGACGGCCTGCACGGTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGtt                                                                                                                                         >  6:230429/1‑149 (MQ=255)
               cGGCCTGCACGGTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCt                                                                                                                                   <  5:83921/149‑1 (MQ=255)
                 gCCTGCACGGTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTc                                                                                                                                 <  1:179805/149‑1 (MQ=255)
                   cTGCACGGTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGa                                                                                                                               >  7:231590/1‑149 (MQ=255)
                   cTGCACGGTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGa                                                                                                                               >  2:155775/1‑149 (MQ=255)
                        cGGTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACg                                                                                                                          <  5:230429/149‑1 (MQ=255)
                         ggTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACgg                                                                                                                         >  4:293826/1‑149 (MQ=255)
                              ccACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGgcgcgc                                                                                                                    <  1:155775/149‑1 (MQ=255)
                               cACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGgcgcg                                                                                                                     >  5:289045/1‑147 (MQ=255)
                               cACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGgcgcg                                                                                                                     >  8:197758/1‑147 (MQ=255)
                                aCGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCgcag                                                                                                                  <  6:289045/149‑1 (MQ=255)
                                        ggcggcCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACg                                                                                                                          <  1:171002/133‑1 (MQ=255)
                                        ggcggcCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACg                                                                                                                          >  2:171002/1‑133 (MQ=255)
                                            gcCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCg                                                                                                      <  6:281876/149‑1 (MQ=255)
                                                  agTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTc                                                                                                <  6:251090/149‑1 (MQ=255)
                                                   gTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTcc                                                                                               <  8:231590/149‑1 (MQ=255)
                                                    tGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCt                                                                                              <  1:152722/149‑1 (MQ=255)
                                                      gCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCa                                                                                            >  4:157780/1‑149 (MQ=255)
                                                             ggCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAAt                                                                                     <  7:197758/149‑1 (MQ=255)
                                                                gcCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCagcagc                                                                                                              <  4:93385/121‑1 (MQ=255)
                                                                gcCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCagcagc                                                                                                              >  3:93385/1‑121 (MQ=255)
                                                                gcCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTACGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATaca                                                                                  >  6:204773/1‑149 (MQ=255)
                                                                    aaCAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGcc                                                                              >  3:21292/1‑149 (MQ=255)
                                                                    aaCAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGcc                                                                              >  7:56127/1‑149 (MQ=255)
                                                                       agcagcTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTAc                                                                           >  7:206735/1‑149 (MQ=255)
                                                                       agcagcTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTAc                                                                           <  3:157780/149‑1 (MQ=255)
                                                                       agcagcTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTAc                                                                           <  2:300462/149‑1 (MQ=255)
                                                                         cagcTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGc                                                                         >  7:253707/1‑149 (MQ=255)
                                                                         cagcTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTAGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGc                                                                         >  1:184055/1‑149 (MQ=255)
                                                                           gcTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGccc                                                                       <  8:56127/149‑1 (MQ=255)
                                                                             ttCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAg                                                                     >  1:193967/1‑149 (MQ=255)
                                                                               cTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGGCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAg                                                                     <  5:204773/147‑1 (MQ=255)
                                                                               cTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAg                                                                     <  2:307472/147‑1 (MQ=255)
                                                                                           ggcGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCGCGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCg                                                        <  8:296453/146‑1 (MQ=255)
                                                                                           gTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCGCGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCg                                                        <  8:300649/148‑1 (MQ=255)
                                                                                           gTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCg                                                        <  4:21292/148‑1 (MQ=255)
                                                                                              gTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAAt                                                                                     <  1:275421/116‑1 (MQ=255)
                                                                                              gTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAAt                                                                                     >  2:275421/1‑116 (MQ=255)
                                                                                                  gTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGtt                                                                  >  4:247007/1‑131 (MQ=255)
                                                                                                  gTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGtt                                                                  <  3:247007/131‑1 (MQ=255)
                                                                                                  gTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTgg                                                >  2:179752/1‑149 (MQ=255)
                                                                                                   tGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTc                                                         >  5:245216/1‑139 (MQ=255)
                                                                                                   tGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTc                                                         <  6:245216/139‑1 (MQ=255)
                                                                                                        aTGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGc                                          <  3:293826/149‑1 (MQ=255)
                                                                                                            cGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGtt                                       >  5:142882/1‑148 (MQ=255)
                                                                                                             gTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTg                                      >  7:58047/1‑148 (MQ=255)
                                                                                                              tCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGaa                                    <  8:206735/149‑1 (MQ=255)
                                                                                                              tCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGaa                                    <  2:193967/149‑1 (MQ=255)
                                                                                                                         gTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGAATGCTTCGtt                           >  1:213021/1‑147 (MQ=255)
                                                                                                                          tAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGAATGCTTCGTTGAg                        <  8:253707/149‑1 (MQ=255)
                                                                                                                                   ccaGCGTGCCGCCCAGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGAATGCTTCGTTGAGGCCCAGaa                >  3:21700/3‑148 (MQ=255)
                                                                                                                                    aaGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGAATGCTTCGTTGAGGCCCAGAATa              <  6:142882/149‑1 (MQ=255)
                                                                                                                                     aGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGAATGCTTCGTTGAGGCCCAGAATa              <  8:58047/148‑1 (MQ=255)
                                                                                                                                         tGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGAATGCTTCGTTGAGGCCCAGAATAGGGTc         <  2:184055/149‑1 (MQ=255)
                                                                                                                                                cagcACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGAATGCTTCGTTGAGGCCCAGAATAGGGTCGCGGGGt  <  4:21700/147‑1 (MQ=255)
                                                                                                                                                  gcACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGAATGCTTCGTTGAGGCCCAGAATAGGGTCGCGGGGt  <  7:117878/147‑1 (MQ=255)
                                                                                                                                                  gcACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGAATGCTTCGTTGAGGCCCAGAATAGGGTCGCGGGGt  >  8:117878/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGCCGGTACCACCGACGGCCTGCACGGTGACCACGCGGCCGGCGGCCAGCAGTGGCGACTCGGCGCCGAACAGCAGCTTCTGCACAGCCTGGTCGTAGGTGGCGATGCCGTCGATCGGCAGGTAGCCGCGCGAAGCGTGCTGGGCGGCACGCTGGGTTTCGGCTTCGATCACGGCGCGCAGCAGCGGAATGCGGCCTTCCTCATTGCAATACACGCCTACGCCCAGGTTGACCTTGTCGGTACGTGGGTCGGCGTTGAATGCTTCGTTGAGGCCCAGAATAGGGTCGCGGGGT  >  NC_002947/2234284‑2234576

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