Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,821,134 G→C 6.7% L3352L (CTC→CTG)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,821,1340GC 6.7% 242.5 / 6.5 75L3352L (CTC→CTG) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (32/38);  new base C (3/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

CAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCATCGCCTCGCCGCT  >  NC_002947/4820992‑4821280
                                                                                                                                                |                                                                                                                                                  
cAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAg                                                                                                                                              <  2:169933/149‑1 (MQ=255)
 aaGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAgg                                                                                                                                             <  8:359955/149‑1 (MQ=255)
 aaGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAgg                                                                                                                                             >  2:314000/1‑149 (MQ=255)
    ggCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TTGCTCAGCGAGGCac                                                                                                                                          <  8:124530/149‑2 (MQ=255)
      ccTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAg                                                                                                                                        <  1:399444/149‑1 (MQ=255)
               ccTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTcg                                                                                                                               <  1:22308/149‑1 (MQ=255)
               ccTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCGGAGAAAGGCCAAGCCGGCCTcg                                                                                                                               >  6:301731/1‑149 (MQ=255)
                cTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGc                                                                                                                                                 >  5:121224/1‑130 (MQ=255)
                cTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGc                                                                                                                                                 <  6:121224/130‑1 (MQ=255)
                cTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTcgc                                                                                                                              >  2:112791/1‑149 (MQ=255)
                 tCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTcgcg                                                                                                                             >  2:251726/1‑149 (MQ=255)
                  cGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCgcgc                                                                                                                            >  4:354828/1‑149 (MQ=255)
                   gaagaaATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCt                                                                                                                           >  4:387491/1‑149 (MQ=255)
                   gaagaaATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCt                                                                                                                           >  2:19901/1‑149 (MQ=255)
                    aagaaATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTc                                                                                                                          <  4:85169/149‑1 (MQ=255)
                      gaaATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCt                                                                                                                        <  7:366382/149‑1 (MQ=255)
                        aaTCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGc                                                                                                                      <  5:375066/149‑1 (MQ=255)
                             ttCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTtgtg                                                                                                                 <  2:353291/149‑1 (MQ=255)
                                 tcaACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAgg                                                                                                             <  4:218666/149‑1 (MQ=255)
                                 tcaACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑AGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAgg                                                                                                             <  2:34789/149‑1 (MQ=255)
                                   aaCAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGc                                                                                                           <  6:52061/149‑1 (MQ=255)
                                    aCAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGcc                                                                                                          >  4:393874/1‑149 (MQ=255)
                                            cgGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGgcg                                                                                                  <  6:429526/149‑1 (MQ=255)
                                             ggCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGgcgc                                                                                                 <  1:251726/149‑1 (MQ=255)
                                               cACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCgctg                                                                                               <  3:354828/149‑1 (MQ=255)
                                                  cAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTc                                                                                            >  6:27495/1‑149 (MQ=255)
                                                     gCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTg                                                                                         >  1:383863/1‑149 (MQ=255)
                                                     gCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTg                                                                                         <  1:112791/149‑1 (MQ=255)
                                                      cATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTgg                                                                                        <  2:5139/149‑1 (MQ=255)
                                                       aTCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTggg                                                                                       >  1:171605/1‑149 (MQ=255)
                                                        tCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGc                                                                                      <  3:393874/149‑1 (MQ=255)
                                                              aGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCAtcgtc                                                                                <  5:378068/149‑1 (MQ=255)
                                                              aGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGACGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCAtcgtc                                                                                >  7:366448/1‑149 (MQ=255)
                                                               ggATCTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGAGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCAtcgtcg                                                                               >  4:175717/1‑149 (MQ=255)
                                                                gATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCgg                                                                              >  5:303567/1‑149 (MQ=255)
                                                                gATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCgg                                                                              >  1:306450/1‑149 (MQ=255)
                                                                   gtggtgGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGcc                                                                            <  3:239403/148‑1 (MQ=255)
                                                                   gtggtgGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGcc                                                                            >  4:239403/1‑148 (MQ=255)
                                                                    tggtggTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCa                                                                           <  3:387491/148‑1 (MQ=255)
                                                                     ggtggtTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAgg                                                                         >  6:143094/1‑149 (MQ=255)
                                                                     ggtggtTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAgg                                                                         >  6:401780/1‑149 (MQ=255)
                                                                     ggtggtTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAgg                                                                         >  1:113018/1‑149 (MQ=255)
                                                                       tggtTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGcc                                                                       >  6:361661/1‑149 (MQ=255)
                                                                         gtTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCcgc                                                                     >  4:223700/1‑149 (MQ=255)
                                                                              tcatcaTGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAg                                                                >  3:425176/1‑149 (MQ=255)
                                                                               catcatGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAACAGGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAg                                                                >  1:325651/1‑149 (MQ=255)
                                                                                atcatGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTc                                                              >  2:187855/1‑149 (MQ=255)
                                                                                 tcatGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTcc                                                             >  6:41014/1‑149 (MQ=255)
                                                                                   atGAACCAGTAGCGTTCCTCTGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTg                                                           <  6:303567/149‑1 (MQ=255)
                                                                                   atGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTg                                                           <  4:295432/149‑1 (MQ=255)
                                                                                    tGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCCTTCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTg                                                           >  1:114414/1‑149 (MQ=255)
                                                                                    tGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTgg                                                          <  5:27495/149‑1 (MQ=255)
                                                                                    tGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTgg                                                          <  1:19901/149‑1 (MQ=255)
                                                                                           gTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCgctg                                                                                               >  1:3348/1‑105 (MQ=255)
                                                                                           gTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCgctg                                                                                               <  2:3348/105‑1 (MQ=255)
                                                                                           gTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCa                                                   >  4:170089/1‑149 (MQ=255)
                                                                                                 ttCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCa                                             <  5:296475/149‑1 (MQ=255)
                                                                                                   ccTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAg                                            <  3:223700/148‑1 (MQ=255)
                                                                                                     tCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGtt                                         >  8:376171/1‑149 (MQ=255)
                                                                                                        gCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGcc                                      <  5:41014/149‑1 (MQ=255)
                                                                                                         ccACGCTCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGccc                                     <  6:255542/149‑1 (MQ=255)
                                                                                                             gcgcACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGTAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCt                                 >  5:353075/1‑149 (MQ=255)
                                                                                                                cACCAGTCGCAGGGTGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCAGGCTtgt                              <  8:366448/149‑1 (MQ=255)
                                                                                                                   cAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGga                           <  7:33149/149‑1 (MQ=255)
                                                                                                                            gTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTgca                  <  4:393117/149‑1 (MQ=255)
                                                                                                                            gTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTgca                  <  7:376171/149‑1 (MQ=255)
                                                                                                                            gTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTgca                  <  5:143094/149‑1 (MQ=255)
                                                                                                                                  gctttGCC‑TTGCTCAGCGAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTtgtg                             >  7:153119/5‑133 (MQ=255)
                                                                                                                                  gctttGCC‑TTGCTCAGCGAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTtgtg                             <  8:153119/129‑1 (MQ=255)
                                                                                                                                  gctttGCC‑TTGCTCAGCGAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCAtc             >  1:423106/5‑149 (MQ=255)
                                                                                                                                    ggTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGTCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCAtcgc           <  3:175717/149‑1 (MQ=255)
                                                                                                                                     ttGCC‑TTGCTCAGCGAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCAtcgcc          >  2:329172/2‑149 (MQ=255)
                                                                                                                                        cc‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCAtcgcctcg       >  8:78898/1‑149 (MQ=255)
                                                                                                                                              cTGTGCAAGCCAAAGCCGGCCTCGCGCTCGTGCTTGTGCAGGGCCTGCAGGGGCTGCTCCTGGGCATCGTCGGCCAGGCGGCGCCTGTCCTGGTAGTCCACCTCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCATCGCCTCGCCGCt  <  1:314000/149‑1 (MQ=255)
                                                                                                                                              atGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCATCGCCTcgccgc   <  5:313909/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCATCGCCTCGCCGCT  >  NC_002947/4820992‑4821280

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