Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,824,152 C→T 7.8% E2346E (GAG→GAA pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,824,1520CT 7.8% 181.1 / 4.3 64E2346E (GAG→GAApvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base T (3/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAG  >  NC_002947/4824008‑4824299
                                                                                                                                                |                                                                                                                                                   
gcgcAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGc                                                                                                                                                 <  8:201474/149‑1 (MQ=255)
gcgcAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGc                                                                                                                                                 <  4:121/149‑1 (MQ=255)
  gcAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAGTTCGCCa                                                                                                                                               <  7:276310/149‑1 (MQ=255)
      tACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGc                                                                                                                                           >  2:106891/1‑149 (MQ=255)
        ccACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCt                                                                                                                                          >  3:299752/1‑148 (MQ=255)
        ccACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCt                                                                                                                                          <  8:251929/148‑1 (MQ=255)
         cACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGa                                                                                                                                        <  1:106891/149‑1 (MQ=255)
                    cacacaACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGccc                                                                                                                             >  2:80814/1‑149 (MQ=255)
                    cacacaACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGccc                                                                                                                             >  6:5734/1‑149 (MQ=255)
                     acacaACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCt                                                                                                                            >  8:126254/1‑149 (MQ=255)
                           caGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGACTGACCACCACGCCCTTGGGtt                                                                                                                      >  1:267700/1‑149 (MQ=255)
                             gCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTg                                                                                                                    >  8:252467/1‑149 (MQ=255)
                               ggCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGCGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGcc                                                                                                                  <  2:302212/149‑1 (MQ=255)
                                     cacaGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGa                                                                                                            <  6:146509/149‑1 (MQ=255)
                                      acaGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGttt                                                                                                                     <  6:227216/139‑1 (MQ=255)
                                      acaGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGttt                                                                                                                     >  5:227216/1‑139 (MQ=255)
                                      acaGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAg                                                                                                           >  1:138855/1‑149 (MQ=255)
                                       caGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGc                                                                                                          >  7:278326/1‑149 (MQ=255)
                                         ggCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCt                                                                                                        <  7:252467/149‑1 (MQ=255)
                                         ggCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGACTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCt                                                                                                        <  2:267700/149‑1 (MQ=255)
                                            gctcgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAg                                                                                                     <  8:278326/149‑1 (MQ=255)
                                             ctcgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAgg                                                                                                    >  1:26183/1‑149 (MQ=255)
                                              tcgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGATGCGCATGACGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCAATGGCTGACAACCACGCCCTTGGGTTTGCCGGTGGAGCCTCAGgt                                                                                                   >  7:114939/1‑149 (MQ=255)
                                              ccgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAgg                                                                                                    <  4:163975/147‑1 (MQ=255)
                                               cgctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGgtg                                                                                                  >  5:148230/1‑149 (MQ=255)
                                                 ctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTa                                                                                                <  2:26183/149‑1 (MQ=255)
                                                 ctGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTa                                                                                                <  5:100455/149‑1 (MQ=255)
                                                      aGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCa                                                                                           <  6:110084/149‑1 (MQ=255)
                                                       gCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAg                                                                                          <  6:158668/149‑1 (MQ=255)
                                                             aaaTTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGc                                                                                    <  1:80814/149‑1 (MQ=255)
                                                              aaTTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGc                                                                                    <  2:115232/148‑1 (MQ=255)
                                                               aTTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAGTTCGCCATgcact                                                                                                                                         <  7:284932/94‑5 (MQ=255)
                                                               aTTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAGTTCGCCATgcact                                                                                                                                         >  8:284932/1‑90 (MQ=255)
                                                                  gATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGt                                                                               >  5:131922/1‑149 (MQ=255)
                                                                   aTCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGgtgt                                                                                                 >  8:28637/1‑130 (MQ=255)
                                                                   aTCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGgtgt                                                                                                 <  7:28637/130‑1 (MQ=255)
                                                                       aaTAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGtgctgc                                                                          <  5:5734/149‑1 (MQ=255)
                                                                       aaTAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGtgct                                                                            >  2:236003/1‑147 (MQ=255)
                                                                            aaaTGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAg                                                                     <  1:236003/149‑1 (MQ=255)
                                                                              aTGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGc                                                                   >  3:288416/1‑149 (MQ=255)
                                                                               tGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCt                                                                  >  5:30581/1‑149 (MQ=255)
                                                                                  aGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAg                                                               <  8:11589/149‑1 (MQ=255)
                                                                                  aGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAg                                                               >  3:27184/1‑149 (MQ=255)
                                                                                  aGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAg                                                               <  3:59704/149‑1 (MQ=255)
                                                                                       gCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTg                                                           <  1:215255/148‑1 (MQ=255)
                                                                                         aaTCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTcc                                                        <  7:126254/149‑1 (MQ=255)
                                                                                                gCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGAAGGCTGAGGTTGTCCAGCGGTg                                                 >  8:255147/1‑149 (MQ=255)
                                                                                                     gcATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCg                                            <  2:138855/149‑1 (MQ=255)
                                                                                                      cATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGc                                           >  1:54972/1‑149 (MQ=255)
                                                                                                          ccGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTg                                         >  5:20218/1‑147 (MQ=255)
                                                                                                          ccGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTg                                         >  1:169977/1‑147 (MQ=255)
                                                                                                          ccGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGt                                        <  4:27184/148‑1 (MQ=255)
                                                                                                                  ggcgATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGgccgc                                 >  5:176564/1‑147 (MQ=255)
                                                                                                                         aCCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCTGCGGTGCATCGCTGTAGGCCGCCAGGCTgg                         <  1:25305/148‑1 (MQ=255)
                                                                                                                          ccGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTaa                                                                                      <  6:236212/86‑1 (MQ=255)
                                                                                                                          ccGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTaa                                                                                      >  5:236212/1‑86 (MQ=255)
                                                                                                                            gCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTggcgg                      >  1:192530/1‑148 (MQ=255)
                                                                                                                                 gCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCACGTAAGACTGGTGCTCAGTAAAAATAAAGTTGTCAAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCg                 >  4:81306/1‑148 (MQ=25)
                                                                                                                                     tGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTc            <  6:148230/149‑1 (MQ=255)
                                                                                                                                      ccATGGCGAATTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCACGTAAGACTCCTTATCATTAAAAAAAAAAATGTCCATCGGTGAATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTc            >  6:111865/2‑148 (MQ=255)
                                                                                                                                       cATGGCGAGTTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCACGTAAGACTCCTTATCATTAAAAAAAAAATTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTc            >  6:6304/1‑147 (MQ=255)
                                                                                                                                           gCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCtt             <  2:210437/142‑1 (MQ=255)
                                                                                                                                           gCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTAGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCtt             >  1:210437/1‑142 (MQ=255)
                                                                                                                                               aCTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAg  <  6:176564/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCGCAGTACCACGCGCGCACCACACAACAGCGGCACCCACAGGCGCTCGCTGGCAGCGTCGAAATTGATCGAATAGAAATGCAGCTCGCAATCATCGCTGCGCATGCCGAATGCGGCGATCACCGCCTGGCAGTGCATGGCGAACTCGCCATGGCTGACCACCACGCCCTTGGGTTTGCCGGTGGAGCCTGAGGTGTAGATCAGGTAAGCCTGGTGCTGCGGCAGGCTGAGGTTGTCCAGCGGTGCATCGCTGTAGGCCGCCAGGCTGGCGGCATCGTCTTCCAGGCTCCAG  >  NC_002947/4824008‑4824299

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: