Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,820,817 C→G 5.7% S3458T (AGC→ACC)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,820,8170CG 5.7% 227.3 / 2.7 70S3458T (AGC→ACC) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (37/29);  new base G (2/2);  total (39/31)
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.73e-01

CTGTAGCGCGACAGCACCAGGGCCCAGGCCGCCTGGGCGAAAGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATCGCGGTAACGCGGTGGCAC  >  NC_002947/4820679‑4820965
                                                                                                                                          |                                                                                                                                                    
cTGTAGCGCGACAGCACCAGGGCCCAGGCCGCCTGGGCGAAAGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTg                                                                                                                                            <  6:174798/149‑1 (MQ=255)
     gcgcGACAGCACCAGGGCCCAGGCCGCCTGGGCGAAAGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCg                                                                                                                                                 >  2:285918/1‑139 (MQ=255)
     gcgcGACAGCACCAGGGCCCAGGCCGCCTGGGCGAAAGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCg                                                                                                                                                 <  1:285918/139‑1 (MQ=255)
                ccAGGGCCCAGGCCGCCTGGTCGATAGTGTTGACTGGGAGCTGTTGGGCCTGGGCCAGTTCGCGCAGAGCCACGCCATCATTCATCTCAAGGCGGGTATAACAATGGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAgctgg                                                                                                                            <  4:171585/149‑1 (MQ=255)
                 cAGGGCCCAGGCCGCCTGGGCGAAAGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAgcgggc                                                                                                                           >  8:85109/1‑149 (MQ=255)
                 cAGGGCCCAGGCCGACTGGGCGAAAGTGTTGAAGGTGAGCTGATGGGCCTGGGCAAGTTCGCGCAGACGCACGCCATCATTCATCTCAATGCGGGTATAAAAATCGCCGACGACCATGCCGCTCCCGGCGTGGTCATGGCCCAgcgg                                                                                                                             >  8:169914/1‑147 (MQ=255)
                     gCCCAGGCCGCCTGGGCGAAAGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTc                                                                                                                       >  6:248761/1‑149 (MQ=255)
                     gCCCAGGCCGCCTGGGCGAAAGTGTTGACGGTGAGCTGAAGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTc                                                                                                                       >  7:59635/1‑149 (MQ=255)
                                 tGGGCGAAAGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCAGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGGTCGGGATGGc                                                                                                           >  8:122977/1‑149 (MQ=255)
                                      gAAAGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTgg                                                                                                      >  1:190720/1‑149 (MQ=255)
                                      gAAAGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTgg                                                                                                      <  7:85109/149‑1 (MQ=255)
                                       aaaGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGGTGGCTGTGGc                                                                                                     <  5:206961/149‑1 (MQ=255)
                                       aaaGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTgg                                                                                                      >  5:267159/1‑148 (MQ=255)
                                        aaGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCa                                                                                                    >  4:187056/1‑149 (MQ=255)
                                           tgtTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGt                                                                                                 >  6:72966/1‑149 (MQ=255)
                                             ttGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTc                                                                                               >  2:130932/1‑149 (MQ=255)
                                             ttGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTc                                                                                               <  5:72966/149‑1 (MQ=255)
                                                aCGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGaa                                                                                            <  5:194323/149‑1 (MQ=255)
                                                 cGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAg                                                                                           <  6:23460/149‑1 (MQ=255)
                                                   gTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGcc                                                                                         >  1:185933/1‑149 (MQ=255)
                                                   gTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGcc                                                                                         <  6:4540/149‑1 (MQ=255)
                                                          gATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAgg                                                                                  >  6:19200/1‑149 (MQ=255)
                                                               gcctgggccAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGc                                                                             <  6:267159/149‑1 (MQ=255)
                                                                       cAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCAc                                                                                                   >  1:189083/1‑119 (MQ=255)
                                                                       cAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCAc                                                                                                   <  2:189083/119‑1 (MQ=255)
                                                                          ttCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAgg                                                                  >  3:19776/1‑149 (MQ=255)
                                                                          ttCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAgg                                                                  >  2:201920/1‑149 (MQ=255)
                                                                           tCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAgg                                                                  <  6:133721/148‑1 (MQ=255)
                                                                             gcgcAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCAc                                                               >  6:51339/1‑149 (MQ=255)
                                                                             gcgcAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCAc                                                               >  8:38298/1‑149 (MQ=255)
                                                                               gcAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGt                                                             >  3:78578/1‑149 (MQ=255)
                                                                                     gcacgcCATCATACATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGACAGGTTGGCCTGCCACCAGGCACGTGcgtcg                                                       >  5:279107/1‑149 (MQ=255)
                                                                                     acacgcCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATgg                                                                                                            >  7:284077/2‑96 (MQ=255)
                                                                                     acacgcCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATgg                                                                                                            <  8:284077/95‑1 (MQ=255)
                                                                                          ccATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAg                                                  <  5:248761/149‑1 (MQ=255)
                                                                                           catcatTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCCCCAGGCACTTGCGTCGTCCAgg                                                 >  3:13382/1‑149 (MQ=255)
                                                                                                   atCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGc                                          <  4:78578/148‑1 (MQ=255)
                                                                                                   atCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGc                                          <  5:51339/148‑1 (MQ=255)
                                                                                                      tcAAGGCGGGTATAACAATCGCCGACGACAATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTg                                      >  8:219/1‑149 (MQ=255)
                                                                                                        aaGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTg                                      <  3:187056/147‑1 (MQ=255)
                                                                                                                  taACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGtt                                                                                                <  7:214547/79‑1 (MQ=255)
                                                                                                                  taACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGtt                                                                                                >  8:214547/1‑79 (MQ=255)
                                                                                                                  taACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGGCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGAt                          <  6:279107/149‑1 (MQ=255)
                                                                                                                  taACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGAt                          <  1:130932/149‑1 (MQ=255)
                                                                                                                  taACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGAt                          >  3:22230/1‑149 (MQ=255)
                                                                                                                   aaCAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAgg                                                 >  2:262870/1‑125 (MQ=255)
                                                                                                                   aaCAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAgg                                                 <  1:262870/125‑1 (MQ=255)
                                                                                                                    aCAATCGCCGTCGTCCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATg                         <  5:325977/148‑1 (MQ=255)
                                                                                                                    aCAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATg                         <  2:190720/148‑1 (MQ=255)
                                                                                                                    aCAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATg                         <  2:185933/148‑1 (MQ=255)
                                                                                                                      aaTCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTaa                      >  6:186396/1‑149 (MQ=255)
                                                                                                                      aaTCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTaa                      >  7:290118/1‑149 (MQ=255)
                                                                                                                       aTCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAAcc                             >  5:205816/1‑141 (MQ=255)
                                                                                                                       aTCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAAcc                             <  6:205816/141‑1 (MQ=255)
                                                                                                                        tCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATc                    <  1:201920/149‑1 (MQ=255)
                                                                                                                         cGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATcg                   >  5:317077/1‑149 (MQ=255)
                                                                                                                           ccGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATcgcg                 <  5:186396/149‑1 (MQ=255)
                                                                                                                             gacgacCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATCGCgg                >  1:73499/1‑148 (MQ=255)
                                                                                                                              atgacCACGCCGGTAGCGGCATGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGt                        >  6:24340/3‑139 (MQ=255)
                                                                                                                              atgacCACGCCGGTAGCGGCATGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGt                        <  5:24340/137‑1 (MQ=255)
                                                                                                                                  ccACGCCGGTAGCGGCATGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGt                                                             >  2:110452/1‑98 (MQ=255)
                                                                                                                                  ccACGCCGGTAGCGGCATGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGt                                                             <  1:110452/98‑1 (MQ=255)
                                                                                                                                   cATGCCGCTACCGGCGTGGTCATGGCGCTGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATCGCGGTAAcgcg         >  1:176291/1‑149 (MQ=255)
                                                                                                                                    aTGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATCGCGGTAACGCgg        <  5:19200/149‑1 (MQ=255)
                                                                                                                                    aTGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGAGGTAATCGCGGTAACGCgg        >  3:235591/1‑149 (MQ=255)
                                                                                                                                     tGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTGGCTCGGGATGGCTGGGGCACGTTCGAAGGCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTACAGGCCCTGGCGTTGCAGCCAACCGATGTAATCGCGGTAACGCGGt       <  8:59635/149‑1 (MQ=255)
                                                                                                                                      gCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATCGCGGTAACGCGGTg      >  2:164982/1‑149 (MQ=255)
                                                                                                                                        cGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGATGAGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATCGCGGTAACGCTGtcg     >  7:145734/1‑146 (MQ=255)
                                                                                                                                        cGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCACGGGATCGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGAAAACGCGGTAAAGCGGTgg     >  2:290966/1‑148 (MQ=255)
                                                                                                                                         gCTACCGGCGTGGTTATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCTTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATCGCGGTAACGCGGTGGCa   >  5:181155/1‑149 (MQ=255)
                                                                                                                                          cTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTCGCGTTGCAGCCAACCGATGTAATCGCGGTAACGCGGTGGCAc  <  7:38298/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
CTGTAGCGCGACAGCACCAGGGCCCAGGCCGCCTGGGCGAAAGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATCGCGGTAACGCGGTGGCAC  >  NC_002947/4820679‑4820965

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑