Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,169,260 Δ1 bp 5.6% 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.6% 362.7 / 12.6 108intergenic (‑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 (44/57);  new base . (4/2);  total (48/60)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01

Only 100 of 110 total aligned reads displayed.

CACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTG  >  NC_000913/3169113‑3169404
                                                                                                                                                   |                                                                                                                                                 
cacaTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGaa                                                                                                                                                  >  3:257408/1‑149 (MQ=255)
 acaTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAAc                                                                                                                                                 >  4:465620/1‑149 (MQ=255)
     gATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑t                                                                                                                                            >  1:929479/1‑149 (MQ=255)
      aTAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑Tc                                                                                                                                           >  1:197370/1‑149 (MQ=255)
         aaTTTCGCCACGCTGTAGTTGAGCGATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCg                                                                                                                                        >  4:768382/1‑149 (MQ=255)
         aaTTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCg                                                                                                                                        >  1:36545/1‑149 (MQ=255)
          aTTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCgg                                                                                                                                       <  1:1004795/149‑1 (MQ=255)
             tCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGaa                                                                                                                                    <  2:155557/149‑2 (MQ=255)
             tCGCCACGCTGTAGTTGAGCAATTAAGGGCTAAAGCTGTGCGCCAACATTGCCACAGGTATTTATGGGAAATCGAGAATGTTTACATCCAGTGTCAGGTCTGTCAATTTGCCTCCCTGCGAAATTATTCCCGGG‑ACGCGTTCCCGGGaa                                                                                                                                    <  2:882216/149‑2 (MQ=255)
                ccACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGGACGCG‑TCCCGGGAGCgg                                                                                                                                 <  2:634508/149‑1 (MQ=255)
                ccACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCgg                                                                                                                                 >  1:520348/1‑149 (MQ=255)
                    gCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAAc                                                                                                                             <  4:915185/149‑1 (MQ=255)
                       gTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAAt                                                                                                                          <  4:118306/149‑1 (MQ=255)
                         aGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTa                                                                                                                        <  4:493010/149‑1 (MQ=255)
                         aGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTa                                                                                                                        <  4:438038/149‑1 (MQ=255)
                         aGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTa                                                                                                                        <  3:573334/149‑1 (MQ=255)
                          gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTAc                                                                                                                       <  4:897420/149‑1 (MQ=255)
                          gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTAc                                                                                                                       >  3:971196/1‑149 (MQ=255)
                          gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTAc                                                                                                                       <  2:407058/149‑1 (MQ=255)
                          gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTAc                                                                                                                       <  1:758602/149‑1 (MQ=255)
                          gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTAc                                                                                                                       >  1:254712/1‑149 (MQ=255)
                           ttGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACg                                                                                                                      >  2:341020/1‑149 (MQ=255)
                             gAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGa                                                                                                                    >  3:537115/1‑149 (MQ=255)
                             gAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGa                                                                                                                    <  3:63066/149‑1 (MQ=255)
                               gCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGAtt                                                                                                                  >  2:318126/1‑149 (MQ=255)
                                cAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACGATTACGGAttt                                                                                                                 <  3:757853/149‑1 (MQ=255)
                                 aaTTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTGCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTa                                                                                                                <  1:938039/149‑1 (MQ=255)
                                   ttCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACt                                                                                                              >  4:64523/1‑149 (MQ=255)
                                    tCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACtt                                                                                                             <  3:263647/149‑1 (MQ=255)
                                     cAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACttt                                                                                                            <  4:971196/149‑1 (MQ=255)
                                          cTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGa                                                                                                       <  2:543014/149‑1 (MQ=255)
                                           tACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGAt                                                                                                      <  2:197370/149‑1 (MQ=255)
                                           tACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGAt                                                                                                      <  3:834722/149‑1 (MQ=255)
                                            aCAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATc                                                                                                     <  2:696575/149‑1 (MQ=255)
                                              aGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTg                                                                                                   <  2:116732/149‑1 (MQ=255)
                                               gtggTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGGCAGGTGTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGt                                                                                                  <  2:401411/147‑1 (MQ=255)
                                               gCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTCTGCGGATCTGt                                                                                                  <  3:593598/149‑1 (MQ=255)
                                                    gCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGAc                                                                                             <  3:2043/149‑1 (MQ=255)
                                                        cAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCa                                                                                         <  1:89924/149‑1 (MQ=255)
                                                         aaCATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCaa                                                                                        <  1:341020/149‑1 (MQ=255)
                                                            aTCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAttt                                                                                     >  4:88725/1‑149 (MQ=255)
                                                               gCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCaa                                                                                  >  2:19345/1‑149 (MQ=255)
                                                               gCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCaa                                                                                  <  2:956866/149‑1 (MQ=255)
                                                                  acaGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAg                                                                               >  2:965166/1‑149 (MQ=255)
                                                                         ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCa                                                                        <  3:88725/149‑1 (MQ=255)
                                                                         ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCa                                                                        <  2:352558/149‑1 (MQ=255)
                                                                             tGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCt                                                                    >  4:943934/1‑149 (MQ=255)
                                                                               ggAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCttt                                                                  >  1:430878/1‑149 (MQ=255)
                                                                                 aaaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                <  3:640626/149‑1 (MQ=255)
                                                                                  aaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATc                                                               <  2:625101/149‑1 (MQ=255)
                                                                                   aTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGGCAATTTCAACAGAATTCCAGTCTTTATCg                                                              >  1:940948/1‑149 (MQ=255)
                                                                                    tCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGa                                                             >  4:308291/1‑149 (MQ=255)
                                                                                    tCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGa                                                             <  4:257408/149‑1 (MQ=255)
                                                                                     cGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGAc                                                            <  4:511274/149‑1 (MQ=255)
                                                                                     cGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGAc                                                            <  4:537115/149‑1 (MQ=255)
                                                                                      gATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGt                                                                                                  >  2:207332/1‑110 (MQ=255)
                                                                                      gATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGt                                                                                                  <  1:207332/110‑1 (MQ=255)
                                                                                       aTAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACtt                                                          >  4:461882/1‑149 (MQ=255)
                                                                                       aTAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACtt                                                          >  2:341842/1‑149 (MQ=255)
                                                                                           tGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCAcc                                                      >  3:155921/1‑149 (MQ=255)
                                                                                             tttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCt                                                    >  4:220008/1‑149 (MQ=255)
                                                                                             tttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCt                                                    >  4:779100/1‑149 (MQ=255)
                                                                                                 cATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAt                                                <  3:768382/149‑1 (MQ=255)
                                                                                                  aTCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtc                                               >  2:302261/1‑149 (MQ=255)
                                                                                                       gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                          <  4:935270/149‑1 (MQ=255)
                                                                                                       gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                          <  3:64523/149‑1 (MQ=255)
                                                                                                            aGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTAt                                     >  1:44582/1‑149 (MQ=255)
                                                                                                             ggTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTTAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATc                                    <  4:31768/149‑1 (MQ=255)
                                                                                                              gTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                <  3:1001751/88‑1 (MQ=255)
                                                                                                              gTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                >  4:1001751/1‑88 (MQ=255)
                                                                                                              gTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCt                                   >  1:721763/1‑149 (MQ=255)
                                                                                                                   tCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCgg                              >  3:873621/1‑149 (MQ=255)
                                                                                                                   tCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCgg                              >  2:443824/1‑149 (MQ=255)
                                                                                                                      ttttGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGc                                                                                                          >  1:513666/1‑70 (MQ=255)
                                                                                                                      ttttGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGc                                                                                                          <  2:513666/70‑1 (MQ=255)
                                                                                                                      ttttGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgt                           <  1:48826/149‑1 (MQ=255)
                                                                                                                      ttttGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgt                           <  4:166245/149‑1 (MQ=255)
                                                                                                                      ttttGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgt                           >  3:770172/1‑149 (MQ=255)
                                                                                                                             tCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTCGACGTCAATTTCAACAGAATTCCAGCCTTTATCGACTTCACCCTCAATCTCAACCGTATCTTTCGGcgtcaccgtc                    >  2:248679/1‑149 (MQ=255)
                                                                                                                              cccTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgtca                   <  2:990887/149‑1 (MQ=255)
                                                                                                                               ccTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgtcac                  >  3:437960/1‑149 (MQ=255)
                                                                                                                                    gAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                               <  1:246821/131‑1 (MQ=255)
                                                                                                                                    gAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                               >  2:246821/1‑131 (MQ=255)
                                                                                                                                     aaaTTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGtt            <  4:603367/149‑1 (MQ=255)
                                                                                                                                       attattCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTcc          <  3:779100/149‑1 (MQ=255)
                                                                                                                                       attattCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTcc          <  3:943934/149‑1 (MQ=255)
                                                                                                                                        ttattCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGt                                                                      <  4:499199/88‑1 (MQ=255)
                                                                                                                                        ttattCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGt                                                                      >  3:499199/1‑88 (MQ=255)
                                                                                                                                         tattCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAg        <  1:956850/149‑1 (MQ=255)
                                                                                                                                         cgctCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAg        >  4:426050/4‑149 (MQ=255)
                                                                                                                                          gctCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGc       >  1:288516/3‑149 (MQ=255)
                                                                                                                                           ttCCCGGGAACGCG‑TCCCCGGAGCGGTAACAATTACGGATGTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCg      >  1:312188/1‑149 (MQ=255)
                                                                                                                                            tCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGAc                                                            >  4:567962/1‑94 (MQ=255)
                                                                                                                                            tCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGAc                                                            <  3:567962/94‑1 (MQ=255)
                                                                                                                                            tCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGt     >  4:541570/1‑149 (MQ=255)
                                                                                                                                             cccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                          <  1:262676/111‑1 (MQ=255)
                                                                                                                                             cccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                          >  2:262676/1‑111 (MQ=255)
                                                                                                                                             cccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGtt    <  2:952700/149‑1 (MQ=255)
                                                                                                                                             cccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGtt    <  1:934115/149‑1 (MQ=255)
                                                                                                                                               cGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTtg  >  3:177224/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTG  >  NC_000913/3169113‑3169404

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