Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,169,260 Δ1 bp 9.5% intergenic (‑29/+24) ygiV ← / ← ygiW transcriptional repressor for mcbR biofilm gene/hydrogen peroxide and cadmium resistance periplasmic protein; stress‑induced OB‑fold protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,169,2600A. 9.5% 379.5 / 31.7 126intergenic (‑29/+24)ygiV/ygiWtranscriptional repressor for mcbR biofilm gene/hydrogen peroxide and cadmium resistance periplasmic protein; stress‑induced OB‑fold protein
Reads supporting (aligned to +/- strand):  ref base A (59/55);  new base . (8/4);  total (67/59)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 127 total aligned reads displayed.

CCACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGT  >  NC_000913/3169112‑3169405
                                                                                                                                                    |                                                                                                                                                  
ccACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTCTCAGGTTTGTCATTATGCCTCCCTGCGAAATTATTCCCGGGaa                                                                                                                                                   >  2:1444472/1‑150 (MQ=255)
 cacaTGATAAATATCGCCTCGCTGTAGTTGAGCAATTCTGGGCTACAGCGGTGCGGCACCATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAAc                                                                                                                                                  <  2:1162896/150‑1 (MQ=255)
   caTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAAcgc                                                                                                                                                >  1:1320583/1‑150 (MQ=255)
      gATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGGTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTc                                                                                                                                            >  1:784606/1‑150 (MQ=255)
        tAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTccc                                                                                                                                          <  2:876346/150‑1 (MQ=255)
           aTTTCGCCACGCTGTAGTTGAGCAAATCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGACGCGTTCCCggg                                                                                                                                       <  2:1296572/150‑1 (MQ=255)
              tCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGc                                                                                                                                    <  1:1574467/150‑1 (MQ=255)
               cGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGTGAAATTATTCCCGGGAACGCGTCCCGGGAGCg                                                                                                                                   >  2:49482/1‑150 (MQ=255)
                gCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCgg                                                                                                                                  >  1:780895/1‑150 (MQ=255)
                gCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCgg                                                                                                                                  <  2:686824/150‑1 (MQ=255)
                 ccACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGACGCGTTCCCggg                                                                                                                                       >  1:169273/1‑144 (MQ=255)
                 ccACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGACGCGTTCCCggg                                                                                                                                       <  2:169273/144‑1 (MQ=255)
                  cACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGACGCGTTCCCGGGaataa                                                                                                                                  >  2:935327/1‑144 (MQ=255)
                  cACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGACGCGTTCCCGGGaataa                                                                                                                                  <  1:935327/148‑5 (MQ=255)
                     gCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACa                                                                                                                             >  1:1664338/1‑150 (MQ=255)
                         tAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTa                                                                                                                         >  1:1285262/1‑150 (MQ=255)
                          aGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTAc                                                                                                                        >  2:1719991/1‑150 (MQ=255)
                          aGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTAc                                                                                                                        >  1:231846/1‑150 (MQ=255)
                             tGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGa                                                                                                                     >  1:99370/1‑150 (MQ=255)
                              gAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGAt                                                                                                                    >  1:1470236/1‑150 (MQ=255)
                                gCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGAttt                                                                                                                  <  1:841383/150‑1 (MQ=255)
                                 cAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACCTCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACCATTACGGATTTa                                                                                                                 >  1:1323871/1‑150 (MQ=255)
                                 cAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTa                                                                                                                 >  1:1582338/1‑150 (MQ=255)
                                  aaTTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTAc                                                                                                                >  2:873135/1‑150 (MQ=255)
                                     tCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACttt                                                                                                             <  2:888910/150‑1 (MQ=255)
                                      cAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGACCGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                            >  2:343183/1‑150 (MQ=255)
                                       aGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGt                                                                                                                                 >  2:108505/1‑128 (MQ=255)
                                       aGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGt                                                                                                                                 <  1:108505/128‑1 (MQ=255)
                                       aGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGc                                                                                                           >  2:1496926/1‑150 (MQ=255)
                                        gggCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCg                                                                                                          >  2:1548956/1‑150 (MQ=255)
                                         ggCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCgg                                                                                                         >  1:1347728/1‑150 (MQ=255)
                                            tACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCTGGAGCGGTAACAATTACGGATTTACTTTGCTGAta                                                                                                      >  2:422073/1‑149 (MQ=255)
                                            tACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATc                                                                                                      <  1:873135/150‑1 (MQ=255)
                                              cAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTg                                                                                                    <  2:231846/150‑1 (MQ=255)
                                                 cGGTGCGGCAACATCGTCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                 >  2:467917/1‑150 (MQ=255)
                                                 cGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                 >  1:779899/1‑150 (MQ=255)
                                                 cGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                 <  1:462887/150‑1 (MQ=255)
                                                   gTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGa                                                                                               <  1:529070/150‑1 (MQ=255)
                                                   gTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGa                                                                                               <  1:31016/150‑1 (MQ=255)
                                                     gCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACg                                                                                             <  1:1302169/150‑1 (MQ=255)
                                                     gCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACg                                                                                             <  2:1706383/150‑1 (MQ=255)
                                                       ggCAAGATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTc                                                                                           <  2:1288123/150‑1 (MQ=255)
                                                       ggCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTc                                                                                           <  1:191896/150‑1 (MQ=255)
                                                          aaCATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAt                                                                                        <  1:1247837/150‑1 (MQ=255)
                                                           aCATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAtt                                                                                       <  2:837200/150‑1 (MQ=255)
                                                               cGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGGCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCaa                                                                                   <  1:1556833/150‑1 (MQ=255)
                                                                 ccccaGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACa                                                                                 <  2:230051/147‑1 (MQ=255)
                                                                     aGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAAt                                                                             >  1:225470/1‑150 (MQ=255)
                                                                     aGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAAt                                                                             <  2:459193/150‑1 (MQ=255)
                                                                        tATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTcc                                                                          >  1:386320/1‑150 (MQ=255)
                                                                        tATTGATGGGAAATCAATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTcc                                                                          <  2:1017415/150‑1 (MQ=255)
                                                                           tGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGt                                                                       >  2:638102/1‑150 (MQ=255)
                                                                            gATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTc                                                                      <  2:1096936/150‑1 (MQ=255)
                                                                               gggAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCACTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATAACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCttt                                                                   <  2:889349/150‑1 (MQ=255)
                                                                                  taaTCGTTAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTTCGAAATTATTGCCTGGAAGGCGGCCGGGGATCGGTAACAATTACGGATTTACTTTGCGGTTCTGTTTGACGTCAATTTCATCAGAATTCCAGTCTTTATc                                                                <  1:1105586/149‑1 (MQ=255)
                                                                                  aaaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATc                                                                <  2:418424/150‑1 (MQ=255)
                                                                                   aaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCg                                                               <  1:1591486/150‑1 (MQ=255)
                                                                                    aTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGa                                                              <  1:1548956/150‑1 (MQ=255)
                                                                                      cGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACt                                                            <  1:1496926/150‑1 (MQ=255)
                                                                                            tGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCAccc                                                      >  1:32259/1‑150 (MQ=255)
                                                                                            tGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCAccc                                                      >  2:1423024/1‑150 (MQ=255)
                                                                                            tGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCAccc                                                      >  2:1460407/1‑150 (MQ=255)
                                                                                             gTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCt                                                     <  2:509889/150‑1 (MQ=255)
                                                                                             gTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCt                                                     >  2:613040/1‑150 (MQ=255)
                                                                                              tttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTg                                                    <  1:1574788/150‑1 (MQ=255)
                                                                                               ttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGa                                                   >  2:1362893/1‑150 (MQ=255)
                                                                                                tACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGaa                                                  <  1:613040/150‑1 (MQ=255)
                                                                                                 aCATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGTGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAt                                                 <  1:49482/150‑1 (MQ=255)
                                                                                                  cATCGAGTGTCAGGTTTGTCATTTTGACTCCCTGCGTTATTATTCCCGGGAAAGCGTCGCGGGAGCGTTAAGAATTACGGATTTACTTTGCGTATCTGTTTGACGTCAATTTCAACAGAATTCGAGTCTTTATCGACTTCAACCTGAAtc                                                <  1:952845/150‑1 (MQ=255)
                                                                                                  cATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtc                                                <  2:1320583/150‑1 (MQ=255)
                                                                                                    tCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtctc                                              <  1:638102/150‑1 (MQ=255)
                                                                                                       aGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGTGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACATAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                           <  1:1358834/150‑1 (MQ=255)
                                                                                                          gtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGt                                        >  2:742020/1‑150 (MQ=255)
                                                                                                          gtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGt                                        >  1:634022/1‑150 (MQ=255)
                                                                                                              ggTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCt                                    >  1:1612868/1‑150 (MQ=255)
                                                                                                                   gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTAgg                               >  1:1546231/1‑150 (MQ=255)
                                                                                                                     cATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcg                             <  2:1574682/150‑1 (MQ=255)
                                                                                                                     cATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcg                             >  2:624634/1‑150 (MQ=255)
                                                                                                                      aTTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgt                            >  1:526914/1‑150 (MQ=255)
                                                                                                                          tGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcacc                        >  1:1289387/1‑150 (MQ=255)
                                                                                                                           gCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccg                       <  1:1460407/150‑1 (MQ=255)
                                                                                                                           gCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAACTCAACCGTATCTTTCGGcgtcaccg                       >  1:1029023/1‑150 (MQ=255)
                                                                                                                            ccTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgt                      >  2:703922/1‑150 (MQ=255)
                                                                                                                            ccTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgt                      >  1:39204/1‑150 (MQ=255)
                                                                                                                                ccTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACg                  >  2:164727/1‑150 (MQ=255)
                                                                                                                                ccTGCGAAATTAATCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACg                  <  1:797319/150‑1 (MQ=255)
                                                                                                                                 cTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCTGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTAACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGc                 >  2:245221/1‑150 (MQ=255)
                                                                                                                                   gCGAAATTATTCCCGGGACCGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCg               >  2:113405/1‑150 (MQ=255)
                                                                                                                                      aaaTTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTc            >  2:1175167/1‑150 (MQ=255)
                                                                                                                                       aaTTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAAGAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTcc           <  2:315630/150‑1 (MQ=255)
                                                                                                                                          cgctCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGc        >  1:352401/4‑150 (MQ=255)
                                                                                                                                           gctCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATc                                     <  2:1009473/118‑1 (MQ=255)
                                                                                                                                           gctCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCGGAATCTCACCCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCg       >  1:986442/3‑150 (MQ=255)
                                                                                                                                           cctCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATc                                     >  1:1009473/3‑120 (MQ=255)
                                                                                                                                           attCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGAGTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCCCCGTCACGCCGGGCCAGCg       >  1:754402/1‑150 (MQ=255)
                                                                                                                                            ctCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGt      >  2:871579/2‑150 (MQ=255)
                                                                                                                                             tCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGtt     <  1:827223/150‑1 (MQ=255)
                                                                                                                                             tCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTCTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGtt     <  2:223954/150‑1 (MQ=255)
                                                                                                                                              cccGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt    >  1:1624514/1‑150 (MQ=255)
                                                                                                                                                cGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATCTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTtgt  <  2:69741/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGT  >  NC_000913/3169112‑3169405

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