Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,126,840:1 +G 6.6% coding (585/960 nt) PP_2744 ← ribose‑phosphate pyrophosphokinase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,126,8401.G 6.6% 248.7 / 10.6 61coding (585/960 nt)PP_2744ribose‑phosphate pyrophosphokinase family protein
Reads supporting (aligned to +/- strand):  ref base . (33/24);  new base G (2/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

GCGCGTAGCAGCGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCG  >  NC_002947/3126697‑3126982
                                                                                                                                                |                                                                                                                                              
gcgcGTAGCAGCGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCcagcag                                                                                                                                           <  8:24687/149‑1 (MQ=255)
 cgcgTAGCAGCGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAgg                                                                                                                                          >  6:387160/1‑149 (MQ=255)
  gcgTAGCAGCGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGt                                                                                                                                         <  5:387160/149‑1 (MQ=255)
  gcgTAGCAGCGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGt                                                                                                                                         >  8:33895/1‑149 (MQ=255)
        cagcGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGCCCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGcc                                                                                                                                    >  7:162656/1‑148 (MQ=255)
        cagcGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGcc                                                                                                                                    <  8:162656/148‑1 (MQ=255)
               tGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTg                                                                                                                            >  3:153651/1‑149 (MQ=255)
                 ccccGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGc                                                                                                                          >  4:208613/1‑149 (MQ=255)
                 ccccGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGc                                                                                                                          >  3:282777/1‑149 (MQ=255)
                   ccGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCg                                                                                                                         >  1:318110/1‑148 (MQ=255)
                       gCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGGCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCGAGCAGTTTCGCCAGGCCCTGGCgcacc                                                                                                                     >  5:199484/1‑145 (MQ=255)
                       gCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCACCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAggg                                                                                                                                <  6:333782/137‑1 (MQ=255)
                       gCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCACCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAggg                                                                                                                                >  5:333782/1‑137 (MQ=255)
                        cTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCACCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGTTTCGCCAGGCCCTGGCgcacct                                                                                                                    <  5:3392/148‑6 (MQ=255)
                           aTCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTc                                                                                                                >  7:321777/1‑149 (MQ=255)
                            tCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTgg                                                                                                                           >  6:127218/1‑137 (MQ=255)
                            tCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTgg                                                                                                                           <  5:127218/137‑1 (MQ=255)
                              aGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGGCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCGAGCAGTTTCGCCAGGCCCTGGCGCACCTGCtggc                                                                                                              <  6:199484/149‑4 (MQ=38)
                                      gAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAAc                                                                                                                     <  7:74079/133‑1 (MQ=255)
                                      gAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAAc                                                                                                                     >  8:74079/1‑133 (MQ=255)
                                          aCGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGt                                                                                                         >  6:272336/1‑141 (MQ=255)
                                          aCGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGt                                                                                                         <  5:272336/141‑1 (MQ=255)
                                               gACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCAc                                                                                               >  7:159340/1‑146 (MQ=255)
                                               gACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCAc                                                                                               <  8:159340/146‑1 (MQ=255)
                                               gACGGTGCTGCCTGCCACGTGACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCTGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGcgcg                                                                                                          <  6:86111/135‑1 (MQ=255)
                                               gACGGTGCTGCCTGCCACGTAACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGcgcg                                                                                                          >  5:86111/1‑135 (MQ=255)
                                                      ctgcctgcCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCgg                                                                                     <  1:313767/149‑1 (MQ=255)
                                                       tgcctgccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGt                                                                                    >  7:94878/1‑149 (MQ=255)
                                                         cctgccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGt                                                                                    <  8:94878/147‑1 (MQ=255)
                                                             ccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACg                                                                              >  3:301841/1‑149 (MQ=255)
                                                                   tACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGt                                                                         >  1:281907/1‑148 (MQ=255)
                                                                         ccAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGcc                                                                  <  4:153651/149‑1 (MQ=255)
                                                                          cAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCt                                                                 <  2:318110/149‑1 (MQ=255)
                                                                           aGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTg                                                                >  3:308116/1‑149 (MQ=255)
                                                                            ggTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGc                                                               >  1:262767/1‑149 (MQ=255)
                                                                            ggTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGc                                                               >  3:85099/1‑149 (MQ=255)
                                                                                        tCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGc                                                   >  3:200074/1‑149 (MQ=255)
                                                                                            gcctgcctGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTAGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCgggg                                               >  7:9043/1‑149 (MQ=255)
                                                                                             cctgcctgCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGc                                              <  8:321777/149‑1 (MQ=255)
                                                                                                 cctgcTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAg                                           <  7:33895/148‑1 (MQ=255)
                                                                                                  ctgctgGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGc                                         <  2:262767/149‑1 (MQ=255)
                                                                                                        gCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTg                                     <  4:282777/147‑1 (MQ=255)
                                                                                                          ggTGTTTCTCCATCAAGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGc                                  >  8:297477/1‑148 (MQ=255)
                                                                                                             gtTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGaa                               >  6:181839/1‑148 (MQ=255)
                                                                                                              tttCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGaa                               >  3:317229/1‑147 (MQ=255)
                                                                                                              tttCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGaa                               >  3:326783/1‑147 (MQ=255)
                                                                                                              tttCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCAAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGaa                               >  4:370357/1‑147 (MQ=255)
                                                                                                                    catcatGTCCACCGAGACCGGTCGCGGGAGGAGGTGCGCCAGGGCCTGGCGAAAGTGCTCGGCGCGTTTGACGCGACGGATATCCGGTTACACGACCGTGACCGCCTGCTCACCGAGTTGCTTGGCCAGGCGTTGGGCGAATAGCt                          <  4:157823/146‑1 (MQ=255)
                                                                                                                    catcatGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCt                          <  4:85099/146‑1 (MQ=255)
                                                                                                                    catcatGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTgg                        <  2:281907/148‑1 (MQ=255)
                                                                                                                        atGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGCGCCAGGCGTTGGGCGAATAGCTGGGTGCa                   >  4:48786/1‑149 (MQ=255)
                                                                                                                                cGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGt           >  3:134530/1‑149 (MQ=255)
                                                                                                                                 gagaCCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGtt          <  5:181839/149‑1 (MQ=255)
                                                                                                                                      cGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGc     >  4:337271/1‑149 (MQ=255)
                                                                                                                                       ggTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCt    >  4:39553/1‑149 (MQ=255)
                                                                                                                                         tCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCg  >  3:125955/1‑149 (MQ=255)
                                                                                                                                         tCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCg  <  3:208613/149‑1 (MQ=255)
                                                                                                                                         tCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCg  <  4:134530/149‑1 (MQ=255)
                                                                                                                                          cgcgCCGAGCAGTTTCGCCAGGCCCTGGCGCACCTGCTGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCa                                            >  2:74978/1‑106 (MQ=21)
                                                                                                                                          cgcgCCGAGCAGTTTCGCCAGGCCCTGGCGCACCTGCTGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCa                                            <  1:74978/106‑1 (MQ=21)
                                                                                                                                             gcCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTc         <  6:253314/138‑1 (MQ=255)
                                                                                                                                             gcCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTc         >  5:253314/1‑138 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GCGCGTAGCAGCGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCG  >  NC_002947/3126697‑3126982

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: