Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,820,812 T→C 6.1% M3460V (ATG→GTG)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,820,8120TC 6.1% 191.8 / 2.1 65M3460V (ATG→GTG) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (33/28);  new base C (2/2);  total (35/30)
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.17e-01

CTGTAGCGCGACAGCACCAGGGCCCAGGCCGCCTGGGCGAAAGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATCGCGGTAACGCGGT  >  NC_002947/4820679‑4820960
                                                                                                                                     |                                                                                                                                                    
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                                                                                          <  5:72966/149‑1 (MQ=255)
                                             ttGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTc                                                                                          >  2:130932/1‑149 (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                                                                                              <  2:189083/119‑1 (MQ=255)
                                                                       cAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCAc                                                                                              >  1:189083/1‑119 (MQ=255)
                                                                          ttCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAgg                                                             >  3:19776/1‑149 (MQ=255)
                                                                          ttCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAgg                                                             >  2:201920/1‑149 (MQ=255)
                                                                           tCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAgg                                                             <  6:133721/148‑1 (MQ=255)
                                                                             gcgcAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCAc                                                          >  8:38298/1‑149 (MQ=255)
                                                                             gcgcAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCAc                                                          >  6:51339/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                                     <  5:51339/148‑1 (MQ=255)
                                                                                                   atCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGc                                     <  4:78578/148‑1 (MQ=255)
                                                                                                      tcAAGGCGGGTATAACAATCGCCGACGACAATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTg                                 >  8:219/1‑149 (MQ=255)
                                                                                                        aaGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTg                                 <  3:187056/147‑1 (MQ=255)
                                                                                                                  taACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGtt                                                                                           >  8:214547/1‑79 (MQ=255)
                                                                                                                  taACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGtt                                                                                           <  7:214547/79‑1 (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                        <  6:205816/141‑1 (MQ=255)
                                                                                                                       aTCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAAcc                        >  5:205816/1‑141 (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)
                                                                                                                                     |                                                                                                                                                    
CTGTAGCGCGACAGCACCAGGGCCCAGGCCGCCTGGGCGAAAGTGTTGACGGTGAGCTGATGGGCCTGGGCCAGTTCGCGCAGACGCACGCCATCATTCATCTCAAGGCGGGTATAACAATCGCCGACGACCATGCCGCTACCGGCGTGGTCATGGCGCAGCGGGCGGTCGCTCGGGATGGCTGTGGCACGTTCGAAGCCGGCCAGGTTGGCCTGCCACCAGGCACGTGCGTCGTCCAGGCCCTGGCGTTGCAGCCAACCGATGTAATCGCGGTAACGCGGT  >  NC_002947/4820679‑4820960

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