Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,561 G→C 8.1% intergenic (‑579/+24) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,5610GC 8.1% 377.6 / 21.5 113intergenic (‑579/+24)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (55/48);  new base C (4/5);  total (59/54)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-02

Only 100 of 116 total aligned reads displayed.

CGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACC  >  NC_000913/987414‑987707
                                                                                                                                                   |                                                                                                                                                  
cGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgc                                                                                                                                                   >  1:1049655/1‑149 (MQ=255)
   aaTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCCGCGGGCGCt                                                                                                                                                >  1:465092/1‑149 (MQ=255)
   aaTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCt                                                                                                                                                <  4:377251/149‑1 (MQ=255)
    aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTg                                                                                                                                               >  1:953810/1‑149 (MQ=255)
     tAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTgg                                                                                                                                              <  2:18423/149‑1 (MQ=255)
      aaGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGc                                                                                                                                             <  4:1074704/149‑1 (MQ=255)
       aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGct                                                                                                                                            <  3:132776/149‑1 (MQ=255)
          aTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtt                                                                                                                                         <  2:129728/149‑1 (MQ=255)
           tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttt                                                                                                                                        <  4:507520/149‑2 (MQ=255)
            gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtttt                                                                                                                                       <  2:701654/149‑3 (MQ=255)
             cGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGaa                                                                                                                                      <  1:100300/149‑1 (MQ=255)
              ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtttttt                                                                                                                                     >  1:40117/1‑145 (MQ=255)
               gATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttttttt                                                                                                                                    >  2:907832/1‑144 (MQ=255)
               gATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttttttt                                                                                                                                    <  2:40117/149‑6 (MQ=255)
                aTTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGAc                                                                                                                                   <  1:796139/149‑1 (MQ=255)
                aTTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAAAAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGAc                                                                                                                                   >  2:713098/1‑149 (MQ=255)
                 ttAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACg                                                                                                                                  <  3:569452/149‑1 (MQ=255)
                    gTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTa                                                                                                                               >  1:155330/1‑149 (MQ=255)
                    gTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTa                                                                                                                               >  2:452440/1‑149 (MQ=255)
                     tttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTaa                                                                                                                              >  3:878103/1‑149 (MQ=255)
                      ttttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTaaa                                                                                                                             <  1:452440/149‑1 (MQ=255)
                       tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAt                                                                                                                            <  2:953810/149‑1 (MQ=255)
                       tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAt                                                                                                                            >  3:918790/1‑149 (MQ=255)
                                 aaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTgg                                                                                                                  <  2:142476/149‑1 (MQ=255)
                                       gCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAc                                                                                                            <  1:344724/149‑1 (MQ=255)
                                           tGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAAAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                                                                        >  1:625629/1‑149 (MQ=255)
                                            gTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGa                                                                                                       <  4:538202/149‑1 (MQ=255)
                                             tAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                                                                      >  1:973831/1‑149 (MQ=255)
                                                tatTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGGGTTACGCGGAGTAc                                                                                                   <  3:78913/149‑1 (MQ=255)
                                                 atTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACg                                                                                                  >  4:119373/1‑149 (MQ=255)
                                                    ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                                                                               <  2:1049655/149‑1 (MQ=255)
                                                      gTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGa                                                                                                       <  3:581029/139‑1 (MQ=255)
                                                      gTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGa                                                                                                       >  4:581029/1‑139 (MQ=255)
                                                       tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGaa                                                                                            <  1:348164/149‑1 (MQ=255)
                                                        aaCAGGATGATCGTGATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAACTTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                                                                           <  1:569703/149‑1 (MQ=255)
                                                            ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                                                                       <  1:425649/149‑1 (MQ=255)
                                                             gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                                                                      >  2:218659/1‑149 (MQ=255)
                                                             gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                                                                      >  2:675341/1‑149 (MQ=255)
                                                               tgatCGTTATGGCGGGGGTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCCAAATTAGAAGCTGGCGTTACGCGGAGTACGGGGGAACGGAatca                                                                                    >  3:813079/1‑149 (MQ=255)
                                                                gatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcac                                                                                   <  1:1085916/149‑1 (MQ=255)
                                                                  tCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGTGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacat                                                                                 <  2:208497/149‑1 (MQ=255)
                                                                    gTTATGGCGGGGTTGCAAATAGATTGCTTTACTTAACAAACAGATAACTTGAAAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTTAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                                                                               >  4:846320/1‑133 (MQ=255)
                                                                    gTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                                                                               <  3:846320/133‑1 (MQ=255)
                                                                     ttATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatcac                                                                              >  2:12041/1‑149 (MQ=255)
                                                                     ttATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAAAAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGGGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGGGTCAGTGGGAACGGGatcacataac                                                                              >  4:140056/1‑149 (MQ=255)
                                                                         ggCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac                                                                          >  3:1026563/1‑149 (MQ=255)
                                                                             gggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTAACAGAAAAGATAAAAAAAGACCCAGCGGACCCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt                                                                      <  3:667514/149‑1 (MQ=255)
                                                                               gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTg                                                                    <  4:126524/149‑1 (MQ=255)
                                                                                ttGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                                                                   <  2:540950/149‑1 (MQ=255)
                                                                                 tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCa                                                                  <  4:634277/149‑1 (MQ=255)
                                                                                    aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                               <  3:945670/149‑1 (MQ=255)
                                                                                    aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGACTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                               >  3:116310/1‑149 (MQ=255)
                                                                                    aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATCAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                               >  1:787840/1‑149 (MQ=255)
                                                                                         gATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGtt                                                          <  1:922751/149‑1 (MQ=255)
                                                                                           ttGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAc                                                        <  1:713098/149‑1 (MQ=255)
                                                                                           ttGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAc                                                        <  3:371237/149‑1 (MQ=255)
                                                                                           ttGCTTGACTCAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAc                                                        >  3:952177/1‑149 (MQ=255)
                                                                                              cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                                                                       <  4:904056/115‑1 (MQ=255)
                                                                                              cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                                                                       >  3:904056/1‑115 (MQ=255)
                                                                                               ttGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTaa                                                    <  3:855773/149‑1 (MQ=255)
                                                                                               ttGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAATTACGTaa                                                    <  3:996883/149‑1 (MQ=255)
                                                                                                    ttAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAAt                                               >  4:202242/1‑149 (MQ=255)
                                                                                                    ttAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAAt                                               >  4:948028/1‑149 (MQ=255)
                                                                                                     tAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATc                                              <  3:984148/149‑1 (MQ=255)
                                                                                                        caaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGa                                           <  3:1069091/149‑1 (MQ=255)
                                                                                                         aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                          >  4:247972/1‑149 (MQ=255)
                                                                                                         aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                          >  3:601999/1‑149 (MQ=255)
                                                                                                         aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                          >  3:184978/1‑149 (MQ=255)
                                                                                                         aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                          >  4:5089/1‑149 (MQ=255)
                                                                                                           acaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGt                                        >  1:518311/1‑149 (MQ=255)
                                                                                                             aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTc                                      >  3:883667/1‑149 (MQ=255)
                                                                                                             aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTc                                      >  3:96849/1‑149 (MQ=255)
                                                                                                             aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTc                                      >  1:656427/1‑149 (MQ=255)
                                                                                                               aTAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAc                                                        <  4:561588/129‑1 (MQ=255)
                                                                                                                 aaCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCaaaa                                  <  4:883667/149‑1 (MQ=255)
                                                                                                                  aCTTGAGAGACAAGATAAAAAAAGTGCCAGCGGGGGCTGGCTCTTGAAGACGCTAAATTCGTAGCTGCCGTTACGCGGAGTACGTGGGACCGGAATCACATCACGTACGTTTTGCTCGCCAGTTCCGTAAGCAATCAGACTTTCAAAAc                                 <  4:614707/149‑1 (MQ=255)
                                                                                                                    ttGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTaa                                                    >  3:20592/1‑128 (MQ=255)
                                                                                                                    ttGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCa                               >  1:166879/1‑149 (MQ=255)
                                                                                                                     tGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCaa                              <  2:380870/149‑1 (MQ=255)
                                                                                                                      gACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAg                             <  3:530405/149‑1 (MQ=255)
                                                                                                                       aCAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGa                            <  3:537990/149‑1 (MQ=255)
                                                                                                                           aaaaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGa                        >  1:1046419/1‑149 (MQ=255)
                                                                                                                              aGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAAc                     >  1:665966/1‑149 (MQ=255)
                                                                                                                              aGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAAc                     <  3:420489/149‑1 (MQ=255)
                                                                                                                                  aaaaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGa                 >  3:317563/1‑149 (MQ=255)
                                                                                                                                   aaaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGaa                >  2:752970/1‑149 (MQ=255)
                                                                                                                                    ttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAAt               >  4:66422/4‑149 (MQ=255)
                                                                                                                                      aaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGTGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGc             >  1:96740/1‑149 (MQ=255)
                                                                                                                                      aaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGc             >  2:940146/1‑149 (MQ=255)
                                                                                                                                         gagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGa          >  1:263356/1‑149 (MQ=255)
                                                                                                                                          agCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGaa         >  3:32308/1‑149 (MQ=255)
                                                                                                                                          agCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGaa         <  2:518311/149‑1 (MQ=255)
                                                                                                                                           gCCAGCGGGCGCTGGCGCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAAc        >  4:255328/1‑149 (MQ=255)
                                                                                                                                            ccAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACa       >  4:1033957/1‑149 (MQ=255)
                                                                                                                                            ccAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACa       <  4:918790/149‑1 (MQ=255)
                                                                                                                                               gCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGtt                                                          <  4:32117/95‑1 (MQ=255)
                                                                                                                                               gCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGtt                                                          >  3:32117/1‑95 (MQ=255)
                                                                                                                                                cGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTAc   <  4:952177/149‑1 (MQ=255)
                                                                                                                                                 ggggGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTAcc  <  4:184978/149‑1 (MQ=255)
                                                                                                                                                 gggCGATGGCTCTGGAAGACGCTAAATTAGAAGCTGGCGTTACGCGCAGTACGTGGGAACGGAATCACATCCCGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAACACCTAAACCTGAATGCGGAACAGTAcc  <  4:143216/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACC  >  NC_000913/987414‑987707

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