Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 779,918 C→G 6.2% A185P (GCC→CCC) ‡ PP_0670 ← bile acid/Na+ symporter family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947779,9180CG 6.2% 209.4 / 2.8 64A185P (GCC→CCC) ‡PP_0670bile acid/Na+ symporter family transporter
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base G (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGAATCACCGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTG  >  NC_002947/779770‑780059
                                                                                                                                                     |                                                                                                                                              
aCCGGAATCACCGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGc                                                                                                                                                >  7:29341/1‑149 (MQ=255)
 ctggAATCACCGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCACCGGGCa                                                                                                                                               <  6:298065/147‑1 (MQ=255)
         aCCGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGgcgc                                                                                                                                        <  8:29341/148‑1 (MQ=255)
         aCCGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGgcgc                                                                                                                                        <  5:90120/148‑1 (MQ=255)
            ggCAGCAACCGGTCGAAAACCCGCACTACCACGTCAATCTGCGAAGCGATCACCACCATCAGGACCTAGGCCATGGCC‑GGCACCGGCAAGCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGTGCCCGGCGCtcgc                                                                                                                                    <  7:110123/149‑4 (MQ=255)
              cagcaAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGt                                                                                                                                    >  3:36138/1‑147 (MQ=255)
               agcaAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAg                                                                                                                                 <  3:290358/149‑1 (MQ=255)
                  aaaCGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCAc                                                                                                                              >  6:199838/1‑149 (MQ=255)
                     cGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACcgc                                                                                                                           <  5:367719/149‑1 (MQ=255)
                        tCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAg                                                                                                                        >  5:390693/1‑149 (MQ=255)
                        tCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAg                                                                                                                        <  4:285674/149‑1 (MQ=255)
                        tCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAg                                                                                                                        <  6:390693/149‑1 (MQ=255)
                          gAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGTCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAgg                                                                                                                       <  5:388640/148‑1 (MQ=255)
                              tCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCCGGGCACCCGCAACCACGCCCAGGCCTCGTTCCAGGTAGCGACTGCCCGCGAGCCCCGGGCACCGGCGCTTGTGAGCACCGCGACGGCCa                                                                                                                   >  1:117054/1‑149 (MQ=255)
                              tCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGt                                                                                                                                    >  4:19366/1‑131 (MQ=255)
                                aCGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATc                                                                                                                 <  4:36138/148‑1 (MQ=255)
                                      tACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAaca                                                                                                          <  5:216951/149‑1 (MQ=255)
                                         cACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGa                                                                                                                                  <  3:125491/122‑1 (MQ=255)
                                         cACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGa                                                                                                                                  >  4:125491/1‑122 (MQ=255)
                                         cACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCAGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGa                                                                                                       <  2:184664/149‑1 (MQ=255)
                                               attCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGTCAACGACGCCCAGGCGTCGTTCCAGGTCGCGTCTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAAc                                                                                                  <  2:20284/146‑1 (MQ=255)
                                                aTCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAg                                                                                                <  5:69975/149‑1 (MQ=255)
                                                    gCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCAc                                                                                             >  5:401985/1‑148 (MQ=255)
                                                                ccaccaTCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAg                                                                                >  2:84177/1‑148 (MQ=255)
                                                                   ccaTCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAAc                                                                            >  1:92030/1‑149 (MQ=255)
                                                                    caTCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACg                                                                           >  8:372761/1‑149 (MQ=255)
                                                                             cAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATg                                                                   >  7:221831/1‑148 (MQ=255)
                                                                                  ccATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGAtt                                                             <  7:372761/149‑1 (MQ=255)
                                                                                    aTGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCa                                                                                                                               >  6:285589/1‑82 (MQ=255)
                                                                                    aTGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCa                                                                                                                               <  5:285589/82‑1 (MQ=255)
                                                                                    aTGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGc                                                           >  1:250411/1‑149 (MQ=255)
                                                                                        cc‑ggcaccggcaACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCAcc                                                       >  7:222812/1‑149 (MQ=255)
                                                                                             caccggcaACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGgc                                                   <  8:221831/149‑1 (MQ=255)
                                                                                                ctggaaCGACGCCTGGGCGTGGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCAc                                                        >  1:330954/5‑141 (MQ=255)
                                                                                                ctggaaCGACGCCTGGGCGTGGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCAc                                                        <  2:330954/137‑1 (MQ=255)
                                                                                                cggcaACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTg                                                                      <  3:269885/127‑1 (MQ=255)
                                                                                                cggcaACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGgctgc                                                >  5:6397/1‑149 (MQ=255)
                                                                                                 ggcaACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGgctgct                                               >  5:99050/1‑149 (MQ=255)
                                                                                                  gcaACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGctgctg                                              <  2:92030/149‑1 (MQ=255)
                                                                                                    aaCCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGc                                                           >  1:380810/1‑134 (MQ=255)
                                                                                                    aaCCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGc                                                           <  2:380810/134‑1 (MQ=255)
                                                                                                      ccACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTGCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCAc                                          >  5:359994/1‑149 (MQ=255)
                                                                                                         cGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGTCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGCTTTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCAccc                                        <  1:84177/149‑1 (MQ=255)
                                                                                                          gCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAg                                      >  8:130887/1‑149 (MQ=255)
                                                                                                                    tgtTCCTGGTCGCGACTTCGCGCGATCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGTGCCATCGTCAAGACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGTTGGAGTTTGCCACCTGGCTGCTGCCACCCAGCATCAGGGc                             <  6:400516/147‑1 (MQ=255)
                                                                                                                      ttCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGc                             >  5:400516/1‑146 (MQ=255)
                                                                                                                        ccAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCa                                       <  5:113715/134‑1 (MQ=255)
                                                                                                                        ccAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCa                                       >  6:113715/1‑134 (MQ=255)
                                                                                                                         cAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGTCACCCAGCATCAGGGCCAGGTa                       <  2:394493/149‑1 (MQ=255)
                                                                                                                            gTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTacac                    >  8:77338/1‑149 (MQ=255)
                                                                                                                             tCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACAcc                   <  6:99050/149‑1 (MQ=255)
                                                                                                                               gcgACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCgg                 <  6:51522/149‑1 (MQ=255)
                                                                                                                               gcgACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCgg                 >  5:51522/1‑149 (MQ=255)
                                                                                                                                 gACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGgca               <  5:199838/149‑1 (MQ=255)
                                                                                                                                  aCTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGgcag              <  2:250411/149‑1 (MQ=255)
                                                                                                                                      cgcgcgAGCGCCGGGCACCGTCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGgaga                                                               <  7:136730/96‑1 (MQ=255)
                                                                                                                                      cgcgcgAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGgaga                                                               >  8:136730/1‑96 (MQ=255)
                                                                                                                                      cgcgcgAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGgcagcagc          <  6:359994/149‑1 (MQ=255)
                                                                                                                                       gcgcgAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGgc                                                   >  7:311408/1‑107 (MQ=255)
                                                                                                                                       gcgcgAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGgc                                                   <  8:311408/107‑1 (MQ=255)
                                                                                                                                         accaGCGCCGGTGCCCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACAGGCAGCAGCAc        >  3:336798/4‑148 (MQ=255)
                                                                                                                                           gAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCAACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAg     >  6:267558/1‑149 (MQ=255)
                                                                                                                                           caGCGCCGGTGCCCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCa                                       <  7:121305/114‑1 (MQ=255)
                                                                                                                                           caGCGCCGGTGCCCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCa                                       >  8:121305/2‑115 (MQ=255)
                                                                                                                                              cgcCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTg  >  7:220583/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
ACCGGAATCACCGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTG  >  NC_002947/779770‑780059

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: ‑