Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 287,067 G→C 6.4% P19R (CCG→CGG)  tauB‑II ← taurine ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947287,0670GC 6.4% 190.8 / 2.0 60P19R (CCG→CGG) tauB‑IItaurine ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base C (2/2);  total (27/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 = 1.00e+00

CGGCACACCATCGAGGGTGATGTGACCGCCACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGGCAGCACTGCTTCGACTTTGCCTTGCTGCTTGAGG  >  NC_002947/286919‑287214
                                                                                                                                                    |                                                                                                                                                   
cGGCAGACCATCGAGGGTGATGTGACCGCCAGTGGGGGCGACGAAGCCGGGGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGTCCCCAGGGCCACCAGCAGTTGGCGTTGTCCCAGGCTCATGTTGATATCGGCCAGCACCg                                                                                                                                                     <  7:140942/149‑1 (MQ=255)
    acacCATCGAGGGTGATGTGACCGCCACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGt                                                                                                                                                 >  5:52986/1‑149 (MQ=255)
      acCATCGAGGGTGATGTGACCGCCACTGGGGGCGACGAAGCCGGCTATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCGGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGc                                                                                                                                               >  1:84651/1‑149 (MQ=255)
         aTCGAGGGTGATGTGACCGCCACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCCCGGTGCTgg                                                                                                                                            <  4:5308/149‑1 (MQ=255)
          tCGAGGGTGATGTGACCGCCACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCCCGGTGCTGgc                                                                                                                                           >  7:59437/1‑149 (MQ=255)
           cGAGGGTGATGTTTCCGCCACTGGGGTCGACGATGACGTAGATGAGGTTGAGCAGTTATGTCTTGCGGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTTGTCCCAGGCTCAGGTTGATATCGGACATCACGCCGGTGCTgta                                                                                                                                           <  8:59437/148‑3 (MQ=255)
              gggTGATGTGACCGCCACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCt                                                                                                                                       >  4:141767/1‑149 (MQ=255)
               ggTGATGTGACCGCCACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTg                                                                                                                                      <  6:177344/149‑1 (MQ=255)
                gTGATGTGACCGCCACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTgg                                                                                                                                     <  3:66361/149‑1 (MQ=255)
                 tGATGTGACCGCCACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTggg                                                                                                                                    >  2:158634/1‑149 (MQ=255)
                         ccGCCACTGGGGGCGACGAAGCCGGCGATGTGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTg                                                                                                                                                >  1:22351/1‑129 (MQ=255)
                         ccGCCACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTg                                                                                                                                                <  2:22351/129‑1 (MQ=255)
                            ccACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCt                                                                                                                         >  3:162959/1‑149 (MQ=255)
                            ccACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCt                                                                                                                         >  8:183049/1‑149 (MQ=255)
                             cccTGGGGGCGACGACGCGGGCGTTGTGGTGGTGCAGTGCGGCCTGGCCGCTGCCGGACGCCCCCAGGGCCACCCGCAGTGGGCGTGGCCCCAGGCTCCGGTTTAAATCGGCCAGCCCCGGGGGGCTGGCGCCTGGGTACTGTGCGCTg                                                                                                                        <  2:264916/147‑1 (MQ=255)
                             cACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCt                                                                                                                         >  4:2450/1‑148 (MQ=255)
                                     gCGACGAAGCCGGCGATGAGGTGGCGCAGTGAGGGCTGGCCGCTGCCGGACGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGGGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTc                                                                                                                <  6:85251/149‑1 (MQ=255)
                                         cgaAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGc                                                                                                            <  6:52986/149‑1 (MQ=255)
                                            aGCCGGCTATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCg                                                                                                         <  2:84651/149‑1 (MQ=255)
                                             gCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGa                                                                                                        <  3:141767/149‑1 (MQ=255)
                                                  cGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAg                                                                                                   <  5:178231/149‑1 (MQ=255)
                                                   gATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGGACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAgg                                                                                                  <  3:28737/149‑1 (MQ=255)
                                                   gATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAgg                                                                                                  >  4:76995/1‑149 (MQ=255)
                                                         gTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCCCGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGa                                                                                            >  5:175871/1‑149 (MQ=255)
                                                            gAGCAGGTAGGTCTTGCCGCTGCGGGACGGCCCCTGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACGGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACg                                                                                          <  6:181829/148‑1 (MQ=255)
                                                             aGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACgg                                                                                         <  4:77066/148‑1 (MQ=255)
                                                               cAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCt                                                                                      >  6:13180/1‑149 (MQ=255)
                                                                 gTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTcc                                                                                    <  3:2450/149‑1 (MQ=255)
                                                                      gTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGAt                                                                                <  4:162959/148‑1 (MQ=255)
                                                                            ccGCTGCCGGAAGGCCCCAGGGCCACCAGCAGGTGGCGTGGGCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGACTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGttt                                                                          <  6:213560/148‑1 (MQ=255)
                                                                               cTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTggaa                                                                      <  8:81432/149‑1 (MQ=255)
                                                                                tGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCTAGGCCATGACGGGCTCCGGATCAGTTTggaag                                                                     <  7:183049/149‑1 (MQ=255)
                                                                                  ccGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGTGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTggaagg                                                                    >  6:100006/1‑148 (MQ=255)
                                                                                  ccGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTc                                                                                                          >  6:45604/1‑110 (MQ=255)
                                                                                  ccGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTc                                                                                                          <  5:45604/110‑1 (MQ=255)
                                                                                    ggAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGAt                                                                                                                      >  5:117693/1‑96 (MQ=255)
                                                                                    ggAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGAt                                                                                                                      <  6:117693/96‑1 (MQ=255)
                                                                                    ggAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGcc                                                                                                <  1:2378/118‑1 (MQ=255)
                                                                                    ggAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGcc                                                                                                >  2:2378/1‑118 (MQ=255)
                                                                                       aGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTGGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTg                                                              <  8:47123/149‑1 (MQ=255)
                                                                                          ccccAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTgg                                                             >  6:71138/1‑147 (MQ=255)
                                                                                          ccccAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTgg                                                             <  5:71138/147‑1 (MQ=255)
                                                                                                      cagcagTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGc                                               >  5:277223/1‑149 (MQ=255)
                                                                                                         cagTTGGCGTGGTGCCAGGCTCAGGTTGATATCGGCCAGCTGCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGc                                               <  4:77642/146‑1 (MQ=255)
                                                                                                         cagTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCACCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGCAACTTGGCGCTGACGTACGGc                                               >  3:77642/1‑146 (MQ=255)
                                                                                                           gTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGagta                                          <  3:127999/149‑1 (MQ=255)
                                                                                                            ttGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGagtag                                         >  3:121225/1‑149 (MQ=255)
                                                                                                              ggCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTc                                       >  3:194443/1‑149 (MQ=255)
                                                                                                                    gTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGgcag                                 <  3:61632/149‑1 (MQ=255)
                                                                                                                     tCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGgcag                                 >  7:185018/1‑148 (MQ=255)
                                                                                                                               aGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGGCAGCACTGCTTCGa                      >  5:124100/1‑149 (MQ=255)
                                                                                                                                 gTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGGCAGCACTGCTTCGACt                    >  6:277257/1‑149 (MQ=255)
                                                                                                                                    gATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGGCAGCACTGCTTCGACTTTg                 >  6:144492/1‑149 (MQ=255)
                                                                                                                                     atatCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGGCAGCACTGCTTCGACTTTg                 <  5:13180/148‑1 (MQ=255)
                                                                                                                                      tatCGGCCAGCACCGGGGTGCTGGTGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGGCAGCACTGCTTCGACTTTGcc               <  5:100006/149‑1 (MQ=255)
                                                                                                                                         cGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCgg                                    <  6:53766/125‑1 (MQ=255)
                                                                                                                                         cGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCgg                                    >  5:53766/1‑125 (MQ=255)
                                                                                                                                             cAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGGCAGCACTGCTTCGACTTTGCCTtgctgc        <  8:185018/149‑1 (MQ=255)
                                                                                                                                              aGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGGCAGCACTGCTTCGACTTTGCCTTgctgct       <  6:124100/149‑1 (MQ=255)
                                                                                                                                                 aCCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGGCAGCACTGCTTCGACTTTGCCTTGCTGCtt      >  4:125253/1‑147 (MQ=255)
                                                                                                                                                   cGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGGCAGCACTGCTTCGACTTTGCCTTGCTGCTTGAgg  <  4:121225/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGCACACCATCGAGGGTGATGTGACCGCCACTGGGGGCGACGAAGCCGGCGATGAGGTTGAGCAGTGAGGTCTTGCCGCTGCCGGAAGGCCCCAGGGCCACCAGCAGTTGGCGTGGTCCCAGGCTCAGGTTGATATCGGCCAGCACCGGGGTGCTGGCGCCTGGGTACTGTGCGCTGATGCGCTCCAGCTCGAGCAAGGCCATGACGGGCTCCGGATCAGTTTGGAAGGAACTTGGCGCTGACGTACGGCGAGTAGTCCGGCAGCACTGCTTCGACTTTGCCTTGCTGCTTGAGG  >  NC_002947/286919‑287214

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