Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,048,569 G→C 5.4% A344G (GCC→GGC)  PP_0907 ← RND transporter membrane fusion protein‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,048,5690GC 5.4% 240.0 / 4.3 74A344G (GCC→GGC) PP_0907RND transporter membrane fusion protein‑like protein
Reads supporting (aligned to +/- strand):  ref base G (33/37);  new base C (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

AGGGCGTTGTCCTGGCCTATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCGGCTGCGGGTTGGGCAGCA  >  NC_002947/1048422‑1048716
                                                                                                                                                     |                                                                                                                                                   
aGGGCGTTGTCCTGGCCTATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGc                                                                                                                                                    <  7:291443/149‑1 (MQ=255)
  ggCGTTGTCCTGGCCTATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTTGGCGa                                                                                                                                                  <  8:115288/149‑1 (MQ=255)
    cGTTGTCCTGGCCTATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGAc                                                                                                                                                 >  7:348879/1‑148 (MQ=255)
      ttGTCCTGGCCTATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGTCAGCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCACTCTTTACCGTGCGCTGTTCGGCGAc                                                                                                                                                 >  1:314658/1‑148 (MQ=255)
           cTGGCCTATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTggg                                                                                                                                         >  7:224397/1‑149 (MQ=255)
           cTGGCCTATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTggg                                                                                                                                         >  2:114611/1‑149 (MQ=255)
           cTGGCCTATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTC‑TTACCGTGCGCTGTTCGGCGACCAGCTggg                                                                                                                                         <  2:266933/148‑1 (MQ=255)
              gCCTATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGtt                                                                                                                                      >  5:319131/1‑149 (MQ=255)
               ccTATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTg                                                                                                                                     <  4:65313/149‑1 (MQ=255)
                cTATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGCGGTAgg                                                                                                                                    >  5:356586/1‑149 (MQ=255)
                 tATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGc                                                                                                                                   >  4:382910/1‑149 (MQ=255)
                   tGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGccc                                                                                                                                 <  5:24285/149‑1 (MQ=255)
                   tGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGccc                                                                                                                                 <  3:268842/149‑1 (MQ=255)
                   tGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGcc                                                                                                                                  <  6:319131/148‑1 (MQ=255)
                      gTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTg                                                                                                                              <  8:224397/149‑1 (MQ=255)
                         ttGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTgtc                                                                                                                            <  8:235695/148‑1 (MQ=255)
                               cACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTa                                                                                                                      >  6:14837/1‑148 (MQ=255)
                                 ccGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTacac                                                                                                                   <  4:373567/149‑1 (MQ=255)
                                   gggTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCg                                                                                                                 >  2:22927/1‑149 (MQ=255)
                                            tttCAGCTGCCCGGCGGTGCCCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTgtc                                                                                                                            >  5:325566/1‑129 (MQ=255)
                                            tttCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTgtc                                                                                                                            >  5:320193/1‑129 (MQ=255)
                                            tttCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTgtc                                                                                                                            <  6:325566/129‑1 (MQ=255)
                                            tttCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGGCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTgtc                                                                                                                            <  6:320193/129‑1 (MQ=255)
                                              tCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATcc                                                                                                      >  7:271441/1‑149 (MQ=255)
                                               cAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATcc                                                                                                      <  8:156140/148‑1 (MQ=255)
                                               cAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCt                                                                                                     >  1:258282/1‑149 (MQ=255)
                                                aGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCtt                                                                                                    <  5:39424/149‑1 (MQ=255)
                                                  cTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGt                                                                                                  <  1:341908/149‑1 (MQ=255)
                                                    gCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTg                                                                                                            <  7:14041/137‑1 (MQ=255)
                                                    gCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTg                                                                                                            >  8:14041/1‑137 (MQ=255)
                                                     cccGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTc                                                                                                 <  1:135040/147‑1 (MQ=255)
                                                           ggTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCg                                                                                         <  5:14837/149‑1 (MQ=255)
                                                             tGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGt                                                                                       <  3:158053/149‑1 (MQ=255)
                                                              gACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGTGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGCCCAGCTGGGGTTGGTCCTTGTCGTCGTACACCGGTTTGCCCTCCTTGTCCTTTCTCGGt                                                                                       >  6:188398/1‑148 (MQ=255)
                                                               accaccTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTg                                                                                      <  8:271441/148‑1 (MQ=255)
                                                                 caccTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAgg                                                                                   >  4:353308/1‑149 (MQ=255)
                                                                       t‑‑tCCCCGGCCTGCAGGCCTTCTCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACa                                                                              <  2:133044/148‑1 (MQ=255)
                                                                            cccGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATc                                                                                                       >  6:369037/1‑120 (MQ=255)
                                                                            cccGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATc                                                                                                       <  5:369037/120‑1 (MQ=255)
                                                                                   tGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCtt                                                                                                    <  4:271266/116‑1 (MQ=255)
                                                                                   tGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCtt                                                                                                    >  3:271266/1‑116 (MQ=255)
                                                                                     cAGGCCTTCGCTGACTACAACCACGCCGTCTCCGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGAACAGCTGGGGTTGGCACTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGtt                                                                   <  5:275068/147‑1 (MQ=255)
                                                                                       ggCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTTCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCa                                                               <  2:314658/149‑1 (MQ=255)
                                                                                                     aCAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAgg                                                  >  6:162962/1‑148 (MQ=255)
                                                                                                      cAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGtt                                                                   >  3:17969/1‑130 (MQ=255)
                                                                                                      cAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGtt                                                                   <  4:17969/130‑1 (MQ=255)
                                                                                                      cAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAgg                                                  <  8:111381/147‑1 (MQ=255)
                                                                                                      cAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAgg                                                  <  1:114611/147‑1 (MQ=255)
                                                                                                      cAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGt                                                 >  8:287199/1‑148 (MQ=255)
                                                                                                       aaCCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTaa                                               >  2:78955/1‑149 (MQ=255)
                                                                                                       aaCCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCAATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTaa                                               >  5:35010/1‑149 (MQ=255)
                                                                                                         ccACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATg                                             >  4:358586/1‑149 (MQ=255)
                                                                                                          cACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATgg                                            <  5:151723/149‑1 (MQ=255)
                                                                                                                   tcGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCg                                   >  5:21608/1‑149 (MQ=255)
                                                                                                                     ggCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATTCAGGGTGTAGGTAATGGCGCTTTCCGGc                                 <  1:22927/149‑1 (MQ=255)
                                                                                                                     ggCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTcc                                    <  3:382910/146‑1 (MQ=255)
                                                                                                                     ggCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCAAAAATTGTACTTTATAGGTGAGGCGAAATAAACCGAGTTGACAAAAAGGGTGTAGGTAATGGCGCTTTACGGc                                 >  8:248635/1‑149 (MQ=255)
                                                                                                                        gcTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGcac                               >  4:106023/1‑148 (MQ=255)
                                                                                                                                 tCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGc                     >  7:68750/1‑149 (MQ=255)
                                                                                                                                   tttACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGCGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCgg                   >  7:176416/1‑149 (MQ=255)
                                                                                                                                     tACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCgg                   <  7:248635/147‑1 (MQ=255)
                                                                                                                                       ccGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCGGCTGc               >  7:12491/1‑149 (MQ=255)
                                                                                                                                         gTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCGGCTGCgg             <  3:353308/149‑1 (MQ=255)
                                                                                                                                         gTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCGGCTGCgg             >  4:193961/1‑149 (MQ=255)
                                                                                                                                         cccagcTGGTCGCCGAACAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAg                                                   <  3:318969/107‑1 (MQ=255)
                                                                                                                                         cccagcTGGTCGCCGAACAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAg                                                   >  4:318969/5‑111 (MQ=255)
                                                                                                                                          tGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCGGCTGCggg            <  7:287199/149‑1 (MQ=255)
                                                                                                                                           cagcTGGTCGCCGAACAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCGGCTGCGGGt           >  5:330674/3‑149 (MQ=255)
                                                                                                                                             gcTGGTCGCCGAACAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCGGCTGCGGGTTg         <  6:330674/149‑1 (MQ=255)
                                                                                                                                              cTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACa                                                           <  2:87668/98‑1 (MQ=255)
                                                                                                                                              cTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACa                                                           >  1:87668/1‑98 (MQ=255)
                                                                                                                                                  tCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCGGCTGCGGGTTGGgcag    >  3:66839/1‑149 (MQ=255)
                                                                                                                                                   cGGCGACCAGCTTGGGTTGGCCCTTGTCTTCCTACACCGCTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGATTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCGGCTGCGGGTTGGgcag    <  6:356586/148‑1 (MQ=255)
                                                                                                                                                   cGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCGGCTGCGGGTTGGgcag    <  3:193961/148‑1 (MQ=255)
                                                                                                                                                    ggCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCGGCTGCGGGTTGGgcagca  <  6:35010/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGGGCGTTGTCCTGGCCTATGCGTATTGCGGCACCGGGTGTGAGTTTCAGCTGCCCGGCGGTGACCACCTGT‑‑TCCCCGGCCTGCAGGCCTTCGCTGACTACAACCACGCCGTCTCGGCGCTCGCCAGTCTTTACCGTGCGCTGTTCGGCGACCAGCTGGGGTTGGCCCTTGTCGTCGTACACCGGTTTGCCATCCTTGTCCTTTCTCGGTGAGGCGACATACACCGAGTTGCCATACAGGGTGTAGGTAATGGCGCTTTCCGGCACCACCACTTGCGGCTGCGGGTTGGGCAGCA  >  NC_002947/1048422‑1048716

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