Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,725,658 C→A 7.9% P436P (CCG→CCT)  opgH ← glucans biosynthesis glucosyltransferase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,725,6580CA 7.9% 186.2 / 5.4 63P436P (CCG→CCT) opgHglucans biosynthesis glucosyltransferase H
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base A (2/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

TCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCA  >  NC_002947/5725512‑5725805
                                                                                                                                                  |                                                                                                                                                   
tCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACgg                                                                                                                                                   >  4:169078/1‑149 (MQ=255)
   tCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGCAAAGgcgc                                                                                                                                                <  7:252620/149‑1 (MQ=255)
   tCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGACGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGgcgc                                                                                                                                                >  6:284071/1‑149 (MQ=255)
     aaGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGcc                                                                                                                                               <  4:57341/148‑1 (MQ=255)
      aGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAg                                                                                                                                             >  6:138252/1‑149 (MQ=255)
         aGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGCAAAGgcgcgt                                                                                                                                              >  6:61610/1‑143 (MQ=255)
         aGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCGTTGCCGGGCAAAGgcgcgt                                                                                                                                              <  5:61610/145‑3 (MQ=255)
           gTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGgcg                                                                                                                                         >  6:24893/1‑148 (MQ=255)
           gTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGAACCGGCGAACGCGCCTTTGCCCGGCACCGGCGCCAGGgcga                                                                                                                                        >  6:301829/1‑148 (MQ=255)
            tttGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCa                                                                                                                                       <  5:294769/149‑1 (MQ=255)
                  ggcAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTc                                                                                                                                 >  2:394257/1‑149 (MQ=255)
                  ggcAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTc                                                                                                                                 >  8:340703/1‑149 (MQ=255)
                    cAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGa                                                                                                                               <  1:95454/149‑1 (MQ=255)
                     aGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGAt                                                                                                                              <  3:241392/149‑1 (MQ=255)
                         cttcGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGaa                                                                                                                           <  7:340703/148‑1 (MQ=255)
                          ttcGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACg                                                                                                                         >  7:78677/1‑149 (MQ=255)
                            cGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACgg                                                                                                                        <  4:66505/148‑1 (MQ=255)
                             gTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACTCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCt                                                                                                                      <  3:169078/149‑1 (MQ=255)
                              tAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCtt                                                                                                                     <  6:202993/149‑1 (MQ=255)
                                    cccGGCAAGTCGTAGTCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGc                                                                                                                <  4:62206/148‑1 (MQ=255)
                                               gTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGt                                                                                                    <  5:284071/149‑1 (MQ=255)
                                                tAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGtt                                                                                                   >  3:78590/1‑149 (MQ=255)
                                                  ggCAATCCACACCCCCCAGCCGTCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTtgt                                                                                                 <  6:269918/149‑1 (MQ=255)
                                                   gCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTtgtg                                                                                                <  3:357390/149‑1 (MQ=255)
                                                     aaTCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGga                                                                                                             >  7:138079/1‑134 (MQ=255)
                                                     aaTCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGga                                                                                                             <  8:138079/134‑1 (MQ=255)
                                                        ccACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGccc                                                                                            >  1:295756/1‑148 (MQ=255)
                                                          acacacCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCag                                                                                         >  6:366715/1‑149 (MQ=255)
                                                                cccAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAAc                                                                                                                          <  7:58817/110‑1 (MQ=255)
                                                                cccAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAAc                                                                                                                          >  8:58817/1‑110 (MQ=255)
                                                                cccAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGc                                                                                   >  8:177253/1‑149 (MQ=255)
                                                                cccAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGc                                                                                   >  5:136916/1‑149 (MQ=255)
                                                                cacAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGACGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCCACGGCGCCAGGGCGCAGTTCTCGAGGAACGGATTCATGAGGAGGATCGCGTTGTGGCCCCAGTAGTGc                                                                                   >  3:285878/3‑149 (MQ=255)
                                                                  cAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGa                                                                                 <  5:24893/149‑1 (MQ=255)
                                                                  cAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGAGGATCGCGTTGTGGCCCCAGTAGTGCGa                                                                                 <  1:18803/149‑1 (MQ=255)
                                                                     ccGGCACGCCGCAGCCGCGCAGCTTCGACGAAGCCGTGGGACAGGACGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGCCGCAGCGCTCGAGGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTc                                                                              <  8:370886/149‑1 (MQ=255)
                                                                       ggCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCAc                                                                            <  7:177253/149‑1 (MQ=255)
                                                                               gcATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGc                                                                    <  2:36578/149‑1 (MQ=255)
                                                                                 aTCAGCGCAGCTTCGTAGTAGTCGTGGGAAAGGTTGTCCCCGGCGAACGAGCCTTTGCCCGGCTACGGCGCCAGGGCGCAGTGCTCGATGATCGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCTACTTCCa                                                                  <  4:285878/149‑1 (MQ=255)
                                                                                    aGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGaa                                                               <  1:394257/149‑1 (MQ=255)
                                                                                    aGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGaa                                                               <  5:366715/149‑1 (MQ=255)
                                                                                    aGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCATTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCCGTAGTGCGACTACCCCAACTGCCAg                                                                 >  1:92974/1‑147 (MQ=255)
                                                                                     gcgcAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGTAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCa                                                                                          <  1:289420/121‑1 (MQ=255)
                                                                                     gcgcAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGTAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCa                                                                                          >  2:289420/1‑121 (MQ=255)
                                                                                         aGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCa                                                          <  5:123623/149‑1 (MQ=255)
                                                                                                      tagTGGGAAAGGATGGCACCGGCGAACGCGCCTTT‑CCCGGCAACGGCGCCATGGCGCAGTGCTCGATGAACGGCTTCATGCCGATGATCGAGTTGTGGCCCCAGTAGTGCGACTCACCACACTGCCAAAAGTTCAGACCGGCAGTGaa                                             >  4:296382/3‑148 (MQ=255)
                                                                                                      tCGTGGGAAAGGATGGCACCGGCGAACGCGCC‑TTGCCCGGCAACGGCGCCAGGGCCCAGTGCTCGATGAACGGCTTCATGCGGATGATAGCGTTGTTGCCCCAGTATTGCGACTCACCCAACTGACAGAAGTGCAGACCGGCACTGaa                                             >  4:291133/1‑148 (MQ=255)
                                                                                                           ggAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTg                                        <  3:380896/149‑1 (MQ=255)
                                                                                                            gAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTgg                                       >  6:82089/1‑149 (MQ=255)
                                                                                                             aaaGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTggg                                      <  2:295756/149‑1 (MQ=255)
                                                                                                              aaGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGc                                     >  2:227145/1‑149 (MQ=255)
                                                                                                                ggATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCg                                   >  3:14199/1‑149 (MQ=255)
                                                                                                                 gATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCg                                   >  1:56351/1‑148 (MQ=255)
                                                                                                                             gAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAg                      >  3:192189/1‑149 (MQ=255)
                                                                                                                             gAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAg                      >  3:186030/1‑149 (MQ=255)
                                                                                                                               tggcgcCGTTGCCGGGCAAAGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACGCGCGGGTGGc                     >  8:248181/3‑148 (MQ=255)
                                                                                                                                   gcCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAAct                >  5:290010/1‑149 (MQ=255)
                                                                                                                                      gttGCCGGGCAAAGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCCTACACGCGGGTAGCGAActgct             <  7:248181/148‑1 (MQ=255)
                                                                                                                                       ttGCCCGGCAACGGCGCCAGGGGGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAAGTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAActgct             <  4:78590/148‑1 (MQ=255)
                                                                                                                                        tGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAActgct             <  7:144459/147‑1 (MQ=255)
                                                                                                                                            cGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCATACACGCGGGTAGCGAACTGCTGCATGc       >  2:66828/1‑149 (MQ=255)
                                                                                                                                              gCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGt     <  6:136916/149‑1 (MQ=255)
                                                                                                                                                 aCGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGTTGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCa  <  5:338642/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCGTCAAGCAGGTTTGGCGGCAGTTCTTCGTAGCTGCCCGGCAAGTCGTAGGCAATCCACACACCCCAGCCGGCACGGCGCATCAGCGCAGCTTCGACGAAGTCGTGGGAAAGGATGGCACCGGCGAACGCGCCTTTGCCCGGCAACGGCGCCAGGGCGCAGTGCTCGATGAACGGCTTCATGCGGATGATCGCGTTGTGGCCCCAGTAGTGCGACTCACCCAACTGCCAGAAGTGCAGACCGGCAGTGAACAGTGGGCCGTACACGCGGGTAGCGAACTGCTGCATGCGTGCA  >  NC_002947/5725512‑5725805

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