Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 312,717 T→C 5.2% Q166R (CAG→CGG)  fdhD ← protein sulfur transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947312,7170TC 5.2% 237.9 / 2.3 77Q166R (CAG→CGG) fdhDprotein sulfur transferase
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base C (2/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

CAGGCCACCCTGGGTGTCGATGCGTTGTCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGACCCCACACAGGCCGCAACCACTGGTGCCGGCC  >  NC_002947/312572‑312866
                                                                                                                                                  |                                                                                                                                                     
caggccaCCCTGGGTGTCGATGCGTTGTCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCC                                                                                                                                                   >  2:502047‑M1/1‑150 (MQ=255)
   gccaCCCTGGGTGTCGATGCGTTGTCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAG                                                                                                                                                >  2:170226‑M1/1‑150 (MQ=255)
      aCCCTGGGTGTCGATGCGTTGTCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGG                                                                                                                                             >  1:299717‑M1/1‑150 (MQ=255)
       cccTGGGTGTCGATGCGTTGTCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGT                                                                                                                                               <  2:626644‑M1/147‑1 (MQ=255)
       cccTGGGTGTCGATGCGTTGTCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGGCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGT                                                                                                                                               >  1:626644‑M1/1‑147 (MQ=255)
                    tGCGTTGTCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGG                                                                                                                                             <  2:243646‑M1/136‑1 (MQ=255)
                    tGCGTTGTCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGG                                                                                                                                             >  1:243646‑M1/1‑136 (MQ=255)
                        ttGTCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGaa                                                                                                                                       >  2:899639/1‑138 (MQ=255)
                        ttGTCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGaa                                                                                                                                       <  1:899639/138‑1 (MQ=255)
                        ttGTCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCCCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATgcg                                                                                                                           >  2:97408/1‑150 (MQ=255)
                         tGTCGCAGCAGGGCGCCGCTCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATgcgc                                                                                                                          <  1:365540/150‑1 (MQ=255)
                           tCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGtcggcctgcgcgagcgc                                                                                                                        >  1:122070‑M1/1‑133 (MQ=255)
                                   gggCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGc                                                                                                                <  1:97408/150‑1 (MQ=255)
                                        ccGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAcc                                                                                                           <  1:587319‑M2/150‑3 (MQ=255)
                                         cGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTGTTCGAACGGAATCAGCAGTTCGCCTTGGCGGTCAATAAGCAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAcca                                                                                                          >  1:870629‑M2/1‑147 (MQ=255)
                                           aTCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGtcggcctgcgcgagcgcatcgacgccttccagc                                                                                                        <  2:299505‑M1/150‑34 (MQ=255)
                                            tCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTcgc                                                                                                                      <  2:638538/135‑1 (MQ=255)
                                            tCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTcgc                                                                                                                      >  1:638538/1‑135 (MQ=255)
                                                tttGTCGAGGGCGTTGTGCCGGCCGATGTCATAGAGGCACTCGTGACTGGCGTCCTGGCGTTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTg                                                                                                   <  1:1173330/150‑1 (MQ=255)
                                                   gTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGtcggcctgcgcgagcgcatcgacgccttccagccactgggc                                                                                                <  2:122070‑M1/150‑42 (MQ=255)
                                                     cGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGtcggcctgcgcgagcgcatcgacgccttccagccactgggcca                                                                                              <  1:170226‑M1/150‑44 (MQ=255)
                                                      gAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAccagccactgg                                                                                                  <  2:738679‑M2/145‑12 (MQ=255)
                                                      gAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAccagccactgg                                                                                                  >  1:738679‑M2/1‑134 (MQ=255)
                                                      gAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGAATGAACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCgg                                                                                             >  2:1002952/1‑150 (MQ=255)
                                                       aGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCGGGCCGACCAGCCACTGGGCCggg                                                                                            <  2:299717/150‑1 (MQ=255)
                                                                tgCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGTCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGgc                                                                                   <  1:502047/150‑1 (MQ=255)
                                                                       cgATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCg                                                                            >  1:566949/1‑150 (MQ=255)
                                                                       cgATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCACGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCg                                                                            >  2:1103458/1‑150 (MQ=255)
                                                                           gTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTTGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCaa                                                                        >  1:389901/1‑150 (MQ=255)
                                                                              ttCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCCGGCCCAGTGGCTGGtcggcctgcgcgagcgcatcgacgccttccagccactgggccag                                                                                             >  1:342002‑M1/1‑82 (MQ=255)
                                                                              ttCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCCGGCCCAGTGGCTGGtcggcctgcgcgagcgcatcgacgccttccagccactgggccag                                                                                             <  2:342002‑M1/126‑45 (MQ=255)
                                                                                   ggCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCa                                                                >  2:392782/1‑150 (MQ=255)
                                                                                      cAAGCAGCAGTTCGCCTTGGCGGTCCATGAAGAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAccagc                                                                                                        <  2:741838‑M2/107‑6 (MQ=255)
                                                                                      cAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAccagc                                                                                                        >  1:741838‑M2/1‑102 (MQ=255)
                                                                                       aaGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATc                                                            >  1:903345/1‑150 (MQ=255)
                                                                                        agcagcAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATcc                                                           <  2:903345/150‑1 (MQ=255)
                                                                                        agcagcAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATcc                                                           >  1:846991/1‑150 (MQ=255)
                                                                                         gcagcaGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAccagccactg                                                                                                   >  2:821750‑M2/1‑99 (MQ=255)
                                                                                         gcagcaGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAccagccactg                                                                                                   <  1:821750‑M2/109‑11 (MQ=255)
                                                                                         gcagcaGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGc                                                                  >  1:608982/1‑142 (MQ=255)
                                                                                         gcagcaGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGc                                                                  <  2:608982/142‑1 (MQ=255)
                                                                                         gcagcaGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCg                                                          >  1:654371/1‑150 (MQ=255)
                                                                                               gTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGt                                                     <  2:7041/149‑1 (MQ=255)
                                                                                               gTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGt                                                     >  1:7041/1‑149 (MQ=255)
                                                                                                    gCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGccgg                                                <  2:1046792/150‑3 (MQ=255)
                                                                                                      cTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCt                                              >  2:988854/1‑150 (MQ=255)
                                                                                                           cGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGt                                         <  2:846991/150‑1 (MQ=255)
                                                                                                                 aTGAACAATGCTGCATGCACGGCACCGCAATGCCGGCCCAGTGGCTGGtcggcctgcgcgagcgcatcgacgccttccagccactg                                                                                                   <  2:73944‑M1/86‑39 (MQ=255)
                                                                                                                 aTGAACAATGCTGCATGCACGGCACCGCAATGCCGGCCCAGTGGCTGGtcggcctgcgcgagcgcatcgacgccttccagccactg                                                                                                   >  1:73944‑M1/1‑48 (MQ=255)
                                                                                                                   gAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGAc                                 <  2:566949/150‑1 (MQ=255)
                                                                                                                    aaCAATGCTGCATGCACGGCACCGCCCTCCTGTCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGAcc                                >  1:338024/1‑150 (MQ=255)
                                                                                                                    aaCAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGc                                                                  <  2:688037/116‑1 (MQ=255)
                                                                                                                    aaCAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGc                                                                  >  1:688037/1‑116 (MQ=255)
                                                                                                                    aaCAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGt                                                     <  2:126342/129‑1 (MQ=255)
                                                                                                                    aaCAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGt                                                     >  1:126342/1‑129 (MQ=255)
                                                                                                                     aCAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGt                                         >  1:536721/1‑140 (MQ=255)
                                                                                                                     aCAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGt                                         <  2:536721/140‑1 (MQ=255)
                                                                                                                     aCAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGAccc                               <  2:1137395/150‑1 (MQ=255)
                                                                                                                      cAATGCTGCATGCACGGAACCGCAATGCTGGCCCAGTGGCAGGAAGGCGTCGATGCGTGCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCCCCCGGCACCACGGCCCGATCCGCCAGTGCCTGTTCCATTTCCTCGAcccc                              <  2:865074/150‑1 (MQ=255)
                                                                                                                       aaTGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGACCCca                             >  1:144133/1‑150 (MQ=255)
                                                                                                                           ctgcATGCACGGCACCGCAATGCTGGCCCAGTAGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGACCCcccaca                         >  1:243421/1‑150 (MQ=255)
                                                                                                                            tgcatgcaCGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGAt                                                                                                                              <  2:146786/48‑1 (MQ=255)
                                                                                                                            tgcatgcaCGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGAt                                                                                                                              >  1:146786/1‑48 (MQ=255)
                                                                                                                              catgcaCGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGACCCCACACAGGc                      <  1:988854/150‑1 (MQ=255)
                                                                                                                                   aCGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGcc                                      >  2:1169576/1‑129 (MQ=255)
                                                                                                                                   aCGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGcc                                      <  1:1169576/129‑1 (MQ=255)
                                                                                                                                   aCGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGACCCCACACAGGCCGCaa                 >  1:1192359/1‑150 (MQ=255)
                                                                                                                                     ggCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAccagccactggg                                                                                                 >  2:732241‑M2/1‑56 (MQ=255)
                                                                                                                                     ggCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAccagccactggg                                                                                                 <  1:732241‑M2/68‑13 (MQ=255)
                                                                                                                                     ggCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCt                                                 <  1:1189685/116‑1 (MQ=255)
                                                                                                                                     ggCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCt                                                 >  2:1189685/1‑116 (MQ=255)
                                                                                                                                             aaTGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGACCCCACACAGGCCGCAACCACTGGTg        <  2:1115003/149‑1 (MQ=255)
                                                                                                                                             aaTGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGACCCCACACAGGCCGCAACCACTGGTg        >  1:1115003/1‑149 (MQ=255)
                                                                                                                                                gCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGACCCCACACAGGCCGCAACCACTGGTGCCgg    <  1:1002952/150‑1 (MQ=255)
                                                                                                                                                 cTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAccagccactgggccggcattgcggtgccgtgcatgcagcattgttcatggaccgccaaggcgaactgctgcttggccgcgaa                           >  1:342002‑M2/1‑44 (MQ=255)
                                                                                                                                                 cTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAccagccactgggccggcattgcggtgccgtgcatgcagcattgttcatggaccgccaaggcgaactgctgcttggccgcgaa                           <  2:342002‑M2/126‑83 (MQ=255)
                                                                                                                                                 cTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGACCCCACACAGGCCGGCACCAGTGGttgcggc   <  1:392782/150‑7 (MQ=255)
                                                                                                                                                  tGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGAccagccactgggccagcattgcggtgccgtgcatgcagcattgttcatggaccgccaaggcgaactgctgcttggccgcgaagacatcggccggcacaacgccctcg  <  1:170226‑M2/150‑108 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CAGGCCACCCTGGGTGTCGATGCGTTGTCGCAGCAGGGCGCCGATCAGTTTGTCGAGGGCGTTGTGCCGGCCGATGTCTTCGCGGCCAAGCAGCAGTTCGCCTTGGCGGTCCATGAACAATGCTGCATGCACGGCACCGCAATGCTGGCCCAGTGGCTGGAAGGCGTCGATGCGCTCGCGCAGGCCGACCAGCCACTGGGCCGGGGGCAGGGGCGCACCCGGCAACACGGCCAGATCCGGCAGTGCCTGCTCCAGTGCCTCGACCCCACACAGGCCGCAACCACTGGTGCCGGCC  >  NC_002947/312572‑312866

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: