Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 288,016 T→C 5.5% Y31C (TAC→TGC)  tauA ← taurine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947288,0160TC 5.5% 214.2 / 2.1 72Y31C (TAC→TGC) tauAtaurine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (34/34);  new base C (2/2);  total (36/36)
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.85e-01

CACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAATGATCT  >  NC_002947/287883‑288162
                                                                                                                                        |                                                                                                                                                  
cACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGCACAGTGGTCTGGTAGGCCACTGTCAgg                                                                                                                                      <  1:182975/149‑1 (MQ=255)
cACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGCCTGGTAGGCCACt                                                                                                                                            <  6:8094/142‑1 (MQ=255)
cACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGCCTACCAGACCACTGt                                                                                                                                          <  5:272524/144‑1 (MQ=255)
 aCGTCGCCGCTGGCCACGGCAGTGTTCACGTGTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTTCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGGACAGTGGTCTGGTAGGCCACTGTCAgg                                                                                                                                      <  1:188465/148‑1 (MQ=255)
 aCGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTc                                                                                                                                    >  1:199723/1‑149 (MQ=255)
  cGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCg                                                                                                                                   <  2:214422/149‑1 (MQ=255)
          cTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTg                                                                                                                           >  2:21016/1‑149 (MQ=255)
               cACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTgaccct                                                                                                                      >  8:275539/1‑144 (MQ=255)
                aCGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGATCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACATGCGCAAGTTTGGCCGGGTCGACATTGTTCTGGTAGGCCACTGTCAGGTCGGATGCATGGGcc                                                                                                                       >  1:64878/1‑148 (MQ=255)
                aCGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCaa                                                                                                                     >  3:10216/1‑149 (MQ=255)
                aCGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGCCTACCAGACCACTGTc                                                                                                                                         >  8:152736/1‑129 (MQ=255)
                aCGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGCCTACCAGACCACTGTc                                                                                                                                         <  7:152736/129‑1 (MQ=255)
                      gTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCt                                                                                                               <  4:162605/149‑1 (MQ=255)
                          tCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGt                                                                                                            <  2:327523/148‑1 (MQ=255)
                           cACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGt                                                                                                            <  3:118074/147‑1 (MQ=255)
                           cACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGt                                                                                                            >  4:118074/1‑147 (MQ=255)
                                 tGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGCCCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGc                                                                                                      >  6:176550/1‑147 (MQ=255)
                                 tGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGc                                                                                                      >  5:119876/1‑147 (MQ=255)
                                   cgcCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAgg                                                                                                  >  7:184375/1‑149 (MQ=255)
                                     ccGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGt                                                                                                >  5:20720/1‑149 (MQ=255)
                                        ttATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAg                                                                                             >  3:186540/1‑149 (MQ=255)
                                          aTCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAgg                                                                                            <  2:64878/148‑1 (MQ=255)
                                            cAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGggcg                                                                                         >  4:259335/1‑149 (MQ=255)
                                             aaaTTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGggcgg                                                                                        >  8:291067/1‑149 (MQ=255)
                                              aaTTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGggcggc                                                                                       <  4:300666/149‑1 (MQ=255)
                                                    cgcCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAg                                                                                 <  6:20720/149‑1 (MQ=255)
                                                     gcCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAgg                                                                                >  2:227253/1‑149 (MQ=255)
                                                      ccAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGCGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGgc                                                                               <  4:186540/149‑1 (MQ=255)
                                                       cAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGgcg                                                                              <  5:123026/149‑1 (MQ=255)
                                                             gATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGc                                                                          <  4:2545/148‑1 (MQ=255)
                                                             gATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCg                                                                         <  3:158075/149‑1 (MQ=255)
                                                             gATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCg                                                                         >  1:177128/1‑149 (MQ=255)
                                                                    gCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTg                                                                  >  5:33555/1‑149 (MQ=255)
                                                                    gCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTg                                                                  >  8:198912/1‑149 (MQ=255)
                                                                     ccTTGCTGGCTTTTTCGTAGGCGCCATCGACCTGCGCCCCTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGGCGGCTGCCTGGGCCAACCAGCGGGGGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGc                                                                 <  5:191143/149‑1 (MQ=255)
                                                                      cTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCg                                                                <  3:259335/149‑1 (MQ=255)
                                                                          cTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGAGTCGACACTGGTCTGGTAGCGCACTGTCAGGTCGGCTGCATGGGCCAACCAGCTGGTGCTGGCGAGCGTCTGGGGGGCGAACAGGCGCAGAGGAGCGGGCGGGAt                                                            >  8:112931/1‑149 (MQ=255)
                                                                           tGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGAt                                                            >  2:48329/1‑148 (MQ=255)
                                                                             gCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATc                                                           <  6:33555/147‑1 (MQ=255)
                                                                              cTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATg                                                        <  8:68396/149‑1 (MQ=255)
                                                                               tttttCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATgg                                                       <  8:184375/149‑1 (MQ=255)
                                                                               tttttCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCAAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATgg                                                       >  5:257292/1‑149 (MQ=255)
                                                                                     gTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGaa                                                  <  3:320198/148‑1 (MQ=255)
                                                                                     gTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAAc                                                 >  4:173072/1‑149 (MQ=255)
                                                                                     gTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAAc                                                 >  2:79684/1‑149 (MQ=255)
                                                                                      tAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAAcc                                                >  3:184007/1‑149 (MQ=255)
                                                                                           gCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCGGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCg                                           >  4:3287/1‑149 (MQ=255)
                                                                                             cATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGt                                          >  5:328462/1‑148 (MQ=255)
                                                                                                 gACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGc                                     <  1:21016/149‑1 (MQ=255)
                                                                                                       cgcCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCgg                                                                        >  4:103381/1‑108 (MQ=255)
                                                                                                       cgcCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCgg                                                                        <  3:103381/108‑1 (MQ=255)
                                                                                                         ccACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCgtg                             <  7:198912/149‑1 (MQ=255)
                                                                                                          cACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCgtgt                            >  6:270550/1‑149 (MQ=255)
                                                                                                          cACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCgtgt                            <  2:177128/149‑1 (MQ=255)
                                                                                                           aCTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGtt                           >  7:127570/1‑149 (MQ=255)
                                                                                                            cTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTg                          <  1:79684/149‑1 (MQ=255)
                                                                                                            cTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTg                          <  7:291067/149‑1 (MQ=255)
                                                                                                             tttGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTgg                         <  4:10216/149‑1 (MQ=255)
                                                                                                                   cGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACg                   <  4:232373/149‑1 (MQ=255)
                                                                                                                    gggTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACga                  <  6:119876/149‑1 (MQ=255)
                                                                                                                     ggTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTc                                <  2:185685/134‑1 (MQ=255)
                                                                                                                     ggTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTc                                >  1:185685/1‑134 (MQ=255)
                                                                                                                        cGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAg              >  7:184104/1‑149 (MQ=255)
                                                                                                                         tgACAGTGGCCTACCAGACCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAg              >  3:318396/2‑149 (MQ=255)
                                                                                                                            aCAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGc            >  3:202559/1‑149 (MQ=255)
                                                                                                                             cAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCt           <  1:48329/149‑1 (MQ=255)
                                                                                                                             cAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCt           <  1:154162/149‑1 (MQ=255)
                                                                                                                                     tGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGATCCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAATGATc   >  5:156515/1‑149 (MQ=255)
                                                                                                                                       gTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAATGAt    <  1:227253/146‑1 (MQ=255)
                                                                                                                                       gTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAATGAt    <  5:270550/146‑1 (MQ=255)
                                                                                                                                       gTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAATGATCt  >  5:135928/1‑148 (MQ=255)
                                                                                                                                       gTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAATGATCt  >  1:113673/1‑148 (MQ=255)
                                                                                                                                        |                                                                                                                                                  
CACGTCGCCGCTGGCCACGGCAGTGATCACGTCTGCGCCGTTATCAAATTTGCGCCAGTCGATGCTGGCCTTGCTGGCTTTTTCGTAGTCGCCATCGACCTGCGCCACTTTGGCCGGGTCGACAGTGGTCTGGTAGGCCACTGTCAGGTCGGCTGCCTGGGCCAACCAGCTGGTGCTGGCGAGGGTCAGGGCGGCGAACAGGCGCAGCGGAGCGTGCGGGATCATGGTTGAACCTCTCGTCAGGCAAGTCGTGTTGGATAACGAGAAGGCTAAATGATCT  >  NC_002947/287883‑288162

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: