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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,169,2600A. 6.5% 315.0 / 14.9 92intergenic (‑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 (48/38);  new base . (3/3);  total (51/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

CACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGT  >  NC_000913/3169113‑3169405
                                                                                                                                                   |                                                                                                                                                  
cacaTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGa                                                                                                                                                    <  1:551951/148‑1 (MQ=255)
 acaTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGaa                                                                                                                                                   >  1:704300/1‑148 (MQ=255)
 acaTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAAc                                                                                                                                                  >  4:481527/1‑149 (MQ=255)
      aTAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑t                                                                                                                                             >  1:421428/1‑148 (MQ=255)
      aTAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑Tc                                                                                                                                            <  3:592983/149‑1 (MQ=255)
            ttCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGa                                                                                                                                      >  4:107666/1‑149 (MQ=255)
             tCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTACTCCCGGG‑ACGCGTTCCCGGGaa                                                                                                                                     <  3:729302/149‑2 (MQ=255)
              cGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGc                                                                                                                                    <  3:107666/149‑1 (MQ=255)
               gCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCg                                                                                                                                   <  3:685942/149‑1 (MQ=255)
               gCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCg                                                                                                                                   <  4:811457/149‑1 (MQ=255)
                 cACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCgg                                                                                                                                  <  2:101262/148‑1 (MQ=255)
                 cACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCgg                                                                                                                                  >  2:194801/1‑148 (MQ=255)
                 cACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGt                                                                                                                                 >  4:509322/1‑149 (MQ=255)
                   cGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTa                                                                                                                                <  2:839584/148‑1 (MQ=255)
                   cGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTa                                                                                                                                <  1:790725/148‑1 (MQ=255)
                    gCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAAc                                                                                                                              <  3:509322/149‑1 (MQ=255)
                        tAGTTGAGCATTTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAAt                                                                                                                           >  1:477625/1‑148 (MQ=255)
                          gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTa                                                                                                                         <  2:477625/148‑1 (MQ=255)
                          gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTAc                                                                                                                        <  3:248522/149‑1 (MQ=255)
                            tGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGAAGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACgg                                                                                                                      <  4:91582/149‑1 (MQ=255)
                              aGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGa                                                                                                                     <  1:817849/148‑1 (MQ=255)
                                  aTTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTAc                                                                                                                >  3:659898/1‑149 (MQ=255)
                                    tCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGG‑ACGCGTTCCCggg                                                                                                                                       <  2:385050/124‑1 (MQ=255)
                                    tCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGG‑ACGCGTTCCCggg                                                                                                                                       >  1:385050/1‑124 (MQ=255)
                                      aGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                            >  4:238090/1‑149 (MQ=255)
                                                   tGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGa                                                                                               >  4:79048/1‑149 (MQ=255)
                                                    gCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGa                                                                                               <  1:200802/148‑1 (MQ=255)
                                                         aaCATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCgg                                                                                                                                        <  3:664029/102‑1 (MQ=255)
                                                         aaCATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCACTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCgg                                                                                                                                        >  4:664029/1‑102 (MQ=255)
                                                          aCATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCaa                                                                                         <  1:526358/148‑1 (MQ=255)
                                                           cATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAtt                                                                                       <  3:813855/149‑1 (MQ=255)
                                                            aTCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAttt                                                                                      >  3:833423/1‑149 (MQ=255)
                                                                  acaGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAg                                                                                >  3:134603/1‑149 (MQ=255)
                                                                    aGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGa                                                                               <  2:421428/148‑1 (MQ=255)
                                                                      gTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGGAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAAtt                                                                            >  4:617194/1‑149 (MQ=255)
                                                                       tATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTCCAACAGAATTc                                                                           <  3:514202/149‑1 (MQ=255)
                                                                         ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                 >  4:550435/1‑125 (MQ=255)
                                                                         ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                 <  3:550435/125‑1 (MQ=255)
                                                                         ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCa                                                                         >  4:394459/1‑149 (MQ=255)
                                                                             tGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCa                                                                                          <  1:111410/128‑1 (MQ=255)
                                                                             tGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCa                                                                                          >  2:111410/1‑128 (MQ=255)
                                                                             tGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCt                                                                     >  4:800472/1‑149 (MQ=255)
                                                                                gAAATCTATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCttt                                                                   >  2:486659/1‑148 (MQ=255)
                                                                                   aTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCg                                                               >  4:408982/1‑149 (MQ=255)
                                                                                   aTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCg                                                               >  3:813577/1‑149 (MQ=255)
                                                                                    tCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCg                                                               >  2:874903/1‑148 (MQ=255)
                                                                                      gATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACt                                                            >  4:22309/1‑149 (MQ=255)
                                                                                       aTAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACtt                                                           >  4:589643/1‑149 (MQ=255)
                                                                                       aTAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACtt                                                           >  4:721981/1‑149 (MQ=255)
                                                                                        tAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACtt                                                           >  2:662200/1‑148 (MQ=255)
                                                                                         aaTGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTc                                                          >  1:612606/1‑148 (MQ=255)
                                                                                         aaTGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTc                                                          >  1:101594/1‑148 (MQ=255)
                                                                                            gTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTaa                                                                                                                               >  2:493463/1‑76 (MQ=255)
                                                                                            gTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTaa                                                                                                                               <  1:493463/76‑1 (MQ=255)
                                                                                              ttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCt                                                     >  2:765685/1‑148 (MQ=255)
                                                                                                 cATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCa                                                         >  2:711912/1‑141 (MQ=255)
                                                                                                 cATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCa                                                         <  1:711912/141‑1 (MQ=255)
                                                                                                 cATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTACCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAt                                                 >  3:236492/1‑149 (MQ=255)
                                                                                                   tCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtct                                               >  3:447639/1‑149 (MQ=255)
                                                                                                     cAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGtt                                                                                                  >  2:60834/1‑96 (MQ=255)
                                                                                                     cAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGtt                                                                                                  <  1:60834/96‑1 (MQ=255)
                                                                                                     cAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtctc                                              >  2:717573/1‑148 (MQ=255)
                                                                                                      aGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCaa                                            <  4:96204/149‑1 (MQ=255)
                                                                                                      aGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCaa                                            <  4:659898/149‑1 (MQ=255)
                                                                                                       gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                           <  3:37065/149‑1 (MQ=255)
                                                                                                       gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                           <  3:252405/149‑1 (MQ=255)
                                                                                                        tgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                           >  2:884897/1‑148 (MQ=255)
                                                                                                           cAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGt                                        <  1:583765/148‑1 (MQ=255)
                                                                                                           cAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGt                                        <  2:101594/148‑1 (MQ=255)
                                                                                                           cAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTa                                       <  3:238090/149‑1 (MQ=255)
                                                                                                             ggTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATc                                     >  4:442161/1‑149 (MQ=255)
                                                                                                               tttGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGAGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCt                                    >  2:318157/1‑148 (MQ=255)
                                                                                                                  gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTc                                 <  1:486659/148‑1 (MQ=255)
                                                                                                                  gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                <  3:481527/149‑1 (MQ=255)
                                                                                                                  gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                <  3:853937/149‑1 (MQ=255)
                                                                                                                  gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                >  4:180886/1‑149 (MQ=255)
                                                                                                                  gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                >  4:192716/1‑149 (MQ=255)
                                                                                                                   tCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCgg                               >  4:675576/1‑149 (MQ=255)
                                                                                                                     aTTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCt                                                                                                     >  4:390463/1‑77 (MQ=255)
                                                                                                                     aTTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCt                                                                                                     <  3:390463/77‑1 (MQ=255)
                                                                                                                              cccTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGa                                                                                               >  1:73875/1‑74 (MQ=255)
                                                                                                                              cccTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGa                                                                                               <  2:73875/74‑1 (MQ=255)
                                                                                                                               ccTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtctc                                              >  4:279922/1‑122 (MQ=255)
                                                                                                                               ccTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtctc                                              <  3:279922/122‑1 (MQ=255)
                                                                                                                                  gCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGc                 >  2:379055/1‑148 (MQ=255)
                                                                                                                                     aaaTTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGt              >  2:431143/1‑148 (MQ=255)
                                                                                                                                        ttattCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                 <  1:415581/94‑1 (MQ=255)
                                                                                                                                        ttattCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                 >  2:415581/1‑94 (MQ=255)
                                                                                                                                            tCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGt      >  4:196539/1‑149 (MQ=255)
                                                                                                                                             cccGGG‑ACGCGTTCCCGGGAGCGGTAACAATTATGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTc                                 >  1:99101/1‑121 (MQ=255)
                                                                                                                                             cccGGG‑ACGCGTTCCCGGGAGCGGTAACAATTATGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTc                                 <  2:99101/121‑1 (MQ=255)
                                                                                                                                              ccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt    >  3:281696/1‑149 (MQ=255)
                                                                                                                                              ccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt    <  3:408982/149‑1 (MQ=255)
                                                                                                                                                 ggAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTtgt  >  2:118666/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGT  >  NC_000913/3169113‑3169405

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: ‑