Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 779,919 A→C 6.4% G184G (GGT→GGG)  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,9190AC 6.4% 191.9 / 4.7 63G184G (GGT→GGG) PP_0670bile acid/Na+ symporter family transporter
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base C (2/2);  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 = 6.16e-01

CCGGAATCACCGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTG  >  NC_002947/779771‑780059
                                                                                                                                                     |                                                                                                                                             
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                                                        <  2:330954/137‑1 (MQ=255)
                                                                                               ctggaaCGACGCCTGGGCGTGGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCAc                                                        >  1:330954/5‑141 (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)
                                                                                                                                                     |                                                                                                                                             
CCGGAATCACCGGCAGCAAACGGTCGAAATCACGCAGTACCACGTCAATCTGCGAAGCGATCACCACCATCAGCACCAAGGCCATGGCC‑GGCACCGGCAACCACGCCCAGGCGTCGTTCCAGGTCGCGACTGCGCGCGAGCGCCGGGCACCGGCGCTGGTGAGCACCGCGAGGGCCATCGGCAACACGATCAACAGCACGAAGGC‑TTCCAGGAACGGGGTGATGGAGATTGCCACCTGGCTGCTGCCACCCAGCATCAGGGCCAGGTACACCGGCAGCAGCACCAGTTG  >  NC_002947/779771‑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: ‑