Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 614,858 G→T 5.5% G234G (GGG→GGT)  entF → enterobactin synthase multienzyme complex component, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913614,8580GT 5.5% 248.6 / 3.8 73G234G (GGG→GGT) entFenterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (38/31);  new base T (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

GCGAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTGGGCTCGGCGGC  >  NC_000913/614710‑615004
                                                                                                                                                    |                                                                                                                                                  
gcgAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACggg                                                                                                                                                    <  3:583216/149‑1 (MQ=255)
  gAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGa                                                                                                                                                  <  2:241817/149‑1 (MQ=255)
   aaaGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGaa                                                                                                                                                 <  3:250430/149‑1 (MQ=255)
   aaaGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGaa                                                                                                                                                 <  1:899266/149‑1 (MQ=255)
      gCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTc                                                                                                                                              >  2:306231/1‑149 (MQ=255)
              tGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTCCCCGTCGGTAAATTccag                                                                                                                                           <  4:359510/144‑3 (MQ=255)
              tGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTCCCCGTCGGTAAATTccag                                                                                                                                           >  3:359510/1‑142 (MQ=255)
                gCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTgg                                                                                                                                    >  4:139481/1‑149 (MQ=255)
                  aGCGTGATGCGGCATTCTGGGCAGAACAGTGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCt                                                                                                                                  <  3:78412/149‑1 (MQ=255)
                       gATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCa                                                                                                                             >  2:884327/1‑149 (MQ=255)
                          gCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACt                                                                                                                          <  2:525387/149‑1 (MQ=255)
                            ggCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACttt                                                                                                                        >  3:746571/1‑149 (MQ=255)
                               aTTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAg                                                                                                                     <  1:97726/149‑1 (MQ=255)
                                ttCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAgg                                                                                                                    >  2:50992/1‑149 (MQ=255)
                                 tCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGgt                                                                                                                   <  3:353148/149‑1 (MQ=255)
                                    gggCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGtgtg                                                                                                                <  4:495432/149‑1 (MQ=255)
                                      gCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTCCCCGTCGGTAAATTcca                                                                                                                                            <  1:433076/119‑2 (MQ=255)
                                      gCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTCCCCGTCGGTAAATTcca                                                                                                                                            >  2:433076/1‑118 (MQ=255)
                                        aGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGc                                                                                                            >  4:97992/1‑149 (MQ=255)
                                        aGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGc                                                                                                            <  2:339043/149‑1 (MQ=255)
                                               cgtcgtCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGa                                                                                                     >  1:761559/1‑149 (MQ=255)
                                                gtcgtcAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGAt                                                                                                    >  4:399786/1‑149 (MQ=255)
                                                   gtcAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTa                                                                                                 >  4:11780/1‑149 (MQ=255)
                                                     cAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGc                                                                                               >  4:68008/1‑149 (MQ=255)
                                                     cAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGc                                                                                               <  3:97992/149‑1 (MQ=255)
                                                      aaCTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGcc                                                                                              >  2:692046/1‑149 (MQ=255)
                                                      aaCTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGcc                                                                                              >  3:703924/1‑149 (MQ=255)
                                                         tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                           <  2:289187/149‑1 (MQ=255)
                                                         tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                           <  4:363666/149‑1 (MQ=255)
                                                         tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                           >  4:429211/1‑149 (MQ=255)
                                                         tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                           >  3:311761/1‑149 (MQ=255)
                                                          gccgccgccCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTg                                                                                          <  2:536309/149‑1 (MQ=255)
                                                                gcccgcGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTgg                                                                                    <  3:68008/149‑1 (MQ=255)
                                                                gcccgcGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTgg                                                                                    <  3:545366/149‑1 (MQ=255)
                                                                       tCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCtt                                                                             >  3:359199/1‑149 (MQ=255)
                                                                           tttCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTgg                                                                         >  3:345224/1‑149 (MQ=255)
                                                                               tcCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTgg                                                                     >  3:32522/1‑149 (MQ=255)
                                                                               tcCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTgg                                                                     <  4:435379/149‑1 (MQ=255)
                                                                               tcCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTgg                                                                     <  4:432426/149‑1 (MQ=255)
                                                                                 cGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTgggg                                                                   <  2:890105/149‑1 (MQ=255)
                                                                                   gCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCg                                                                 >  1:816236/1‑149 (MQ=255)
                                                                                    cACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGa                                                                >  4:353690/1‑149 (MQ=255)
                                                                                     aCCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGAt                                                               <  3:454375/149‑1 (MQ=255)
                                                                                      ccTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGCTTTAGCCCTTGCGCTGGCAGCATTGTGCCTGCGGCGAtt                                                              <  2:761559/149‑1 (MQ=255)
                                                                                       cTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATtg                                                             >  4:529367/1‑149 (MQ=255)
                                                                                             cGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAAt                                                       >  3:415214/1‑149 (MQ=255)
                                                                                                   gcAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATg                                                 >  1:604026/1‑149 (MQ=255)
                                                                                                     aGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGa                                               <  3:485545/149‑1 (MQ=255)
                                                                                                     aGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGa                                               <  4:32522/149‑1 (MQ=255)
                                                                                                     aGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGa                                               <  1:174251/149‑1 (MQ=255)
                                                                                                      gcgcCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGTCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGAc                                              <  1:692046/149‑1 (MQ=255)
                                                                                                      gcgcCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCTTTGTGGCTGGGGCGATTGTGCAATCGTATGGAc                                              >  4:669091/1‑149 (MQ=255)
                                                                                                        gcCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTa                                            >  2:913231/1‑149 (MQ=255)
                                                                                                           tCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTa                                                                                                 >  3:722112/1‑93 (MQ=255)
                                                                                                           tCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTa                                                                                                 <  4:722112/93‑1 (MQ=255)
                                                                                                                     ttCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGAttt                               >  1:821581/1‑149 (MQ=255)
                                                                                                                      tCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTa                              <  4:359199/149‑1 (MQ=255)
                                                                                                                           gcCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCttt                         >  4:306673/1‑149 (MQ=255)
                                                                                                                            ccTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTa                        >  2:150664/1‑149 (MQ=255)
                                                                                                                             cTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGcc                                                                               <  4:435013/93‑1 (MQ=255)
                                                                                                                             cTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGcc                                                                               >  3:435013/1‑93 (MQ=255)
                                                                                                                             cTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTAt                       >  3:377857/1‑149 (MQ=255)
                                                                                                                             cTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTAt                       <  3:529367/149‑1 (MQ=255)
                                                                                                                               gAAACTGTAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTa                              >  2:499791/1‑140 (MQ=255)
                                                                                                                               gAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTa                              <  1:499791/140‑1 (MQ=255)
                                                                                                                                 aaCTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGt                   >  4:771664/1‑149 (MQ=255)
                                                                                                                                   cTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCg                 >  3:786049/1‑149 (MQ=255)
                                                                                                                                    tGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAg                                                                                                >  4:190001/1‑69 (MQ=255)
                                                                                                                                    tGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAg                                                                                                <  3:190001/69‑1 (MQ=255)
                                                                                                                                    tGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGa                >  3:736507/1‑149 (MQ=255)
                                                                                                                                    tGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGa                >  4:870216/1‑149 (MQ=255)
                                                                                                                                        aTTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTgg            <  2:421178/149‑1 (MQ=255)
                                                                                                                                        aTTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTgg            >  2:280275/1‑149 (MQ=255)
                                                                                                                                        aTTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTgg            >  1:338385/1‑149 (MQ=255)
                                                                                                                                                  ggggAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTGGGCTCggcggc  <  2:281601/149‑1 (MQ=255)
                                                                                                                                                  ggggAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTGGGCTCggcggc  <  1:50992/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTGGGCTCGGCGGC  >  NC_000913/614710‑615004

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑