Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,821,130 C→T 6.5% E3354K (GAG→AAG)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,821,1300CT 6.5% 232.3 / 3.5 77E3354K (GAG→AAG) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (35/37);  new base T (3/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

CCTTCGCCAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCATCGCCTCG  >  NC_002947/4820985‑4821275
                                                                                                                                                     |                                                                                                                                                 
ccTTCGCCAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCt                                                                                                                                                <  3:408782/149‑1 (MQ=255)
 cTTCGCCAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATAATAAGATGTGGTGGTTGCTAAACATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTCCC‑TCga                                                                                                                                                 >  1:353291/1‑148 (MQ=255)
 cTTCGCCAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCt                                                                                                                                                >  8:33149/1‑148 (MQ=255)
   tCGCCAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCt                                                                                                                                                >  5:429526/1‑146 (MQ=255)
       cAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAg                                                                                                                                         <  2:169933/149‑1 (MQ=255)
        aaGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAgg                                                                                                                                        <  8:359955/149‑1 (MQ=255)
        aaGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAgg                                                                                                                                        >  2:314000/1‑149 (MQ=255)
           ggCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TTGCTCAGCGAGGCac                                                                                                                                     <  8:124530/149‑2 (MQ=255)
             ccTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAg                                                                                                                                   <  1:399444/149‑1 (MQ=255)
                      ccTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTcg                                                                                                                          <  1:22308/149‑1 (MQ=255)
                      ccTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCGGAGAAAGGCCAAGCCGGCCTcg                                                                                                                          >  6:301731/1‑149 (MQ=255)
                       cTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGc                                                                                                                                            <  6:121224/130‑1 (MQ=255)
                       cTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGc                                                                                                                                            >  5:121224/1‑130 (MQ=255)
                       cTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTcgc                                                                                                                         >  2:112791/1‑149 (MQ=255)
                        tCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTcgcg                                                                                                                        >  2:251726/1‑149 (MQ=255)
                         cGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCgcgc                                                                                                                       >  4:354828/1‑149 (MQ=255)
                          gaagaaATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCt                                                                                                                      >  4:387491/1‑149 (MQ=255)
                          gaagaaATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCt                                                                                                                      >  2:19901/1‑149 (MQ=255)
                           aagaaATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTc                                                                                                                     <  4:85169/149‑1 (MQ=255)
                             gaaATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCt                                                                                                                   <  7:366382/149‑1 (MQ=255)
                               aaTCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGc                                                                                                                 <  5:375066/149‑1 (MQ=255)
                                    ttCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTtgtg                                                                                                            <  2:353291/149‑1 (MQ=255)
                                        tcaACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAgg                                                                                                        <  4:218666/149‑1 (MQ=255)
                                        tcaACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑AGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAgg                                                                                                        <  2:34789/149‑1 (MQ=255)
                                          aaCAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGc                                                                                                      <  6:52061/149‑1 (MQ=255)
                                           aCAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGcc                                                                                                     >  4:393874/1‑149 (MQ=255)
                                                   cgGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGgcg                                                                                             <  6:429526/149‑1 (MQ=255)
                                                    ggCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGgcgc                                                                                            <  1:251726/149‑1 (MQ=255)
                                                      cACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCgctg                                                                                          <  3:354828/149‑1 (MQ=255)
                                                         cAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTc                                                                                       >  6:27495/1‑149 (MQ=255)
                                                            gCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTg                                                                                    <  1:112791/149‑1 (MQ=255)
                                                            gCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTg                                                                                    >  1:383863/1‑149 (MQ=255)
                                                             cATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTgg                                                                                   <  2:5139/149‑1 (MQ=255)
                                                              aTCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTggg                                                                                  >  1:171605/1‑149 (MQ=255)
                                                               tCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑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                                                                         >  1:306450/1‑149 (MQ=255)
                                                                         gATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCgg                                                                         >  5:303567/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:401780/1‑149 (MQ=255)
                                                                              ggtggtTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAgg                                                                    >  6:143094/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                                                     <  1:19901/149‑1 (MQ=255)
                                                                                             tGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTgg                                                     <  5:27495/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             <  5:143094/149‑1 (MQ=255)
                                                                                                                                     gTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTgca             <  7:376171/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)
                                                                                                                                                     |                                                                                                                                                 
CCTTCGCCAAGGGCCTGGTAGACCTCGAAGAAATCGTTCATCAACAGCGAGCGGCACCAGGCATCGATA‑‑AGGATGTGGTGGTTGCTCATCATGAACCAGTAGCGTTCCTCGGCCACGCGCACCAGGCGCAGGTGGAAC‑GGTGCC‑TCGCTGAGCAAGGCAAAGCCGGCCTCGCGCTCCTGCTTGTGCAGGGCCTGCAGGCGCTGCTCCTGGGCATCGTCGGCCAGGCCGCGCCAGTCCTGGTAGTCCACCGCCAGGTTGCCCGGCTTGTGGATGATCTGCAGCATCGCCTCG  >  NC_002947/4820985‑4821275

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