Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,396 T→G 7.0% intergenic (+202/+27) yhdZ → / ← rrfF putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,3960TG 7.0% 260.1 / 6.0 113intergenic (+202/+27)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base T (39/63);  new base G (3/6);  total (43/70)
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.79e-01

Only 100 of 117 total aligned reads displayed.

GATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGC  >  NC_000913/3423249‑3423536
                                                                                                                                                   |                                                                                                                                            
gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACGCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCtt                                                                                                                                             >  1:1127893/1‑149 (MQ=255)
gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGATTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCtt                                                                                                                                             >  2:1012912/1‑149 (MQ=255)
 atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                                                                                                            <  1:935312/149‑1 (MQ=39)
 atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                                                                                                            >  2:502493/1‑149 (MQ=39)
  tttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTcc                                                                                                                                           <  2:635107/149‑1 (MQ=39)
          cTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                                                                                                                   <  2:845257/149‑1 (MQ=37)
           tGGCAGTTTATGGCGGGCGTCCTTCCCGCATCCCTCCGGTCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGATTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                                  <  2:268613/149‑1 (MQ=11)
            ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                 <  1:385232/149‑1 (MQ=37)
             gCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                                <  1:593676/149‑1 (MQ=37)
             gCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                                >  1:514899/1‑149 (MQ=37)
                 tttATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGt                                                                                                                            <  2:360368/149‑1 (MQ=14)
                 tggctGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTGCAACTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGt                                                                                                                            <  1:653952/145‑1 (MQ=14)
                  ttATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                                                                                                           >  2:1119435/1‑149 (MQ=39)
                    ttggctggcGTCCTGCCCGCCACCATCCGGGCCGTTGCTTCGCAACGTTCAAATCAGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                                                                                                                         <  2:1094494/148‑1 (MQ=255)
                    aTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                                                                                                                         <  1:300927/149‑1 (MQ=39)
                     tGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                                                                                                                        <  2:527687/149‑1 (MQ=39)
                       gcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCTCAGTCGGAAGACTGGGCCCTTTGTTTTAtc                                                                                                                      <  2:950685/149‑2 (MQ=14)
                       gcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt                                                                                                                      <  1:825069/149‑1 (MQ=39)
                       gcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt                                                                                                                      >  1:1024565/1‑149 (MQ=39)
                         gggcgTCCTGCCGGCGACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                                                                                                                    <  2:788020/149‑1 (MQ=255)
                         gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                                                                                                                    <  1:706925/149‑1 (MQ=39)
                          ggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAAATAAAACAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATCTgt                                                                                                                   <  1:280077/149‑2 (MQ=14)
                           gcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCTGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                                                                                                                  <  2:1106725/149‑1 (MQ=37)
                             gTCCTGCCCGCCACCATCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAGTTGTCCTACTCGGGAGAGTGTTCACCGACAATCAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTGTATTTGATGt                                                                                                                >  1:183286/1‑149 (MQ=255)
                                cTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                                                                                                             >  2:136183/1‑149 (MQ=255)
                                cTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                                                                                                             <  2:114220/149‑1 (MQ=255)
                                 tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                                                                                                            <  2:638021/149‑1 (MQ=255)
                                 tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                                                                                                            <  2:1024565/149‑1 (MQ=255)
                                 tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAACGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                                                                                                            <  1:405751/149‑1 (MQ=255)
                                    ccGACACCCTCCGGGCCGTTTCTTCGCAACGTTCAAATCCGCTCCCGGAGGATTTGTCCTACTCGGCAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCACAGTATTCCGAATGAGCATTTTGTTTTTATTGATGTCTGGCat                                                                                                         >  2:762655/1‑148 (MQ=14)
                                     cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATCTGTTGTTTgtcggt                                                                                                        <  2:422858/149‑6 (MQ=14)
                                        cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCTGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCATCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc                                                                                                     <  1:1070176/149‑1 (MQ=255)
                                        cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc                                                                                                     <  2:81895/149‑1 (MQ=255)
                                        cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc                                                                                                     <  2:462205/149‑1 (MQ=255)
                                        cACCCTCCGCGCCGTTGCATCGCAACGTTCAAATCCGCTCCCCGCGGATTTGTCCTACTCGCGAGAGTGTTCACCGACAAACAACACATAATACACAAGGCCAAGTCGTCCGACTGCGGAGTATGTTTTATTGGAGATCGGGCACttac                                                                                                     >  2:1108825/1‑147 (MQ=2)
                                          cccTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt                                                                                                   <  2:953802/149‑1 (MQ=255)
                                           ccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGTTGTCTGGCAGTTCCCTa                                                                                                  >  2:909396/1‑149 (MQ=255)
                                            cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTGTTATTTGATGTCTGGCAGTTCCGTaa                                                                                                 >  1:403746/1‑148 (MQ=255)
                                               gtgaacactctcccgagtaggacaaatccgccgggagcggatttgTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct                                                                                              >  1:411269‑M2/46‑149 (MQ=255)
                                               cGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct                                                                                              <  2:190482/149‑1 (MQ=255)
                                                ggccagTTCCTTCTCAAAGTTCAACTGCGCTCCCGGCGGATATGTCCTACTCGGGAGCGTGTTCACCGACAAACAACATATAAAACAAAAGGCACAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACtctc                                                                                             <  2:96812/144‑1 (MQ=14)
                                                  gCCGTTGCTTCGCAACGTGAAAATCCGATCCCGGCTTATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGc                                                                                           <  2:856098/149‑1 (MQ=14)
                                                   ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGTCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCa                                                                                          <  2:994445/149‑1 (MQ=255)
                                                    cGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCAt                                                                                         >  1:1147904/1‑149 (MQ=255)
                                                      ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATgg                                                                                       <  1:161368/149‑1 (MQ=255)
                                                      tggCTTCGCAACGCTCAAATCCGCTCCCGGAGGATTTGGCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATgg                                                                                       <  1:1100848/147‑1 (MQ=255)
                                                       tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg                                                                                      <  1:499374/149‑1 (MQ=255)
                                                       tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg                                                                                      <  1:190691/149‑1 (MQ=255)
                                                        gCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATgggg                                                                                     >  1:757048/1‑149 (MQ=255)
                                                           tCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGgagc                                                                                  >  2:557473/1‑148 (MQ=255)
                                                              cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                               >  1:703432/1‑149 (MQ=255)
                                                               aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAACGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAcccc                                                                              <  1:642647/149‑1 (MQ=255)
                                                                  gTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcaca                                                                           >  2:53128/1‑149 (MQ=255)
                                                                    acaaatccgccgggagcggatttgTCCTACTCGGGAGAGTGTTAACCGCCAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGCGACCCCACACt                                                                         >  1:681723‑M2/25‑149 (MQ=255)
                                                                       AATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGACCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGATCCCTAATCTCGCATGGGGAGTCCCAACACTAcc                                                                      <  1:768220‑M1/149‑1 (MQ=255)
                                                                          CCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAATACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGATTTATTTGATGTCTGGCAGTTCCGTACTCTCGCATGGGGAGACACCACACTAACATc                                                                   <  1:213805‑M1/149‑1 (MQ=255)
                                                                           CGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCg                                                                  <  1:551644‑M1/149‑1 (MQ=255)
                                                                              tCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGCCCCCACACTACCATCGgcg                                                               >  1:289844/1‑149 (MQ=255)
                                                                                  ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAc                                                           >  2:167305/1‑149 (MQ=255)
                                                                                      gatttgTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTTTTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGCCCCCCCCCTACCAACGGCGCTACGGCg                                                       >  2:589755‑M2/7‑149 (MQ=255)
                                                                                       atttgTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGt                                                      <  2:468425‑M2/144‑1 (MQ=255)
                                                                                       atttgTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGt                                                      >  2:705113‑M2/6‑149 (MQ=255)
                                                                                         ttgTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACCACCATCGGCGCTACGGCGttt                                                    <  1:932968‑M2/146‑1 (MQ=255)
                                                                                           gTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCa                                                  <  1:232542‑M2/148‑1 (MQ=255)
                                                                                           gTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCa                                                  <  1:502493‑M2/148‑1 (MQ=255)
                                                                                            tCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTATTATTTGATGAATGGCATATCCCTACTCTCGCATGGTTCGACACAACACTACAATCGCAGCTACGGCGattaaa                                                 >  1:297056/1‑146 (MQ=14)
                                                                                              cTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACtt                                               >  1:1029978/1‑149 (MQ=255)
                                                                                              cTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACtt                                               <  1:681165/149‑1 (MQ=255)
                                                                                                aCTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCCCACTACCATCGGCGCTACGGCGTGTCACTTCt                                             >  1:473391/1‑149 (MQ=255)
                                                                                                aCTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCAGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                             <  1:970593/149‑1 (MQ=255)
                                                                                                  tCGGGAGAGTGTTCACCGACAAACCACAGATAAAACAAAAGTCCCAGTCTTCCGACTGAGCGTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTTCGGCGTTTCACTTCTGa                                           <  2:1065551/149‑1 (MQ=255)
                                                                                                    gggAGAGTGTTCACCTACAAACAACATATAAAACACAAGGCCCAGTCTTCCGACTTAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGt                                         <  2:597557/149‑1 (MQ=255)
                                                                                                     ggtgagTGTTCACCGACAAACAACTCATAAACCTAATGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCTTGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGtt                                        <  2:403746/146‑1 (MQ=14)
                                                                                                      gagaGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCTCTTCTGAGTTc                                       >  1:808796/1‑149 (MQ=255)
                                                                                                      gagaGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                       >  2:863188/1‑149 (MQ=255)
                                                                                                         agTGTTCAACGACAAACAACAGATAAAACAAAATGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc                                    <  2:213805/149‑1 (MQ=255)
                                                                                                             ttCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                                <  2:593674/149‑1 (MQ=39)
                                                                                                              tCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATggg                               >  1:376074/1‑149 (MQ=39)
                                                                                                                  cGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCa                           >  1:44534/1‑149 (MQ=39)
                                                                                                                   gACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg                          <  2:851857/149‑1 (MQ=39)
                                                                                                                     caaacaaCTGATAACACAAACGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                        <  2:703432/149‑1 (MQ=255)
                                                                                                                       aacaacAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTgg                      <  2:514899/149‑1 (MQ=39)
                                                                                                                       aacaacAGATAAAAACAAAGGGCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTgg                      <  1:641797/149‑1 (MQ=255)
                                                                                                                         caacaGATAAACCAAAAGGCCCAGTCTTCCGACTGAGCCTTTGTTTTATTTGATGTCTGGCAGTACCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac                   <  2:1154657/149‑1 (MQ=255)
                                                                                                                         caacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCGTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTgggc                    >  1:32871/1‑148 (MQ=255)
                                                                                                                          aacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac                   <  1:559284/149‑1 (MQ=39)
                                                                                                                              caTAAAACAAAAGGCCCCTGATTCCGACTGATCCTTTTGTTTTATATGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTCCCATCGGCGCTACTGCGTTTCACTTCTGAGTTCGTAATTTTGTCAGGTGGGaccacc               <  2:508081/148‑1 (MQ=14)
                                                                                                                               aTAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACGCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGCCCCCcg              >  2:876197/1‑149 (MQ=255)
                                                                                                                                tAAAACAAAAGGCCCAGTCTTCCGACTGAGACTTATGTTTTATATGATGACAGGCAGTTCCATACTCTCCCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcgc             >  2:570680/1‑149 (MQ=11)
                                                                                                                                    aCAAAAGGCCCAGTCTTCTTACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACAGTACCCTCGGCGCTACGGCGTTTCACTTCTGAGTGCGTCATGGGGTCAGGTGGGACCACCGCGCTa         <  2:231170/149‑1 (MQ=14)
                                                                                                                                    aCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTa         <  1:136183/149‑1 (MQ=39)
                                                                                                                                      aaaaGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtt                                                     >  1:697010/1‑103 (MQ=11)
                                                                                                                                      aaaaGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGtt                                                     <  2:697010/103‑1 (MQ=11)
                                                                                                                                      aaaaGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACt       >  2:113209/1‑149 (MQ=38)
                                                                                                                                        aaGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCAt                                                                    <  2:864682/86‑1 (MQ=255)
                                                                                                                                        aaGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCAt                                                                    >  1:864682/1‑86 (MQ=255)
                                                                                                                                         cggCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGATCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgcc    <  1:23516/148‑1 (MQ=37)
                                                                                                                                           gCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgccgc  <  1:53128/149‑1 (MQ=37)
                                                                                                                                                   ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                              >  1:433449/1‑97 (MQ=40)
                                                                                                                                                   ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                              <  2:433449/97‑1 (MQ=40)
                                                                                                                                                   |                                                                                                                                            
GATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGC  >  NC_000913/3423249‑3423536

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: