Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,169,260 Δ1 bp 5.1% 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. 5.1% 532.5 / 20.2 158intergenic (‑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 (81/68);  new base . (5/3);  total (87/71)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

Only 100 of 162 total aligned reads displayed.

CCACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTG  >  NC_000913/3169112‑3169406
                                                                                                                                                    |                                                                                                                                                   
ccACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGa                                                                                                                                                     <  4:882107/149‑1 (MQ=255)
ccACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGa                                                                                                                                                     >  2:731118/1‑149 (MQ=255)
 cacaTGATTAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATAGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCTGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGaa                                                                                                                                                    <  4:1235032/149‑1 (MQ=255)
 cacaTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGaa                                                                                                                                                    <  1:1019722/149‑1 (MQ=255)
  acaTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAAc                                                                                                                                                   >  4:549014/1‑149 (MQ=255)
     tGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAAcgcg                                                                                                                                                >  2:369560/1‑149 (MQ=255)
       aTAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑Tc                                                                                                                                             >  4:464658/1‑149 (MQ=255)
                 ccACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCgg                                                                                                                                   >  4:568933/1‑149 (MQ=255)
                  cACGCTGTAGTTGAGCAATTCAGGGCTACAGCAGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGt                                                                                                                                  >  2:684500/1‑149 (MQ=255)
                   aCGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTa                                                                                                                                 >  3:1481747/1‑149 (MQ=255)
                   aCGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTa                                                                                                                                 <  3:464658/149‑1 (MQ=255)
                              gAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGa                                                                                                                      <  4:1384014/149‑1 (MQ=255)
                                 cAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGAttt                                                                                                                   >  1:1010403/1‑149 (MQ=255)
                                  aaTGCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTa                                                                                                                  >  1:1297641/1‑149 (MQ=255)
                                   aTTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTa                                                                                                                  <  2:244623/148‑1 (MQ=255)
                                   aTTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTa                                                                                                                  >  1:244623/1‑148 (MQ=255)
                                    ttCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACt                                                                                                                >  4:1004337/1‑149 (MQ=255)
                                      cAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACttt                                                                                                              <  1:82824/149‑1 (MQ=255)
                                       aGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                             >  2:1059596/1‑149 (MQ=255)
                                        gggCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGc                                                                                                            <  2:850321/149‑1 (MQ=255)
                                         ggCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCg                                                                                                           <  3:1329485/149‑1 (MQ=255)
                                            tACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTAc                                                                                                                         >  1:799410/1‑132 (MQ=255)
                                            tACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTAc                                                                                                                         <  2:799410/132‑1 (MQ=255)
                                            tACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATGACGGATTTACTTTGCGGAt                                                                                                        >  4:1012097/1‑149 (MQ=255)
                                             aCAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTAc                                                                                                                 <  2:84026/139‑1 (MQ=255)
                                             aCAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGAGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATc                                                                                                       <  3:1025883/149‑1 (MQ=255)
                                                 cGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGtt                                                                                                   >  3:345926/1‑149 (MQ=255)
                                                  ggTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                  <  4:715138/149‑1 (MQ=255)
                                                  ggTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                  <  2:545393/149‑1 (MQ=255)
                                                      cGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACg                                                                                              <  1:1363126/149‑1 (MQ=255)
                                                          aaCATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCaa                                                                                          <  3:864617/149‑1 (MQ=255)
                                                             aTCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAttt                                                                                       >  1:851901/1‑149 (MQ=255)
                                                               cGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCa                                                                                     <  2:1450094/149‑1 (MQ=255)
                                                                gCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCaa                                                                                    >  4:1190937/1‑149 (MQ=255)
                                                                gCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCaa                                                                                    <  1:934020/149‑1 (MQ=255)
                                                                 ccACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAAc                                                                                   >  4:105154/1‑149 (MQ=255)
                                                                  cacaGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACa                                                                                  >  3:34385/1‑149 (MQ=255)
                                                                    caGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGa                                                                                >  2:834731/1‑149 (MQ=255)
                                                                     aGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGaa                                                                               <  4:601880/149‑1 (MQ=255)
                                                                       gTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAAtt                                                                             >  4:382345/1‑149 (MQ=255)
                                                                       ctATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTAATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAAtt                                                                             >  2:1280985/2‑149 (MQ=255)
                                                                        tATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTc                                                                            >  3:1225098/1‑149 (MQ=255)
                                                                          ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCa                                                                          <  1:573261/149‑1 (MQ=255)
                                                                          ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCa                                                                          <  1:1146216/149‑1 (MQ=255)
                                                                            gATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACg                                                                                                                        <  2:877824/101‑1 (MQ=255)
                                                                              tGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCt                                                                      >  3:592438/1‑149 (MQ=255)
                                                                              tGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCt                                                                      <  2:1185567/149‑1 (MQ=255)
                                                                               gggAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCtt                                                                     >  2:609123/1‑149 (MQ=255)
                                                                               gggAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCtt                                                                     >  1:1101626/1‑149 (MQ=255)
                                                                                  aaaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                  >  2:619633/1‑149 (MQ=255)
                                                                                   aaTCGATAACGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATc                                                                 >  1:1298411/1‑149 (MQ=255)
                                                                                    aTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCg                                                                >  2:177359/1‑149 (MQ=255)
                                                                                     tCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGa                                                               >  1:1315681/1‑149 (MQ=255)
                                                                                       gATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTATCAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACt                                                             >  2:1011244/1‑149 (MQ=255)
                                                                                         tAATGTTTACATCCCGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGG‑ACGCG‑TTCCcgg                                                                                                                                        <  1:1366269/71‑1 (MQ=255)
                                                                                         tAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGG‑ACGCG‑TTCCcgg                                                                                                                                        >  2:1366269/1‑71 (MQ=255)
                                                                                            tGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCAcc                                                        >  1:1279736/1‑149 (MQ=255)
                                                                                             gTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCAccc                                                       >  4:166096/1‑149 (MQ=255)
                                                                                             gTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCAccc                                                       >  3:955780/1‑149 (MQ=255)
                                                                                              tttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCt                                                      >  4:48497/1‑149 (MQ=255)
                                                                                                 aCATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTTCAGTCTATATCGACTTCACGCTGaa                                                   >  3:762548/1‑149 (MQ=255)
                                                                                                 aCATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGaa                                                   <  2:343302/149‑1 (MQ=255)
                                                                                                 aCATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGaa                                                   >  1:1054145/1‑149 (MQ=255)
                                                                                                  cATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAt                                                  <  1:172167/149‑1 (MQ=255)
                                                                                                     ccAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtctc                                               >  2:1133026/1‑149 (MQ=255)
                                                                                                      cAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTc                                                                       >  1:1241383/1‑124 (MQ=255)
                                                                                                        gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                            <  2:559831/149‑1 (MQ=255)
                                                                                                        gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                            <  3:266503/149‑1 (MQ=255)
                                                                                                        gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                            <  3:185612/149‑1 (MQ=255)
                                                                                                         tgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAcc                                           <  3:950261/149‑1 (MQ=255)
                                                                                                         tgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAcc                                           >  1:336144/1‑149 (MQ=255)
                                                                                                          gtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAg                                                                         >  3:1000669/1‑118 (MQ=255)
                                                                                                           tCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGt                                         >  2:1316630/1‑149 (MQ=255)
                                                                                                              ggTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGaa                                                                               >  2:1231939/1‑108 (MQ=255)
                                                                                                               gcttGTCATTTTGCCTCCCTGCGATATTATTCCGGGGAAAGCG‑TCCCGGGAGCGGTCACAATTACGGATTTACTTTGCGGATCTGTTTGGCGTCAATTTCAACAGAATTCCAGTCTTTATCGCCTTCACCCTGAATCTCAACCGTATCt                                     >  4:864095/3‑149 (MQ=255)
                                                                                                               gTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCt                                     >  3:1476051/1‑149 (MQ=255)
                                                                                                               gTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCt                                     >  1:86783/1‑149 (MQ=255)
                                                                                                                   gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                 <  1:689224/149‑1 (MQ=255)
                                                                                                                   gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                 <  4:1476051/149‑1 (MQ=255)
                                                                                                                      aTTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcg                              >  1:694629/1‑149 (MQ=255)
                                                                                                                         ttGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtca                           <  2:137456/149‑1 (MQ=255)
                                                                                                                           gCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcacc                         >  4:945519/1‑149 (MQ=255)
                                                                                                                             cTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgt                       >  4:1118393/1‑149 (MQ=255)
                                                                                                                              tCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCt                                                                      <  2:727684/101‑1 (MQ=255)
                                                                                                                              tCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgtc                      <  4:1243443/149‑1 (MQ=255)
                                                                                                                                ccTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgtcac                    <  1:177359/149‑1 (MQ=255)
                                                                                                                                  tGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                            <  3:774570/123‑1 (MQ=255)
                                                                                                                                  tGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGc                  >  3:995042/1‑149 (MQ=255)
                                                                                                                                  tGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGc                  >  4:63121/1‑149 (MQ=255)
                                                                                                                                     gAAATTATTCTCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGt               >  3:749021/1‑149 (MQ=255)
                                                                                                                                     gAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGAGGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGt               >  3:1033335/1‑149 (MQ=255)
                                                                                                                                     gAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGt               >  4:1358611/1‑149 (MQ=255)
                                                                                                                                      aaaTTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCTGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGtt              <  2:1297641/149‑1 (MQ=255)
                                                                                                                                       aaTTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCTGCGTGACGGTGACGCCgaaa             <  4:1314694/149‑4 (MQ=255)
                                                                                                                                          cgctCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAg          >  3:1426469/4‑149 (MQ=255)
                                                                                                                                             tCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGt       >  2:564018/1‑149 (MQ=255)
                                                                                                                                             tCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGt       >  3:1300524/1‑149 (MQ=255)
                                                                                                                                             tCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGt       >  1:1430369/1‑149 (MQ=255)
                                                                                                                                               ccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt     >  1:746117/1‑149 (MQ=255)
                                                                                                                                               ccGGGAACGCG‑TCCCGGGAGCGGCAACAATTACGTATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATTGACTTCATCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt     <  4:7737/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTG  >  NC_000913/3169112‑3169406

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