Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,658,592 C→T 20.2% 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,5920CT20.2% 228.6 / 44.1 108A175V (GCC→GTC) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (42/43);  new base T (14/8);  total (57/51)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

Only 100 of 111 total aligned reads displayed.

GCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCT  >  NC_000913/1658493‑1658690
                                                                                                   |                                                                                                  
gCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGACAATATGACCCGCTCTTCGCCATCAGc                                                                                                    >  1:1972555/1‑100 (MQ=255)
 ccTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGcc                                                                                                   >  1:3611643/1‑100 (MQ=255)
  cTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCt                                                                                                  >  1:1805920/1‑100 (MQ=255)
      agaaATCCGTCGAACAATATTTCAATGATGCTGTATATATTCAGTACTCTTCGGGTATTGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggc                                                                                              <  1:2408176/100‑1 (MQ=255)
      agaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggc                                                                                              >  1:3780002/1‑100 (MQ=255)
      agaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggc                                                                                              >  1:3516548/1‑100 (MQ=255)
       gaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcg                                                                                             >  1:3211341/1‑100 (MQ=255)
         aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGt                                                                                           >  1:3075360/1‑100 (MQ=255)
         aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGt                                                                                           <  1:2011706/100‑1 (MQ=255)
          aaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTc                                                                                          <  1:4057623/100‑1 (MQ=255)
                cGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGc                                                                                    >  1:2085062/1‑100 (MQ=255)
                 gAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGcc                                                                                   >  1:2974286/1‑100 (MQ=255)
                   aCAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtg                                                                                 <  1:61615/100‑1 (MQ=255)
                   aCAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtg                                                                                 >  1:3950788/1‑100 (MQ=255)
                    cAATATCGCAATGAGGCTGTATATATTCAGCACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtga                                                                                <  2:3604290/100‑1 (MQ=255)
                     aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgat                                                                               >  2:861339/1‑100 (MQ=255)
                     aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgat                                                                               <  1:3819293/100‑1 (MQ=255)
                        atGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGaa                                                                            >  1:1127728/1‑100 (MQ=255)
                         tggtagtGAGGCGGTATATATTCAGCACTCTTCGAGGATCGTTGGCGGCAATATGTCCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtgat                                                                               <  2:430921/90‑1 (MQ=255)
                         tGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAAc                                                                           <  2:3010857/100‑1 (MQ=255)
                         tGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAAc                                                                           <  1:2326135/100‑1 (MQ=255)
                          ggtaaTGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACt                                                                          <  2:784888/97‑1 (MQ=255)
                          ggCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACt                                                                          >  2:4018561/1‑100 (MQ=255)
                           gCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTg                                                                         <  1:241237/100‑1 (MQ=255)
                              aTGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTa                                                                      >  2:2187708/1‑100 (MQ=255)
                              aTGAAGCTGTATATATTCAGTACTCTTCGGGGATCGTTTGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTAAAACGCCTCATGAGCt                                                                          <  2:3613107/96‑1 (MQ=255)
                                gAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcg                                                                    <  1:787729/100‑1 (MQ=255)
                                gAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcg                                                                    <  1:2047628/100‑1 (MQ=255)
                                  ggCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTGGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcggc                                                                  >  1:1000598/1‑100 (MQ=255)
                                  ggCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcggc                                                                  <  1:1064950/100‑1 (MQ=255)
                                    cTGTATATATTCAGTACTCTTTGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcggcgg                                                                <  1:1931685/100‑1 (MQ=255)
                                    cTGTATATATTCAGTACTCTTCGGGGATCGTGGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcggcgg                                                                <  2:483881/100‑1 (MQ=255)
                                    cTGTATATATTCAGTACTCTTCGGGGATCGGTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcggcgg                                                                <  2:3117272/100‑1 (MQ=255)
                                     tGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGt                                                               >  2:2884436/1‑100 (MQ=255)
                                     tGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGt                                                               >  2:1136984/1‑100 (MQ=255)
                                     tGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGt                                                               >  2:3210407/1‑100 (MQ=255)
                                      gTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGtt                                                              >  1:24976/1‑100 (MQ=255)
                                       tatataTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATTACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                             >  2:1545749/1‑100 (MQ=255)
                                       tatataTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                             <  1:2635689/100‑1 (MQ=255)
                                         tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCAc                                                           >  1:542065/1‑100 (MQ=255)
                                         tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCAc                                                           >  1:731834/1‑100 (MQ=255)
                                          atatTCAGTACTCTTTGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACt                                                          <  2:337915/100‑1 (MQ=255)
                                          atatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACt                                                          <  1:375131/100‑1 (MQ=255)
                                          atatTCAGTACCCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACt                                                          <  2:251839/100‑1 (MQ=255)
                                            atTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCa                                                        >  2:2322675/1‑100 (MQ=255)
                                             ttCAGTACTCTTCGGGGATCGTTGGAGGCAATTTGACCCGTTCTTCGCCATCAAGCTCGGCGGTCAGAGGCCTGAGGAACTGCTACGGAGGTTCGATCaa                                                       >  2:1932220/1‑100 (MQ=255)
                                               cAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAAcc                                                     >  1:305390/1‑100 (MQ=255)
                                                aGTACTCTTCGGGGATCGTTCGCGGCAATATAACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCCGATGAACTGCTACGGCGGTTCACCCAACCa                                                    >  2:1934128/1‑100 (MQ=255)
                                                  tACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGt                                                  <  1:1588254/100‑1 (MQ=255)
                                                  tACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGt                                                  <  1:1587337/100‑1 (MQ=255)
                                                   aCTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTa                                                 <  1:1320765/100‑1 (MQ=255)
                                                   aCTCTTCGGGGATCGTTAGCGGAAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATTAACTGCTCCGGCAGTTCACTCCACCAGt                                                  >  2:2855667/1‑99 (MQ=255)
                                                    ctctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTAt                                                >  1:919921/1‑100 (MQ=255)
                                                      ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATgg                                              >  2:2505173/1‑100 (MQ=255)
                                                      ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATgg                                              >  2:3668058/1‑100 (MQ=255)
                                                      ctTCGGCGATCGTTGGCGGCCGTATGTCGCGCTCGTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAAGGCGGTTCACTCAACCAGTATg                                               <  2:3664560/99‑1 (MQ=255)
                                                       ttCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGc                                             <  1:3454465/100‑1 (MQ=255)
                                                          ggggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTcc                                          >  2:1540442/1‑100 (MQ=255)
                                                           gggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCt                                         <  1:2648835/100‑1 (MQ=255)
                                                            ggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTa                                        <  2:3156838/100‑1 (MQ=255)
                                                             gATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTAc                                       >  2:1480860/1‑100 (MQ=255)
                                                              aTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACa                                      <  2:3663612/100‑1 (MQ=255)
                                                               tCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAg                                     <  1:735178/100‑1 (MQ=255)
                                                               tCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAg                                     <  2:457247/100‑1 (MQ=255)
                                                               tCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAg                                     <  2:2138021/100‑1 (MQ=255)
                                                                cGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGc                                    <  2:1394209/100‑1 (MQ=255)
                                                                 gTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCa                                   >  2:288827/1‑100 (MQ=255)
                                                                 gTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCa                                   >  1:4123568/1‑100 (MQ=255)
                                                                 gTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCa                                   >  1:872174/1‑100 (MQ=255)
                                                                 gTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCa                                   >  2:1495776/1‑100 (MQ=255)
                                                                 gTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCa                                   >  1:617795/1‑100 (MQ=255)
                                                                  ttGGCGTCAATATGCCCCCCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCCACGGCGGTTCACTCAACCAGTATGGCTCCTACa                                      <  2:3799052/96‑1 (MQ=255)
                                                                   tGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACt                                 >  2:1151087/1‑100 (MQ=255)
                                                                    ggcggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTg                                >  2:3401617/1‑100 (MQ=255)
                                                                     gcggcAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgc                               >  1:2769121/1‑100 (MQ=255)
                                                                       ggcAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgcgc                             >  1:2058244/1‑100 (MQ=255)
                                                                         cAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCaa                           >  1:509912/1‑100 (MQ=255)
                                                                              tGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTc                      >  1:2964050/1‑100 (MQ=255)
                                                                               gACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc                     <  1:2525938/100‑1 (MQ=255)
                                                                                 cccGCTCTTCGCCATCAGCCTCGGCGGTCAAGCGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtg                   >  1:361889/1‑100 (MQ=255)
                                                                                 cccGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtg                   >  1:328848/1‑100 (MQ=255)
                                                                                  ccGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACTGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgt                  <  1:2856924/100‑1 (MQ=255)
                                                                                   cGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgtg                 >  1:497212/1‑100 (MQ=255)
                                                                                    gCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGc                <  2:508157/100‑1 (MQ=255)
                                                                                    gCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGc                <  2:1127721/100‑1 (MQ=255)
                                                                                         tCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGc           <  1:12585/100‑1 (MQ=255)
                                                                                         tCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGc           <  2:919915/100‑1 (MQ=255)
                                                                                           gCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGccc         >  1:1138687/1‑100 (MQ=255)
                                                                                           gCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGccc         >  1:1974250/1‑100 (MQ=255)
                                                                                            ccATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCt        >  2:1992688/1‑100 (MQ=255)
                                                                                            ccATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCt        >  2:3291665/1‑100 (MQ=255)
                                                                                            ccATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCt        <  1:2485677/100‑1 (MQ=255)
                                                                                              aTCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTac      <  1:2152670/100‑1 (MQ=255)
                                                                                               tCAGCCTCGGCGGTCAAACGCCTTATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCGAATTTCCTGTGCCATGCCCTaca     <  2:1656194/100‑1 (MQ=255)
                                                                                               tCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTaca     <  1:1083555/100‑1 (MQ=255)
                                                                                               tCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTac      >  2:1979664/1‑99 (MQ=255)
                                                                                                 aGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACAcc   <  1:1371257/100‑1 (MQ=255)
                                                                                                 aGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACAcc   >  2:751177/1‑100 (MQ=255)
                                                                                                  gCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCt  >  2:310224/1‑100 (MQ=255)
                                                                                                  accTCGGCGGTCAAACGCCCGATTAACTGCTACCGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACAcc   >  2:3745889/2‑99 (MQ=255)
                                                                                                   |                                                                                                  
GCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCT  >  NC_000913/1658493‑1658690

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: