Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,821,138 A→G 6.7% L3351S (TTG→TCG)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,821,1380AG 6.7% 209.6 / 4.4 74L3351S (TTG→TCG) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base G (3/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

CAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCATCGCCTCGCCGCTGTTCC  >  NC_002947/4820992‑4821285
                                                                                                                                                    |                                                                                                                                                   
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‑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                                                                              >  1:113018/1‑149 (MQ=255)
                                                                     ggtggtTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAgg                                                                              >  6:401780/1‑149 (MQ=255)
                                                                     ggtggtTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAgg                                                                              >  6:143094/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                                                               <  1:19901/149‑1 (MQ=255)
                                                                                    tGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTgg                                                               <  5:27495/149‑1 (MQ=255)
                                                                                           gTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCgctg                                                                                                    <  2:3348/105‑1 (MQ=255)
                                                                                           gTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCgctg                                                                                                    >  1:3348/1‑105 (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                       <  5:143094/149‑1 (MQ=255)
                                                                                                                            gTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTgca                       <  7:376171/149‑1 (MQ=255)
                                                                                                                            gTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTgca                       <  4:393117/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)
                                                                                                                                                   cAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCATCGCCTCGCCGCTGTTcc  >  7:398091/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAAGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCATCGCCTCGCCGCTGTTCC  >  NC_002947/4820992‑4821285

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