Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,821,760 C→A 6.1% E3144* (GAG→TAG)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,821,7600CA 6.1% 212.2 / 3.0 66E3144* (GAG→TAG) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base A (2/2);  total (34/32)
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.68e-01

CCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACGCCGTAGCCCAGGCCCTT  >  NC_002947/4821616‑4821907
                                                                                                                                                |                                                                                                                                                    
ccAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAgg                                                                                                                                                  <  4:100596/149‑1 (MQ=255)
 cAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGAAGAGCACCGCCATACACCTGGCCATCGAAGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAggg                                                                                                                                                 >  4:289751/1‑149 (MQ=255)
      gcgcgATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCTCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTgg                                                                                                                                            <  8:111627/149‑1 (MQ=255)
      gcgcgATGGTCTGCGGGTCGTAGCGCGCGGCGCTGAAGGTCCAGCGCAGGTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTgg                                                                                                                                            <  8:47784/149‑1 (MQ=255)
         cgATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAAc                                                                                                                                         <  8:387291/149‑1 (MQ=255)
         cgATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTACAGGGGCTGGaa                                                                                                                                          >  7:173056/1‑148 (MQ=255)
         cgATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTACAGGGGCTGGaa                                                                                                                                          >  7:178286/1‑148 (MQ=255)
          gATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACa                                                                                                                                        <  5:294406/149‑1 (MQ=255)
                   ctgGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCg                                                                                                                                  <  5:277786/146‑1 (MQ=255)
                   ctgGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCAc                                                                                                                               >  7:233422/1‑149 (MQ=255)
                       tCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGATCAAGGCCGCACTGTc                                                                                                                           <  3:271716/149‑1 (MQ=255)
                          tAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGTAACGGCGCATCGATGTCATGCGCCAACCTGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGaa                                                                                                                        <  5:284562/149‑1 (MQ=255)
                            gcgcTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAAct                                                                                                                      >  3:33689/1‑149 (MQ=255)
                              gcTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAActgc                                                                                                                    <  4:110265/149‑1 (MQ=255)
                                tCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAActgctg                                                                                                                  >  8:247954/1‑149 (MQ=255)
                                  ggCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGt                                                                                                                >  5:346463/1‑149 (MQ=255)
                                  ggCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGt                                                                                                                >  8:185132/1‑149 (MQ=255)
                                        gAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACt                                                                                                          >  7:32736/1‑149 (MQ=255)
                                        gAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACt                                                                                                          >  3:210367/1‑149 (MQ=255)
                                         aaGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTACAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAAc                                                                                                           <  8:173056/147‑1 (MQ=255)
                                         aaGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTACAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAAc                                                                                                           <  8:178286/147‑1 (MQ=255)
                                           ggTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGc                                                                                                       >  5:225530/1‑149 (MQ=255)
                                             tCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGccc                                                                                                     <  7:185132/149‑1 (MQ=255)
                                               cAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAg                                                                                                   >  1:329170/1‑149 (MQ=255)
                                                aGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTgg                                                                                                        >  8:199977/1‑143 (MQ=255)
                                                aGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTgg                                                                                                        <  7:199977/143‑1 (MQ=255)
                                                aGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAgg                                                                                                  <  6:225530/149‑1 (MQ=255)
                                                   gcAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAg                                                                                               <  4:121763/149‑1 (MQ=255)
                                                   gcAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAg                                                                                               <  8:233422/149‑1 (MQ=255)
                                                    cAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGt                                                                                              <  6:351967/149‑1 (MQ=255)
                                                          cAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAgg                                                                                        <  7:247954/149‑1 (MQ=255)
                                                          cAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAgg                                                                                        >  4:298286/1‑149 (MQ=255)
                                                          cAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAg                                                                                         <  2:329170/148‑1 (MQ=255)
                                                           agagCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGt                                                                                       >  7:236846/1‑149 (MQ=255)
                                                              gCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGAt                                                                                    <  5:416648/149‑1 (MQ=255)
                                                               cACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATg                                                                                   >  2:65501/1‑149 (MQ=255)
                                                                          acCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCgg                                                                        >  3:78049/1‑149 (MQ=255)
                                                                                      aCGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCAt                                                            <  4:210367/149‑1 (MQ=255)
                                                                                        gCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCt                                                          >  7:73358/1‑149 (MQ=255)
                                                                                          tCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGc                                                        >  1:161985/1‑149 (MQ=255)
                                                                                            aaCCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTc                                                      <  2:188692/149‑1 (MQ=255)
                                                                                            aaCCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTc                                                      >  6:340717/1‑149 (MQ=255)
                                                                                            aaCCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTc                                                      <  4:33689/149‑1 (MQ=255)
                                                                                                  ttGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGccc                                                >  7:194176/1‑149 (MQ=255)
                                                                                                    gggCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCgcgc                                                  >  1:98618/1‑145 (MQ=255)
                                                                                                    gggCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCgcgc                                                  <  2:98618/145‑1 (MQ=255)
                                                                                                               cATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGGTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGc                                   >  2:165367/1‑148 (MQ=255)
                                                                                                                       tCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAgcg                          >  6:368394/1‑149 (MQ=255)
                                                                                                                       tCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAgcg                          <  3:298286/149‑1 (MQ=255)
                                                                                                                         aTGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCTCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAgcg                          >  4:111779/1‑147 (MQ=255)
                                                                                                                         aTGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAgcg                          >  4:19760/1‑147 (MQ=255)
                                                                                                                           gcgcCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAg                       >  3:175300/1‑148 (MQ=255)
                                                                                                                             gcCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTAc                    >  6:89493/1‑149 (MQ=255)
                                                                                                                               cAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACGc                  <  8:73358/149‑1 (MQ=255)
                                                                                                                                aaCCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACg                   >  4:258863/1‑147 (MQ=255)
                                                                                                                                  cccGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACGCCg                >  5:17451/1‑148 (MQ=255)
                                                                                                                                   ccGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACGCCg                >  5:293938/1‑147 (MQ=255)
                                                                                                                                   ccGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACGCCGTa              <  4:78049/149‑1 (MQ=255)
                                                                                                                                     ggCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACGCCGTAGc            <  6:204346/149‑1 (MQ=255)
                                                                                                                                      gCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACGCCGTAGcc           <  4:175300/149‑1 (MQ=255)
                                                                                                                                          aGGCCTCCAGGGGCTGGAACAAGGCCGCTCTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACGCCGTAGCCCAgg       <  8:194176/149‑1 (MQ=255)
                                                                                                                                             ccTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACGCCGTAGCCCAGGccc    <  3:258863/149‑1 (MQ=255)
                                                                                                                                              cTCCAGGGGCTGGTACAAGGCCGCACTGTCGAACTGCTTGTCGTACTGGCCCAGGTACTTGAATGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCTGCCAGGTAGCGCAGTACGCCGTAGCCCAGGccc    <  1:254330/149‑1 (MQ=255)
                                                                                                                                              cTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACGCCGTAGCCCAGGccc    >  2:254330/1‑148 (MQ=255)
                                                                                                                                              cTCCAGGGGCTGGAAAAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCGCGCGCCCGGCGGCATCCGCAAGGTAGCGAAGTACGCCGTAGCCCATGCCCt   >  3:332602/1‑149 (MQ=255)
                                                                                                                                               tCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCACAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACGCCGTAGACCAGGCCCtt  >  7:231374/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCAGACGCGCGATGGTCTGCTGGTCGTAGCGCTCGGCGCTGAAGGTCCAGCGCAGTTGCAGAGCACCGCCATACACCTGGCCATCGACGCTCAACCAGTTGGGCAACGGCGCATCGAGGTCATGCGCCAACCCGGCCGAGGCCTCCAGGGGCTGGAACAAGGCCGCACTGTCGAACTGCTGGTCGAACTGGCCCAGGTAGTTGAAGGTGATGCGCGCCTGCGGCAAGGCCGCCATCTGCTCGCGCCCGGCGGCATCCGCCAGGTAGCGCAGTACGCCGTAGCCCAGGCCCTT  >  NC_002947/4821616‑4821907

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: