Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 614,859 G→A 5.5% E235K (GAA→AAA)  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,8590GA 5.5% 239.7 / 3.2 73E235K (GAA→AAA) entFenterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (38/31);  new base A (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 = 5.98e-01

GAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTGGGCTCGGCGGC  >  NC_000913/614712‑615004
                                                                                                                                                   |                                                                                                                                                 
gAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGa                                                                                                                                                  <  2:241817/149‑1 (MQ=255)
 aaaGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGaa                                                                                                                                                 <  1:899266/149‑1 (MQ=255)
 aaaGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGaa                                                                                                                                                 <  3:250430/149‑1 (MQ=255)
    gCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTc                                                                                                                                              >  2:306231/1‑149 (MQ=255)
            tGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTCCCCGTCGGTAAATTccag                                                                                                                                           >  3:359510/1‑142 (MQ=255)
            tGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTCCCCGTCGGTAAATTccag                                                                                                                                           <  4:359510/144‑3 (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                                                                                               <  3:97992/149‑1 (MQ=255)
                                                   cAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGc                                                                                               >  4:68008/1‑149 (MQ=255)
                                                    aaCTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGcc                                                                                              >  2:692046/1‑149 (MQ=255)
                                                    aaCTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGcc                                                                                              >  3:703924/1‑149 (MQ=255)
                                                       tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                           >  4:429211/1‑149 (MQ=255)
                                                       tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                           <  4:363666/149‑1 (MQ=255)
                                                       tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                           >  3:311761/1‑149 (MQ=255)
                                                       tGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCtt                                                                                           <  2:289187/149‑1 (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:432426/149‑1 (MQ=255)
                                                                             tcCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTgg                                                                     <  4:435379/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                                               <  1:174251/149‑1 (MQ=255)
                                                                                                   aGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGa                                               <  4:32522/149‑1 (MQ=255)
                                                                                                   aGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGa                                               <  3:485545/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                                                                                                <  3:190001/69‑1 (MQ=255)
                                                                                                                                  tGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAg                                                                                                >  4:190001/1‑69 (MQ=255)
                                                                                                                                  tGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGa                >  3:736507/1‑149 (MQ=255)
                                                                                                                                  tGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGa                >  4:870216/1‑149 (MQ=255)
                                                                                                                                      aTTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTgg            >  2:280275/1‑149 (MQ=255)
                                                                                                                                      aTTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTgg            >  1:338385/1‑149 (MQ=255)
                                                                                                                                      aTTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTgg            <  2:421178/149‑1 (MQ=255)
                                                                                                                                                ggggAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTGGGCTCggcggc  <  1:50992/149‑1 (MQ=255)
                                                                                                                                                ggggAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTGGGCTCggcggc  <  2:281601/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAAAGCGAAGCCTGGCAGCGTGATGCGGCATTCTGGGCAGAACAGCGTCGTCAACTGCCGCCGCCCGCGTCACTTTCTCCGGCACCTTTACCGGGGCGCAGCGCCTCGGCAGATATTCTGCGCCTGAAACTGGAATTTACCGACGGGGAATTCCGCCAGCTGGCTACGCAACTTTCAGGTGTGCAGCGTACCGATTTAGCCCTTGCGCTGGCAGCCTTGTGGCTGGGGCGATTGTGCAATCGTATGGACTACGCCGCCGGATTTATCTTTATGCGTCGACTGGGCTCGGCGGC  >  NC_000913/614712‑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: ‑