Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,658,592 C→T A175V (GCC→GTC)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,658,5920CT99.2% 382.7 / NA 119A175V (GCC→GTC) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (59/58);  total (59/60)
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

Only 100 of 119 total aligned reads displayed.

CTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACA  >  NC_000913/1658495‑1658687
                                                                                                 |                                                                                               
cTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCt                                                                                               >  1:152110/1‑100 (MQ=255)
  gaagaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcg                                                                                             >  1:3229478/1‑100 (MQ=255)
    agaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggc                                                                                           >  2:481681/1‑100 (MQ=255)
     gaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggcg                                                                                          >  1:511573/1‑100 (MQ=255)
     gaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggcg                                                                                          >  2:2301422/1‑100 (MQ=255)
      gatctcGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggcgg                                                                                         <  1:2710341/95‑1 (MQ=255)
       aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGt                                                                                        <  2:80450/100‑1 (MQ=255)
       aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGt                                                                                        <  2:1205701/100‑1 (MQ=255)
       aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGt                                                                                        >  2:3013568/1‑100 (MQ=255)
        aaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTc                                                                                       >  2:1723156/1‑100 (MQ=255)
         aCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCa                                                                                      <  2:3316183/100‑1 (MQ=255)
          ccGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCaa                                                                                     <  2:3201655/100‑1 (MQ=255)
          ccGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCaa                                                                                     <  2:305422/100‑1 (MQ=255)
          ccGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCa                                                                                      <  2:2625152/99‑1 (MQ=255)
          ccGTCGAACAATATGGCAATGACGCTCTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGt                                                                                        <  1:2044153/97‑1 (MQ=255)
           cGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggcggagatc                                                                                    >  2:2710384/1‑95 (MQ=255)
               gAACAATATGGCAATGAGGCTGAATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATCACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGc                                                                                 >  1:2026446/1‑99 (MQ=255)
                aaCAGTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCt                                                                               <  1:1721963/100‑1 (MQ=255)
                aaCAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCt                                                                               <  2:3139329/100‑1 (MQ=255)
                aaCAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCt                                                                               <  2:495366/100‑1 (MQ=255)
                 aCAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtg                                                                              >  1:202509/1‑100 (MQ=255)
                  cAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGTGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtga                                                                             <  2:2226885/100‑1 (MQ=255)
                  cAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtga                                                                             <  2:2141569/100‑1 (MQ=255)
                  cAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtga                                                                             <  2:1213168/100‑1 (MQ=255)
                  cAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtga                                                                             <  2:2141823/100‑1 (MQ=255)
                   aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtgat                                                                            <  2:2226714/100‑1 (MQ=255)
                   aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtgat                                                                            <  2:185726/100‑1 (MQ=255)
                     tatGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtgatga                                                                          >  1:1006021/1‑100 (MQ=255)
                      atGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGaa                                                                         >  2:3318367/1‑100 (MQ=255)
                        ggCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACt                                                                       <  2:1461497/100‑1 (MQ=255)
                          cAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGc                                                                     >  2:1637731/1‑100 (MQ=255)
                          cAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGACGGCAATATGACCCGCTCTACGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGc                                                                     >  1:1900728/1‑100 (MQ=255)
                             tGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAAccgccgt                                                                  >  2:3403599/1‑94 (MQ=255)
                                ggCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAcggc                                                               >  1:3025651/1‑100 (MQ=255)
                                   tGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGt                                                            >  2:259037/1‑100 (MQ=255)
                                   tGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGt                                                            <  2:1408278/100‑1 (MQ=255)
                                     tatataTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                          <  2:1295238/100‑1 (MQ=255)
                                     tatataTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                          <  2:2878795/100‑1 (MQ=255)
                                     tatataTTCAGTACTCGTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                          <  1:188160/100‑1 (MQ=255)
                                       tatatTCAGTACTCTTCGGGGATCGTTGGCGGTAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCAc                                                        >  2:1077247/1‑100 (MQ=255)
                                       tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCAc                                                        <  1:1065618/100‑1 (MQ=255)
                                        atatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACt                                                       >  1:1962374/1‑100 (MQ=255)
                                         gatTCTGGACTTTTTGGGTATCGTTGGGGGCAATATGACCCGCGCTTCGCCGTCAGTGTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACt                                                       <  1:1546687/98‑1 (MQ=255)
                                          atTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGAACCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTTCGGCGGTTCACTCa                                                     >  1:497222/1‑100 (MQ=255)
                                           ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCaa                                                    >  1:992611/1‑100 (MQ=255)
                                            tCAGTAGTCTTCGGAGATCGTTGGTGGCAATATGTCCCGCTCTGCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAAc                                                   <  1:475208/100‑1 (MQ=255)
                                            ctaGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAAc                                                   <  2:35567/98‑1 (MQ=255)
                                              aGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAAcag                                                 >  1:35567/1‑98 (MQ=255)
                                               gTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCa                                                 >  1:2367048/1‑99 (MQ=255)
                                                tcctctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGt                                               <  1:1085181/98‑1 (MQ=255)
                                                 aCTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTa                                              <  2:1322741/100‑1 (MQ=255)
                                                 aCTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTa                                              >  2:1389381/1‑100 (MQ=255)
                                                 aCTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTa                                              <  2:1416288/100‑1 (MQ=255)
                                                 aCTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTa                                              <  2:1616267/100‑1 (MQ=255)
                                                   tctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATg                                            <  2:2527861/100‑1 (MQ=255)
                                                    ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATgg                                           >  2:1326741/1‑100 (MQ=255)
                                                       cGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATgg                                           >  1:2952781/1‑97 (MQ=255)
                                                       cGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTc                                        <  1:1224501/100‑1 (MQ=255)
                                                         gggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCt                                      >  1:85997/1‑100 (MQ=255)
                                                         gggATCGTTGGCGGCAATAAGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTcc                                       <  1:381102/99‑1 (MQ=255)
                                                          ggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTa                                     >  1:2008734/1‑100 (MQ=255)
                                                           gATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTAc                                    <  2:3382048/100‑1 (MQ=255)
                                                           gATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTAc                                    <  2:152112/100‑1 (MQ=255)
                                                            aTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACa                                   <  2:2471318/100‑1 (MQ=255)
                                                             tCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAg                                  >  2:3010036/1‑100 (MQ=255)
                                                              cGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGc                                 >  2:939053/1‑100 (MQ=255)
                                                              cGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGc                                 <  2:16918/100‑1 (MQ=255)
                                                               gTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCa                                >  1:2688234/1‑100 (MQ=255)
                                                                ttGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCAc                               <  2:1329495/100‑1 (MQ=255)
                                                                 tGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACt                              >  1:686725/1‑100 (MQ=255)
                                                                  ggggacAATAAGACCCGCTCTTAGACATCAGGCTCGTCGGTCAAACACCTGATGATCTGCTACGGCGGTTCTCTCAACCAGTATGGCTCCTACAGc                                 <  1:3371635/91‑1 (MQ=255)
                                                                  ggcggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTg                             >  2:2166538/1‑100 (MQ=255)
                                                                   gcggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgc                            >  1:2837549/1‑100 (MQ=255)
                                                                    cggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgag                           >  1:2416783/1‑98 (MQ=255)
                                                                       cAATATGTCCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCaa                        <  1:2581837/100‑1 (MQ=255)
                                                                        aaTATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCaaa                       >  2:3437876/1‑100 (MQ=255)
                                                                        aaTATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCaaa                       >  1:1286068/1‑100 (MQ=255)
                                                                         atatGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAt                      <  2:1687935/100‑1 (MQ=255)
                                                                          tatGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAtt                     >  1:3293260/1‑100 (MQ=255)
                                                                          tatGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAtt                     >  1:1644927/1‑100 (MQ=255)
                                                                            tGACCCGCTCTTCGCCATCAGTCTCGGCGGTCACACGCCTGATGAACTCCTACGGCGGTTGACTCACCCGGGATGGCTCCCACAGCACTGAGCaa                        >  1:3001561/1‑95 (MQ=255)
                                                                            tGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTc                   >  2:647159/1‑100 (MQ=255)
                                                                             gACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCCACCAGTCTGGCTCCTACAGCACTGCGCACATTTcc                  >  1:610789/1‑100 (MQ=255)
                                                                             gACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCCGTATGGCTACTACAGCACTGCGCAAATTTcc                  >  1:2531918/1‑100 (MQ=255)
                                                                             gACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc                  >  2:3024894/1‑100 (MQ=255)
                                                                             gACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATAGCTCCTACAGCACTGCGCAAATTTcc                  >  1:3229633/1‑100 (MQ=255)
                                                                              aCCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCt                 >  2:421505/1‑100 (MQ=255)
                                                                               cccGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtg                >  1:260278/1‑100 (MQ=255)
                                                                               cccGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtg                <  2:2868429/100‑1 (MQ=255)
                                                                                  gCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGc             <  1:3270242/100‑1 (MQ=255)
                                                                                   ctctTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGcc            >  1:2991311/1‑100 (MQ=255)
                                                                                     ctTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCAt          >  2:193515/1‑100 (MQ=255)
                                                                                       tCGCCATTAGTATCGGCGGTCAAACGCCCGAAGAACTGATACGGCGGTTCACTCAACCTGTATGGCTCGTACAGCATTGCGCAAATTTCCTGTGCCATGc        <  1:1747590/100‑1 (MQ=255)
                                                                                         gCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGccc      >  2:1429915/1‑100 (MQ=255)
                                                                                         gCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGccc      <  1:672213/100‑1 (MQ=255)
                                                                                          ccATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCt     <  1:1077234/100‑1 (MQ=255)
                                                                                           cATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTa    >  1:1917281/1‑100 (MQ=255)
                                                                                           cATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTa    >  1:601979/1‑100 (MQ=255)
                                                                                            aTCAGTCTCGGCGCTCAAACGCCTGATGAACTGCTACGCCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTac   >  1:2144587/1‑100 (MQ=255)
                                                                                             tCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCTAATTTCCTGTGCTATGCCCTaca  <  1:525001/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACA  >  NC_000913/1658495‑1658687

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: