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

CACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGT  >  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)
      aTAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGt                                                                                                                                             >  1:421428/1‑148 (MQ=255)
      aTAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTc                                                                                                                                            <  3:592983/149‑1 (MQ=255)
            ttCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGa                                                                                                                                      >  4:107666/1‑149 (MQ=255)
             tCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTACTCCCGGGACGCGTTCCCGGGaa                                                                                                                                     <  3:729302/149‑2 (MQ=255)
              cGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGc                                                                                                                                    <  3:107666/149‑1 (MQ=255)
               gCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCg                                                                                                                                   <  3:685942/149‑1 (MQ=255)
               gCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCg                                                                                                                                   <  4:811457/149‑1 (MQ=255)
                 cACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCgg                                                                                                                                  <  2:101262/148‑1 (MQ=255)
                 cACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCgg                                                                                                                                  >  2:194801/1‑148 (MQ=255)
                 cACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGt                                                                                                                                 >  4:509322/1‑149 (MQ=255)
                   cGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTa                                                                                                                                <  2:839584/148‑1 (MQ=255)
                   cGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTa                                                                                                                                <  1:790725/148‑1 (MQ=255)
                    gCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAAc                                                                                                                              <  3:509322/149‑1 (MQ=255)
                        tAGTTGAGCATTTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAAt                                                                                                                           >  1:477625/1‑148 (MQ=255)
                          gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTa                                                                                                                         <  2:477625/148‑1 (MQ=255)
                          gTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTAc                                                                                                                        <  3:248522/149‑1 (MQ=255)
                            tGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGAAGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACgg                                                                                                                      <  4:91582/149‑1 (MQ=255)
                              aGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGa                                                                                                                     <  1:817849/148‑1 (MQ=255)
                                  aTTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTAc                                                                                                                >  3:659898/1‑149 (MQ=255)
                                    tCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGACGCGTTCCCggg                                                                                                                                       <  2:385050/124‑1 (MQ=255)
                                    tCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGACGCGTTCCCggg                                                                                                                                       >  1:385050/1‑124 (MQ=255)
                                      aGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                            >  4:238090/1‑149 (MQ=255)
                                                   tGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGa                                                                                               >  4:79048/1‑149 (MQ=255)
                                                    gCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGa                                                                                               <  1:200802/148‑1 (MQ=255)
                                                         aaCATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCgg                                                                                                                                        <  3:664029/102‑1 (MQ=255)
                                                         aaCATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCACTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCgg                                                                                                                                        >  4:664029/1‑102 (MQ=255)
                                                          aCATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCaa                                                                                         <  1:526358/148‑1 (MQ=255)
                                                           cATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAtt                                                                                       <  3:813855/149‑1 (MQ=255)
                                                            aTCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAttt                                                                                      >  3:833423/1‑149 (MQ=255)
                                                                  acaGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAg                                                                                >  3:134603/1‑149 (MQ=255)
                                                                    aGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGa                                                                               <  2:421428/148‑1 (MQ=255)
                                                                      gTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGGAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAAtt                                                                            >  4:617194/1‑149 (MQ=255)
                                                                       tATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTCCAACAGAATTc                                                                           <  3:514202/149‑1 (MQ=255)
                                                                         ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                 >  4:550435/1‑125 (MQ=255)
                                                                         ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                 <  3:550435/125‑1 (MQ=255)
                                                                         ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCa                                                                         >  4:394459/1‑149 (MQ=255)
                                                                             tGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCa                                                                                          <  1:111410/128‑1 (MQ=255)
                                                                             tGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCa                                                                                          >  2:111410/1‑128 (MQ=255)
                                                                             tGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCt                                                                     >  4:800472/1‑149 (MQ=255)
                                                                                gAAATCTATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCttt                                                                   >  2:486659/1‑148 (MQ=255)
                                                                                   aTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCg                                                               >  4:408982/1‑149 (MQ=255)
                                                                                   aTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCg                                                               >  3:813577/1‑149 (MQ=255)
                                                                                    tCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCg                                                               >  2:874903/1‑148 (MQ=255)
                                                                                      gATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACt                                                            >  4:22309/1‑149 (MQ=255)
                                                                                       aTAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACtt                                                           >  4:589643/1‑149 (MQ=255)
                                                                                       aTAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACtt                                                           >  4:721981/1‑149 (MQ=255)
                                                                                        tAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACtt                                                           >  2:662200/1‑148 (MQ=255)
                                                                                         aaTGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTc                                                          >  1:612606/1‑148 (MQ=255)
                                                                                         aaTGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTc                                                          >  1:101594/1‑148 (MQ=255)
                                                                                            gTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTaa                                                                                                                               >  2:493463/1‑76 (MQ=255)
                                                                                            gTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTaa                                                                                                                               <  1:493463/76‑1 (MQ=255)
                                                                                              ttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCt                                                     >  2:765685/1‑148 (MQ=255)
                                                                                                 cATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCa                                                         >  2:711912/1‑141 (MQ=255)
                                                                                                 cATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCa                                                         <  1:711912/141‑1 (MQ=255)
                                                                                                 cATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTACCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAt                                                 >  3:236492/1‑149 (MQ=255)
                                                                                                   tCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtct                                               >  3:447639/1‑149 (MQ=255)
                                                                                                     cAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGtt                                                                                                  >  2:60834/1‑96 (MQ=255)
                                                                                                     cAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGtt                                                                                                  <  1:60834/96‑1 (MQ=255)
                                                                                                     cAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtctc                                              >  2:717573/1‑148 (MQ=255)
                                                                                                      aGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCaa                                            <  4:96204/149‑1 (MQ=255)
                                                                                                      aGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCaa                                            <  4:659898/149‑1 (MQ=255)
                                                                                                       gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                           <  3:37065/149‑1 (MQ=255)
                                                                                                       gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                           <  3:252405/149‑1 (MQ=255)
                                                                                                        tgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                           >  2:884897/1‑148 (MQ=255)
                                                                                                           cAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGt                                        <  1:583765/148‑1 (MQ=255)
                                                                                                           cAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGt                                        <  2:101594/148‑1 (MQ=255)
                                                                                                           cAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTa                                       <  3:238090/149‑1 (MQ=255)
                                                                                                             ggTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATc                                     >  4:442161/1‑149 (MQ=255)
                                                                                                               tttGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGAGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCt                                    >  2:318157/1‑148 (MQ=255)
                                                                                                                  gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTc                                 <  1:486659/148‑1 (MQ=255)
                                                                                                                  gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                <  3:481527/149‑1 (MQ=255)
                                                                                                                  gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                <  3:853937/149‑1 (MQ=255)
                                                                                                                  gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                >  4:180886/1‑149 (MQ=255)
                                                                                                                  gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                >  4:192716/1‑149 (MQ=255)
                                                                                                                   tCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCgg                               >  4:675576/1‑149 (MQ=255)
                                                                                                                     aTTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCt                                                                                                     >  4:390463/1‑77 (MQ=255)
                                                                                                                     aTTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCt                                                                                                     <  3:390463/77‑1 (MQ=255)
                                                                                                                              cccTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGa                                                                                               >  1:73875/1‑74 (MQ=255)
                                                                                                                              cccTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGa                                                                                               <  2:73875/74‑1 (MQ=255)
                                                                                                                               ccTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtctc                                              >  4:279922/1‑122 (MQ=255)
                                                                                                                               ccTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAAtctc                                              <  3:279922/122‑1 (MQ=255)
                                                                                                                                  gCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGc                 >  2:379055/1‑148 (MQ=255)
                                                                                                                                     aaaTTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGt              >  2:431143/1‑148 (MQ=255)
                                                                                                                                        ttattCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                 <  1:415581/94‑1 (MQ=255)
                                                                                                                                        ttattCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                 >  2:415581/1‑94 (MQ=255)
                                                                                                                                            tCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGt      >  4:196539/1‑149 (MQ=255)
                                                                                                                                             cccGGGACGCGTTCCCGGGAGCGGTAACAATTATGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTc                                 >  1:99101/1‑121 (MQ=255)
                                                                                                                                             cccGGGACGCGTTCCCGGGAGCGGTAACAATTATGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTc                                 <  2:99101/121‑1 (MQ=255)
                                                                                                                                              ccGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt    >  3:281696/1‑149 (MQ=255)
                                                                                                                                              ccGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt    <  3:408982/149‑1 (MQ=255)
                                                                                                                                                 ggAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTtgt  >  2:118666/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGT  >  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: