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

CCACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTGGT  >  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)
           aTTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCgg                                                                                                                                           <  3:444065/149‑1 (MQ=255)
           aTTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCgg                                                                                                                                           >  2:452536/1‑149 (MQ=255)
            tttCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCggg                                                                                                                                          >  2:83233/1‑149 (MQ=255)
                gCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCg                                                                                                                                      >  4:413197/1‑149 (MQ=255)
                 ccACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCgg                                                                                                                                     >  2:148972/1‑149 (MQ=255)
                       tGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACaa                                                                                                                               >  2:153243/1‑149 (MQ=255)
                       tGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACaa                                                                                                                               <  1:654101/149‑1 (MQ=255)
                        gTAGTCGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAAt                                                                                                                              >  3:552792/1‑149 (MQ=255)
                              gAGCAATTCAGGGCTACAGCGGTGCGTCAACATCGCCAGAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGGTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGa                                                                                                                        <  2:392077/149‑1 (MQ=255)
                              gAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGa                                                                                                                        <  1:652230/149‑1 (MQ=255)
                              gAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGa                                                                                                                        <  2:610621/149‑1 (MQ=255)
                                    ttCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACt                                                                                                                  >  1:489746/1‑149 (MQ=255)
                                    ttCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACt                                                                                                                  >  1:687591/1‑149 (MQ=255)
                                       aGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                               >  3:766966/1‑149 (MQ=255)
                                         ggCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCg                                                                                                             <  4:19938/149‑1 (MQ=255)
                                            tACAGCGGTGCGGCAACATCGCCGCAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGAt                                                                                                          <  2:271560/149‑1 (MQ=255)
                                                gCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGt                                                                                                      <  1:436933/149‑1 (MQ=255)
                                                gCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGt                                                                                                      <  3:669416/149‑1 (MQ=255)
                                                  ggTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                    <  4:798671/149‑1 (MQ=255)
                                                  ggTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGttt                                                                                                    <  1:238798/149‑1 (MQ=255)
                                                    tGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCTGGAGCGGTAACAATTACGGATTCACTTTGCGGATCTGTTTGa                                                                                                  <  3:472027/149‑1 (MQ=255)
                                                       ggCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACTCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGt                                                                                               <  2:688535/149‑1 (MQ=255)
                                                       ggCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGt                                                                                               <  2:68734/149‑1 (MQ=255)
                                                             aTCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAttt                                                                                         >  1:558171/1‑149 (MQ=255)
                                                              tCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTc                                                                                        <  3:190053/149‑1 (MQ=255)
                                                                gCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCaa                                                                                      >  4:800950/1‑149 (MQ=255)
                                                                gCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCaa                                                                                      <  3:224570/149‑1 (MQ=255)
                                                                 ccACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAAc                                                                                     <  1:148972/149‑1 (MQ=255)
                                                                   acaGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAg                                                                                   >  3:639040/1‑149 (MQ=255)
                                                                   acaGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAg                                                                                   >  1:508673/1‑149 (MQ=255)
                                                                      ggTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAAt                                                                                >  2:504734/1‑149 (MQ=255)
                                                                          ttGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCa                                                                            >  2:116693/1‑149 (MQ=255)
                                                                            gATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGt                                                                          <  1:83233/149‑1 (MQ=255)
                                                                               gggAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                               >  1:42242/1‑109 (MQ=255)
                                                                               gggAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTg                                                                                                               <  2:42242/109‑1 (MQ=255)
                                                                                ggAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCttt                                                                      >  1:508827/1‑149 (MQ=255)
                                                                                 gAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTa                                                                     >  2:350409/1‑149 (MQ=255)
                                                                                  aaaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                    <  1:153243/149‑1 (MQ=255)
                                                                                  aaaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                    <  3:800950/149‑1 (MQ=255)
                                                                                  aaaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTAt                                                                    <  2:320255/149‑1 (MQ=255)
                                                                                   aaTCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATc                                                                   <  4:179892/149‑1 (MQ=255)
                                                                                       gATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAtt                                                                                          >  1:418089/1‑122 (MQ=255)
                                                                                       gATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAAtt                                                                                          <  2:418089/122‑1 (MQ=255)
                                                                                       gATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACt                                                               >  3:262239/1‑149 (MQ=255)
                                                                                        aTAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGc                                                                                                              <  4:804165/101‑1 (MQ=255)
                                                                                        aTAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGc                                                                                                              >  3:804165/1‑101 (MQ=255)
                                                                                          aaTGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCa                                                            >  3:209345/1‑149 (MQ=255)
                                                                                            tGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCAcc                                                          >  3:720436/1‑149 (MQ=255)
                                                                                              tttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGt                                                                                               <  4:746498/110‑1 (MQ=255)
                                                                                              tttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGt                                                                                               >  3:746498/1‑110 (MQ=255)
                                                                                              tttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCt                                                        >  1:164760/1‑149 (MQ=255)
                                                                                              tttACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCt                                                        >  4:652098/1‑149 (MQ=255)
                                                                                                tACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGa                                                      >  3:428296/1‑149 (MQ=255)
                                                                                                       aGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCaa                                               <  3:413197/149‑1 (MQ=255)
                                                                                                        gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                              <  2:687591/149‑1 (MQ=255)
                                                                                                        gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                              <  1:688531/149‑1 (MQ=255)
                                                                                                        gtgtCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGACCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAAc                                              <  2:489746/149‑1 (MQ=255)
                                                                                                           tCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGt                                           >  4:218143/1‑149 (MQ=255)
                                                                                                            cAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTa                                          >  3:658071/1‑149 (MQ=255)
                                                                                                                   gTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCg                                   >  3:152572/1‑149 (MQ=255)
                                                                                                                      aTTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcg                                >  3:450860/1‑149 (MQ=255)
                                                                                                                      aTTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcg                                >  1:523173/1‑149 (MQ=255)
                                                                                                                            ccTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccg                          >  3:408557/1‑149 (MQ=255)
                                                                                                                             cTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgt                         >  3:271125/1‑149 (MQ=255)
                                                                                                                             cTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgt                         >  4:360829/1‑149 (MQ=255)
                                                                                                                              tCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGaa                                                     <  3:768156/120‑1 (MQ=255)
                                                                                                                              tCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGaa                                                     >  4:768156/1‑120 (MQ=255)
                                                                                                                              tCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtgacggtg                        >  1:477886/1‑148 (MQ=255)
                                                                                                                              tCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGcgtcaccgtc                        >  2:781017/1‑149 (MQ=255)
                                                                                                                                     gAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGt                 >  3:7170/1‑149 (MQ=255)
                                                                                                                                      aaaTTATTCCCCGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTGACGGTGACGCCgaa                <  2:477886/149‑3 (MQ=255)
                                                                                                                                       aaTTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTc               >  3:311393/1‑149 (MQ=255)
                                                                                                                                           attCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGc           <  2:558171/149‑1 (MQ=255)
                                                                                                                                            ctCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCg          >  2:744007/2‑149 (MQ=255)
                                                                                                                                            ctCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCg          >  2:242602/2‑149 (MQ=255)
                                                                                                                                             tCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTc               >  3:50036/1‑143 (MQ=255)
                                                                                                                                             tCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTc               <  4:50036/143‑1 (MQ=255)
                                                                                                                                             tCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTc               <  4:19292/143‑1 (MQ=255)
                                                                                                                                             tCCCGGGACGCGTTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTc               >  3:19292/1‑143 (MQ=255)
                                                                                                                                               ccGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt       <  4:311393/149‑1 (MQ=255)
                                                                                                                                               ccGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt       <  4:658071/149‑1 (MQ=255)
                                                                                                                                               ccGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt       <  4:720436/149‑1 (MQ=255)
                                                                                                                                               ccGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt       >  3:625954/1‑149 (MQ=255)
                                                                                                                                               ccGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGttt       >  2:392506/1‑149 (MQ=255)
                                                                                                                                                    aaCGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTGGt  <  4:766966/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCACATGATAAATTTCGCCACGCTGTAGTTGAGCAATTCAGGGCTACAGCGGTGCGGCAACATCGCCACAGGTATTGATGGGAAATCGATAATGTTTACATCCAGTGTCAGGTTTGTCATTTTGCCTCCCTGCGAAATTATTCCCGGGAACGCGTCCCGGGAGCGGTAACAATTACGGATTTACTTTGCGGATCTGTTTGACGTCAATTTCAACAGAATTCCAGTCTTTATCGACTTCACCCTGAATCTCAACCGTATCTTTCGGCGTCACCGTCACGCCGTTCCAGCGTTTGTGGT  >  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: