Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,772,222 G→T 7.2% R2349S (CGC→AGC)  pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,772,2220GT 7.2% 217.9 / 4.9 69R2349S (CGC→AGC) pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (31/33);  new base T (2/3);  total (33/36)
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.64e-01

TCTGCGCGACACTGAAACGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATATGGTGGTAGGTCCAGATCAGTTGATGGCGCTGATC  >  NC_002947/4772075‑4772369
                                                                                                                                                   |                                                                                                                                                    
tCTGCGCGACACTGAAACGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGgcg                                                                                                                                                      >  5:379500/1‑148 (MQ=255)
    cgcgACACTGAAACGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGcagc                                                                                                                                                  <  6:177872/148‑1 (MQ=255)
    cgcgACACTGAAACGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGcagc                                                                                                                                                  <  2:248718/148‑1 (MQ=255)
       gACACTGAAACGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGcagccagc                                                                                                                                              <  7:174062/149‑1 (MQ=255)
        acacTGAAACGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCTgcgccagca                                                                                                                                             <  5:95262/149‑2 (MQ=255)
         cacTGAAACGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCTgcgccagcag                                                                                                                                            >  6:169365/1‑147 (MQ=255)
         cacTGAAACGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCTgcgccagcag                                                                                                                                            <  5:169365/149‑3 (MQ=255)
          acTGAAACGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGa                                                                                                                                           <  8:124296/149‑1 (MQ=255)
            tGAAACGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATg                                                                                                                                         >  6:17954/1‑149 (MQ=255)
              aaaCGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATg                                                                                                                                         <  4:127739/147‑1 (MQ=255)
               aaCGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAg                                                                                                                                      <  1:12070/149‑1 (MQ=255)
                     tGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCTgcgccagcagg                                                                                                                                           >  8:281383/1‑135 (MQ=255)
                     tGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCTgcgccagcagg                                                                                                                                           <  7:281383/138‑4 (MQ=255)
                     tGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTa                                                                                                                                       >  6:118737/1‑142 (MQ=255)
                     tGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTa                                                                                                                                       <  5:118737/142‑1 (MQ=255)
                     tGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCg‑g                                                                                                                               >  7:347459/1‑149 (MQ=255)
                     tGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCg‑g                                                                                                                               >  7:223254/1‑149 (MQ=255)
                        agagGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑Gt                                                                                                                              >  8:263147/1‑147 (MQ=255)
                         gagGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTaa                                                                                                                            <  3:120795/148‑1 (MQ=255)
                          agGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAgg                                                                                                                          >  1:152709/1‑149 (MQ=255)
                            ggCCTGGTGGCCGCTGCCGCGCTCTCCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGca                                                                                                                        >  7:233322/1‑149 (MQ=255)
                               cTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCg‑g                                                                                                                               >  7:193096/1‑139 (MQ=255)
                               cTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCg‑g                                                                                                                               <  8:193096/139‑1 (MQ=255)
                               cTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGcacca                                                                                                                     >  3:102945/1‑149 (MQ=255)
                                    ggCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCGCCTGCAGGCTGGCCAGCGGCTGCTTCCAGAACGTTTCGCGGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTg                                                                                                                <  8:211464/149‑1 (MQ=255)
                                       cGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTc                                                                                                              >  6:306538/1‑148 (MQ=255)
                                            ccGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGc                                                                                                          <  5:99592/147‑1 (MQ=255)
                                            ccGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGt                                                                                                        >  2:127613/1‑149 (MQ=255)
                                                      cGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGaa                                                                                              >  7:319157/1‑149 (MQ=255)
                                                       gCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGaa                                                                                              >  7:247641/1‑148 (MQ=255)
                                                          ggCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATaa                                                                                           >  8:97246/1‑148 (MQ=255)
                                                           gCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACa                                                                                         <  7:263147/149‑1 (MQ=255)
                                                               aGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGcc                                                                                     <  7:11862/149‑1 (MQ=255)
                                                                 gCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAg                                                                                   <  4:101232/149‑1 (MQ=255)
                                                                   aTCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACTTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGATCGTCTGCCCCGAATAACAGGCCAGGa                                                                                 <  8:223254/149‑1 (MQ=255)
                                                                   aTCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAgg                                                                                  >  6:365407/1‑148 (MQ=255)
                                                                        ccAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCg                                                                            >  5:374949/1‑149 (MQ=255)
                                                                          aggtaggtCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCGCTTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGc                                                                           <  1:97021/149‑1 (MQ=255)
                                                                           ggtaggtCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGccc                                                                         <  8:347459/149‑1 (MQ=255)
                                                                           ggtaggtCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGccc                                                                         <  8:247641/149‑1 (MQ=255)
                                                                               ggtCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATc                                                                      <  2:152709/148‑1 (MQ=255)
                                                                                 tCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAAc                                                                   >  6:322947/1‑149 (MQ=255)
                                                                                   ggCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACt                                                                  <  1:127613/148‑1 (MQ=255)
                                                                                       ccTGCAGGCTGGCAAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGATGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACAACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATAAAATGGGAc                                                             >  8:175756/1‑149 (MQ=255)
                                                                                        cTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGAc                                                             <  7:97246/148‑1 (MQ=255)
                                                                                         tGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACg                                                            >  8:172387/1‑148 (MQ=255)
                                                                                           cAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTa                                                          >  7:255954/1‑148 (MQ=255)
                                                                                           cAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTAc                                                         <  7:172387/149‑1 (MQ=255)
                                                                                                 ggCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGc                                                   <  8:255954/149‑1 (MQ=255)
                                                                                                     aGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACCTCGCCCATCAACTGCGACGTACTCCAGCCATc                                               <  4:102945/149‑1 (MQ=255)
                                                                                                       ctgctgctTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGa                                             <  3:332103/149‑1 (MQ=255)
                                                                                                             ctTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGa                                       >  2:37575/1‑149 (MQ=255)
                                                                                                              ttCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGa                                       >  1:213999/1‑148 (MQ=255)
                                                                                                                ccAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGata                                     <  5:322947/148‑1 (MQ=255)
                                                                                                                 cAGCACCTTTCGCTGGGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGGGTCTGCCACGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATAtg                                   <  3:276027/149‑1 (MQ=255)
                                                                                                                  aGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGCCTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATAtgg                                  >  5:289465/1‑149 (MQ=255)
                                                                                                                  aGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTc                                                       >  3:86585/1‑128 (MQ=255)
                                                                                                                  aGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTc                                                       <  4:86585/128‑1 (MQ=255)
                                                                                                                  aGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATAtgg                                  >  3:58581/1‑149 (MQ=255)
                                                                                                                    aaCGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATAtggtg                                >  5:272774/1‑149 (MQ=255)
                                                                                                                    aaCGTTTCGCGGTGGCGCGGAGCCTGCGGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCCAATAACAGGCCAGGACTTCGCCCATCAACTGCGCCGTACTCCAGCCATCGCGCAAGTTAtggtg                                <  6:32572/149‑1 (MQ=255)
                                                                                                                                  gcgcGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGATAGATATGGTGGTTGGTCCAGATCa                  >  8:257616/1‑149 (MQ=255)
                                                                                                                                     cgGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATATGGTGGTAGGTCCAGATCAGtt               <  8:233322/149‑1 (MQ=255)
                                                                                                                                        aTCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATATGGTGGTAGGTCCAGATCAGTTGAt            >  8:286414/1‑149 (MQ=255)
                                                                                                                                         tCCTGCTGGCGCAGCCAGCAGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATATGGTGGTAGGTCCAGATCAGTTGATg           <  7:175756/149‑1 (MQ=255)
                                                                                                                                            tgctgGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATAtggtg                                <  8:117363/125‑1 (MQ=255)
                                                                                                                                            tgctgGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATAtggtg                                >  7:117363/1‑125 (MQ=255)
                                                                                                                                            tgctgGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATATGGTGGTAGGTCCAGATCAGTTGATGgcg        >  5:254773/1‑149 (MQ=255)
                                                                                                                                             gctgGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATATGGTGGTAGGTCCAGATCAGTTGATGgcgc       >  7:61084/1‑149 (MQ=255)
                                                                                                                                              cttgagcAGCCAGCCGATGTAGTTGCG‑GTAATGCACCACTCTGGCGAGCGTCTGCCCCGAATAACAGGCCAGGTCTTCGCCCATCAACTGCGACGTAATCCAGCAATCGAGCAAGATATGGTGGTAGGTCCAGATCAGTTGAGGgcgc       <  7:72641/143‑1 (MQ=255)
                                                                                                                                              ctgGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATATGGTGGTAGGTCCAGATCAGTTGATGgcgc       <  8:319157/148‑1 (MQ=255)
                                                                                                                                                 gcgcAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATATGGTGGTAGGTCCAGATCAGTTGATGGCGCTGAt   >  7:110595/1‑149 (MQ=255)
                                                                                                                                                  cgcAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATATGGTGGTAGGTCCAGATCAGTTGATGGCGCTGATc  <  6:374949/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCTGCGCGACACTGAAACGGCTGTAGAGGGCCTGGTGGCCGCTGCCGCGCTCACCGCTGGCGAAGGCATCGGCCAGGTAGGTCGGCTCCTGCAGGCTGGCCAGCTGCTGCTTCCAGAACGTTTCGCTGTGGCGCGGATCCTGCTGGCGCAGCCAGCCGATGTAGTTGCG‑GTAAGGCACCACGCTGTCGAGCGTCTGCCCCGAATAACAGGCCAGGACTTCGCCCATCAACTGCGACGTACTCCAGCCATCGAGCAAGATATGGTGGTAGGTCCAGATCAGTTGATGGCGCTGATC  >  NC_002947/4772075‑4772369

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