Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,169,260 Δ1 bp 7.4% 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. 7.4% 318.1 / 17.7 95intergenic (‑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 (41/47);  new base . (3/4);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

CACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTG  >  NC_000913/3169113‑3169406
                                                                                                                                                   |                                                                                                                                                   
cacaTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGaa                                                                                                                                                    <  3:924409/149‑1 (MQ=255)
 acaTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGaa                                                                                                                                                    >  2:820541/1‑148 (MQ=255)
  caTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAAc                                                                                                                                                   >  1:893118/1‑148 (MQ=255)
     gATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑t                                                                                                                                              >  3:437937/1‑149 (MQ=255)
      aTAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑t                                                                                                                                              >  2:743459/1‑148 (MQ=255)
       tAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑Tcc                                                                                                                                            >  4:20647/1‑149 (MQ=255)
          aTTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGG‑ACGCG‑ttcccg                                                                                                                                         <  2:651316/148‑6 (MQ=255)
            ttCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCggg                                                                                                                                        <  1:275408/148‑1 (MQ=255)
            ttCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGG‑ACGCGTTCCCggg                                                                                                                                        <  1:61748/148‑1 (MQ=255)
             tCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAg                                                                                                                                      <  3:166288/149‑1 (MQ=255)
                ccACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGG‑ACGCGTTCCCGGGa                                                                                                                                       >  4:41969/1‑145 (MQ=255)
                ccACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGG‑ACGCGTTCCCGGGa                                                                                                                                       <  3:41969/145‑1 (MQ=255)
                 cACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGt                                                                                                                                  >  4:763279/1‑149 (MQ=255)
                 cACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGt                                                                                                                                  >  3:278845/1‑149 (MQ=255)
                          gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTa                                                                                                                          <  2:422240/148‑1 (MQ=255)
                          gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTAc                                                                                                                         <  3:45973/149‑1 (MQ=255)
                           ttGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTAc                                                                                                                         >  2:640018/1‑148 (MQ=255)
                              aGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGa                                                                                                                      <  1:743459/148‑1 (MQ=255)
                               gCAATTCAGGACTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGAtt                                                                                                                    <  3:470625/149‑1 (MQ=255)
                                 aaTTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGAttt                                                                                                                   <  1:592139/148‑1 (MQ=255)
                                   ttCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACt                                                                                                                >  4:847059/1‑149 (MQ=255)
                                     cAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACttt                                                                                                              <  4:425358/149‑1 (MQ=255)
                                      aGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTTTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                             >  4:716027/1‑149 (MQ=255)
                                      aGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                             >  4:26729/1‑149 (MQ=255)
                                       gggCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                             >  1:548858/1‑148 (MQ=255)
                                       gggCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                             <  1:867368/148‑1 (MQ=255)
                                       gggCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGc                                                                                                            <  3:847059/149‑1 (MQ=255)
                                            aCAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATc                                                                                                       <  3:904831/149‑1 (MQ=255)
                                             cAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCt                                                                                                      <  4:427335/149‑1 (MQ=255)
                                                cGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGtt                                                                                                   <  3:716027/149‑1 (MQ=255)
                                                 ggTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGtt                                                                                                   <  1:673973/148‑1 (MQ=255)
                                                 ggTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGtt                                                                                                   <  2:775971/148‑1 (MQ=255)
                                                       gCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGt                                                                                             >  1:825316/1‑148 (MQ=255)
                                                             tCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAttt                                                                                       >  2:241164/1‑148 (MQ=255)
                                                             tCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTAAGGATTTACTTTGCGGATCTGTTTGACGTCAAttt                                                                                       >  1:53976/1‑148 (MQ=255)
                                                               gCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGAGGTCAATTTCaa                                                                                    <  3:763279/149‑1 (MQ=255)
                                                                ccACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCaa                                                                                    <  2:499010/148‑1 (MQ=255)
                                                                ccACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAAc                                                                                   <  3:843188/149‑1 (MQ=255)
                                                                 cacaGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACa                                                                                  <  4:441731/149‑1 (MQ=255)
                                                                  acaGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAg                                                                                 >  3:308272/1‑149 (MQ=255)
                                                                   caGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAg                                                                                 >  2:195419/1‑148 (MQ=255)
                                                                    aGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGa                                                                                <  2:786437/148‑1 (MQ=255)
                                                                         ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAATGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCa                                                                          <  4:779375/149‑1 (MQ=255)
                                                                         ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTcc                                                                           <  1:820541/148‑1 (MQ=255)
                                                                         ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCa                                                                          <  3:177308/149‑1 (MQ=255)
                                                                           gATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGt                                                                        <  4:624802/149‑1 (MQ=255)
                                                                            aTGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTc                                                                       >  4:899495/1‑149 (MQ=255)
                                                                             tGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCt                                                                      >  3:82063/1‑149 (MQ=255)
                                                                                 aaaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAt                                                                                         <  4:676976/126‑1 (MQ=255)
                                                                                 aaaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAt                                                                                         >  3:676976/1‑126 (MQ=255)
                                                                                 aaaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTa                                                                   <  2:825316/148‑1 (MQ=255)
                                                                                  aaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                  <  1:280407/148‑1 (MQ=255)
                                                                                  aaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATc                                                                 <  3:1986/149‑1 (MQ=255)
                                                                                     cGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGa                                                               >  2:52415/1‑148 (MQ=255)
                                                                                     cGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGa                                                               <  1:305336/148‑1 (MQ=255)
                                                                                     cGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGAc                                                              <  3:675697/149‑1 (MQ=255)
                                                                                       aTAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACtt                                                            >  4:82600/1‑149 (MQ=255)
                                                                                        tAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACtt                                                            >  1:319965/1‑148 (MQ=255)
                                                                                        tAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACtt                                                            >  2:155849/1‑148 (MQ=255)
                                                                                        tAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACtt                                                            >  1:722686/1‑148 (MQ=255)
                                                                                         aaTGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTc                                                           >  1:884246/1‑148 (MQ=255)
                                                                                              ttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCggg                                                                                                                                        >  2:751789/1‑66 (MQ=255)
                                                                                              ttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCggg                                                                                                                                        <  1:751789/66‑1 (MQ=255)
                                                                                                 cATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAt                                                  <  3:26729/149‑1 (MQ=255)
                                                                                                   tCCAGAGTGAGGTTTGTCATTTTGCCTCCCTTAGAAAATATTCACGGGAACGCT‑TTCTGGCAGCTGTACTAGTTACGATTTTAGTTTTTGGATCTGTTGGACGTAAGTTTCAACAGAATTACAATCTTTACAGACCTCACCTTTAAtc                                                 >  2:693918/1‑148 (MQ=255)
                                                                                                     cAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtctc                                               >  1:114750/1‑148 (MQ=255)
                                                                                                          tCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGt                                         >  3:126193/1‑149 (MQ=255)
                                                                                                           cAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGt                                         >  1:310186/1‑148 (MQ=255)
                                                                                                             ggTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAttt                                                                                       <  3:354009/100‑1 (MQ=255)
                                                                                                             ggTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAttt                                                                                       >  4:354009/1‑100 (MQ=255)
                                                                                                               tttGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCt                                     >  2:655122/1‑148 (MQ=255)
                                                                                                                 tGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTc                                  <  3:375140/149‑1 (MQ=255)
                                                                                                                  gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                 <  3:82600/149‑1 (MQ=255)
                                                                                                                  gTCATTTTGCCTCCCTGCGAAATTATGCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCTTCCTCACCCTGAATCTCAACAGTATCTTTCg                                 <  4:487648/149‑1 (MQ=255)
                                                                                                                        ttGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtc                            <  1:155849/148‑1 (MQ=255)
                                                                                                                           ccTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcacc                         >  2:1664/1‑148 (MQ=255)
                                                                                                                             tCCCTGCGAAATTATTCCTGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGa                                                                                                <  2:792287/75‑1 (MQ=255)
                                                                                                                             tCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGa                                                                                                >  1:792287/1‑75 (MQ=255)
                                                                                                                             tCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgtc                      >  4:750027/1‑149 (MQ=255)
                                                                                                                             tCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgt                       >  1:512134/1‑148 (MQ=255)
                                                                                                                                  gCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGc                  >  2:273006/1‑148 (MQ=255)
                                                                                                                                    gAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCTACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGt               <  4:317880/149‑1 (MQ=255)
                                                                                                                                    gAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                            <  3:530572/120‑1 (MQ=255)
                                                                                                                                    gAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                            >  4:530572/1‑120 (MQ=255)
                                                                                                                                     aaaTTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGt               <  2:884246/148‑1 (MQ=255)
                                                                                                                                        ttattCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTcc            >  1:368804/1‑148 (MQ=255)
                                                                                                                                        ttattCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCa           <  4:82063/149‑1 (MQ=255)
                                                                                                                                         tattCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCa           >  1:944930/1‑148 (MQ=255)
                                                                                                                                           ctCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGa                                                               >  4:483301/2‑94 (MQ=255)
                                                                                                                                           ctCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGa                                                               <  3:483301/93‑1 (MQ=255)
                                                                                                                                             cccGGGAACGCG‑TCCCGGGAGCGGTAACAATTAAGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGt       <  2:53976/148‑1 (MQ=255)
                                                                                                                                              ccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt     <  3:899495/149‑1 (MQ=255)
                                                                                                                                              ccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGtt      <  1:52415/148‑1 (MQ=255)
                                                                                                                                              ccGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGtt      >  2:954444/1‑148 (MQ=255)
                                                                                                                                               cGGGAACGCG‑TGCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTtg    >  3:541789/1‑149 (MQ=255)
                                                                                                                                                 ggAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTtgtg  <  3:915988/149‑1 (MQ=255)
                                                                                                                                                 ggAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGGTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTtgt   <  2:548858/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTG  >  NC_000913/3169113‑3169406

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