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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,169,2600A. 7.0% 298.6 / 15.7 86intergenic (‑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 (43/37);  new base . (4/2);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CCACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTGGT  >  NC_000913/3169112‑3169408
                                                                                                                                                    |                                                                                                                                                     
ccACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGa                                                                                                                                                       >  2:652230/1‑149 (MQ=255)
  acaTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCGCAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAAc                                                                                                                                                     >  1:271560/1‑149 (MQ=255)
    aTGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAAcgc                                                                                                                                                   <  3:376598/149‑1 (MQ=255)
           aTTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCgg                                                                                                                                           <  3:444065/149‑1 (MQ=255)
           aTTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCgg                                                                                                                                           >  2:452536/1‑149 (MQ=255)
            tttCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCggg                                                                                                                                          >  2:83233/1‑149 (MQ=255)
                gCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCg                                                                                                                                      >  4:413197/1‑149 (MQ=255)
                 ccACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCgg                                                                                                                                     >  2:148972/1‑149 (MQ=255)
                       tGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACaa                                                                                                                               >  2:153243/1‑149 (MQ=255)
                       tGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACaa                                                                                                                               <  1:654101/149‑1 (MQ=255)
                        gTAGTCGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAAt                                                                                                                              >  3:552792/1‑149 (MQ=255)
                              gAGCAATTCAGGGCTACAGCGGTGCGTCAACATCGCCAGAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGGTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGa                                                                                                                        <  2:392077/149‑1 (MQ=255)
                              gAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGa                                                                                                                        <  1:652230/149‑1 (MQ=255)
                              gAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGa                                                                                                                        <  2:610621/149‑1 (MQ=255)
                                    ttCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACt                                                                                                                  >  1:489746/1‑149 (MQ=255)
                                    ttCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACt                                                                                                                  >  1:687591/1‑149 (MQ=255)
                                       aGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                               >  3:766966/1‑149 (MQ=255)
                                         ggCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCg                                                                                                             <  4:19938/149‑1 (MQ=255)
                                            tACAGCGGTGCGGCAACATCGCCGCAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGAt                                                                                                          <  2:271560/149‑1 (MQ=255)
                                                gCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGt                                                                                                      <  1:436933/149‑1 (MQ=255)
                                                gCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGt                                                                                                      <  3:669416/149‑1 (MQ=255)
                                                  ggTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                    <  4:798671/149‑1 (MQ=255)
                                                  ggTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                    <  1:238798/149‑1 (MQ=255)
                                                    tGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCTGGAGCGGTAACAATTACGGATTCACTTTGCGGATCTGTTTGa                                                                                                  <  3:472027/149‑1 (MQ=255)
                                                       ggCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACTCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGt                                                                                               <  2:688535/149‑1 (MQ=255)
                                                       ggCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGt                                                                                               <  2:68734/149‑1 (MQ=255)
                                                             aTCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAttt                                                                                         >  1:558171/1‑149 (MQ=255)
                                                              tCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTc                                                                                        <  3:190053/149‑1 (MQ=255)
                                                                gCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCaa                                                                                      >  4:800950/1‑149 (MQ=255)
                                                                gCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCaa                                                                                      <  3:224570/149‑1 (MQ=255)
                                                                 ccACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAAc                                                                                     <  1:148972/149‑1 (MQ=255)
                                                                   acaGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAg                                                                                   >  3:639040/1‑149 (MQ=255)
                                                                   acaGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAg                                                                                   >  1:508673/1‑149 (MQ=255)
                                                                      ggTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAAt                                                                                >  2:504734/1‑149 (MQ=255)
                                                                          ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCa                                                                            >  2:116693/1‑149 (MQ=255)
                                                                            gATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGt                                                                          <  1:83233/149‑1 (MQ=255)
                                                                               gggAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                               >  1:42242/1‑109 (MQ=255)
                                                                               gggAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                               <  2:42242/109‑1 (MQ=255)
                                                                                ggAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCttt                                                                      >  1:508827/1‑149 (MQ=255)
                                                                                 gAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTa                                                                     >  2:350409/1‑149 (MQ=255)
                                                                                  aaaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                    <  1:153243/149‑1 (MQ=255)
                                                                                  aaaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                    <  3:800950/149‑1 (MQ=255)
                                                                                  aaaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                    <  2:320255/149‑1 (MQ=255)
                                                                                   aaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATc                                                                   <  4:179892/149‑1 (MQ=255)
                                                                                       gATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAtt                                                                                          >  1:418089/1‑122 (MQ=255)
                                                                                       gATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAtt                                                                                          <  2:418089/122‑1 (MQ=255)
                                                                                       gATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACt                                                               >  3:262239/1‑149 (MQ=255)
                                                                                        aTAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGc                                                                                                              <  4:804165/101‑1 (MQ=255)
                                                                                        aTAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGc                                                                                                              >  3:804165/1‑101 (MQ=255)
                                                                                          aaTGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCa                                                            >  3:209345/1‑149 (MQ=255)
                                                                                            tGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCAcc                                                          >  3:720436/1‑149 (MQ=255)
                                                                                              tttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGt                                                                                               <  4:746498/110‑1 (MQ=255)
                                                                                              tttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGt                                                                                               >  3:746498/1‑110 (MQ=255)
                                                                                              tttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCt                                                        >  1:164760/1‑149 (MQ=255)
                                                                                              tttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCt                                                        >  4:652098/1‑149 (MQ=255)
                                                                                                tACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGa                                                      >  3:428296/1‑149 (MQ=255)
                                                                                                       aGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCaa                                               <  3:413197/149‑1 (MQ=255)
                                                                                                        gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                              <  2:687591/149‑1 (MQ=255)
                                                                                                        gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                              <  1:688531/149‑1 (MQ=255)
                                                                                                        gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑ACCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                              <  2:489746/149‑1 (MQ=255)
                                                                                                           tCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGt                                           >  4:218143/1‑149 (MQ=255)
                                                                                                            cAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTa                                          >  3:658071/1‑149 (MQ=255)
                                                                                                                   gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                   >  3:152572/1‑149 (MQ=255)
                                                                                                                      aTTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcg                                >  3:450860/1‑149 (MQ=255)
                                                                                                                      aTTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcg                                >  1:523173/1‑149 (MQ=255)
                                                                                                                            ccTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccg                          >  3:408557/1‑149 (MQ=255)
                                                                                                                             cTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgt                         >  3:271125/1‑149 (MQ=255)
                                                                                                                             cTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgt                         >  4:360829/1‑149 (MQ=255)
                                                                                                                              tCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGaa                                                     <  3:768156/120‑1 (MQ=255)
                                                                                                                              tCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGaa                                                     >  4:768156/1‑120 (MQ=255)
                                                                                                                              tCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtgacggtg                        >  1:477886/1‑148 (MQ=255)
                                                                                                                              tCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgtc                        >  2:781017/1‑149 (MQ=255)
                                                                                                                                     gAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGt                 >  3:7170/1‑149 (MQ=255)
                                                                                                                                      aaaTTATTCCCCGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTGACGGTGACGCCgaa                <  2:477886/149‑3 (MQ=255)
                                                                                                                                       aaTTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTc               >  3:311393/1‑149 (MQ=255)
                                                                                                                                           attCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGc           <  2:558171/149‑1 (MQ=255)
                                                                                                                                            ctCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCg          >  2:744007/2‑149 (MQ=255)
                                                                                                                                            ctCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCg          >  2:242602/2‑149 (MQ=255)
                                                                                                                                             tCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTc               >  3:50036/1‑143 (MQ=255)
                                                                                                                                             tCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTc               <  4:50036/143‑1 (MQ=255)
                                                                                                                                             tCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTc               <  4:19292/143‑1 (MQ=255)
                                                                                                                                             tCCCGGG‑ACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTc               >  3:19292/1‑143 (MQ=255)
                                                                                                                                               ccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt       <  4:311393/149‑1 (MQ=255)
                                                                                                                                               ccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt       <  4:658071/149‑1 (MQ=255)
                                                                                                                                               ccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt       <  4:720436/149‑1 (MQ=255)
                                                                                                                                               ccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt       >  3:625954/1‑149 (MQ=255)
                                                                                                                                               ccGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt       >  2:392506/1‑149 (MQ=255)
                                                                                                                                                    aaCGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTGGt  <  4:766966/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCG‑TCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTGGT  >  NC_000913/3169112‑3169408

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