Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,563,432 T→A 5.3% I768F (ATC→TTC)  dosP ← oxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,563,4320TA 5.3% 226.8 / 3.6 75I768F (ATC→TTC) dosPoxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator
Reads supporting (aligned to +/- strand):  ref base T (38/32);  new base A (2/2);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAA  >  NC_000913/1563284‑1563578
                                                                                                                                                    |                                                                                                                                                  
gACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGAt                                                                                                                                                    >  3:98696/1‑149 (MQ=255)
             tAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCAtc                                                                                                                                       >  6:108091/1‑149 (MQ=255)
              aaaTGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCAtc                                                                                                                                       <  2:167695/148‑1 (MQ=255)
               aaTGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCAtct                                                                                                                                      >  1:172442/1‑148 (MQ=255)
               aaTGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGAACTTGCGTAGCAtct                                                                                                                                      >  8:291640/1‑148 (MQ=255)
               aaTGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGAACTTGCGTAGCAtct                                                                                                                                      >  8:286504/1‑148 (MQ=255)
                aTGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCAtctc                                                                                                                                     >  2:101977/1‑148 (MQ=255)
                 tGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGc                                                                                                                                          >  5:188958/1‑142 (MQ=255)
                 tGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGc                                                                                                                                          <  6:188958/142‑1 (MQ=255)
                 tGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCa                                                                                                                                    >  3:259817/1‑148 (MQ=255)
                 tGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCa                                                                                                                                    >  3:270165/1‑148 (MQ=255)
                     aaCCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTg                                                                                                                               >  5:180802/1‑149 (MQ=255)
                          gcgcGTGCGGGCGAGAGGAATGTGTCAGATCTGCAGCGGTACCCCGCGGCTCACCCCGCCGCGAATTGCTTCGGCGGGTAGGGGGCGGGAACAGAAATATCCCTGCATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCtt                                                                                                                          <  2:239368/149‑1 (MQ=255)
                              gTGCGGGCGAGATGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGTAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGt                                                                                                                      <  3:85157/149‑1 (MQ=255)
                                 cGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGttt                                                                                                                    <  3:101671/148‑1 (MQ=255)
                                      gagaGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGc                                                                                                              >  4:132836/1‑149 (MQ=255)
                                           gAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCttt                                                                                                                         >  3:90041/1‑133 (MQ=255)
                                           gAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCttt                                                                                                                         <  4:90041/133‑1 (MQ=255)
                                           gAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTcc                                                                                                         >  5:41248/1‑149 (MQ=255)
                                             aTTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGAGAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGc                                                                                                       >  6:264192/1‑149 (MQ=255)
                                              tttGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGcc                                                                                                      <  1:101977/149‑1 (MQ=255)
                                                  tCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCacga                                                                                                  <  3:229452/149‑1 (MQ=255)
                                                    aGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCg                                                                                                                                              <  2:83198/103‑1 (MQ=255)
                                                    aGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCg                                                                                                                                              >  1:83198/1‑103 (MQ=255)
                                                     gATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGCCAGTGGATCTTGCGTAGCATCTCAGATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACgg                                                                                               >  1:97377/1‑149 (MQ=255)
                                                              cGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGAACTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGa                                                                                      <  7:286504/149‑1 (MQ=255)
                                                              cGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGAACTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGa                                                                                      <  7:291640/149‑1 (MQ=255)
                                                              cGGTAACACGCTGCTCATCCAGACTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAAtt                                                                                        >  1:247638/1‑147 (MQ=255)
                                                               ggTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAg                                                                                     <  3:9011/149‑1 (MQ=255)
                                                               ggTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGAGCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAg                                                                                     <  4:98696/149‑1 (MQ=255)
                                                                gTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAgg                                                                                    <  2:172442/149‑1 (MQ=255)
                                                                     acGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTg                                                                               >  7:234601/1‑149 (MQ=255)
                                                                     acGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTg                                                                               >  8:248330/1‑149 (MQ=255)
                                                                     acGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTg                                                                               >  8:325634/1‑149 (MQ=255)
                                                                            tCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATAGCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCAAGACGGTTAAATTGAGGCTTTGACCAATg                                                                        >  6:180820/1‑149 (MQ=255)
                                                                              aTCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCt                                                                      >  6:245661/1‑149 (MQ=255)
                                                                                 cAGCCTTGATTTTGTTCGGCGTGTATGGTGCGGGAAAAGAAATATCCCTGAATTACTCGACAGGTGATCTTGCGTAGCATCTCAATTTGCTCTTTGTTTTCGACTCCTTCCGCCACGACTGTTAAATTGATGCTTTGCCCTATGCTgg                                                                    <  6:201783/148‑1 (MQ=255)
                                                                                 cAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTgg                                                                    <  6:41248/148‑1 (MQ=255)
                                                                                 cAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTgg                                                                    <  5:108091/148‑1 (MQ=255)
                                                                                 cAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGt                                                                   <  8:234601/149‑1 (MQ=255)
                                                                                    ccTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAAt                                                                <  7:325634/149‑1 (MQ=255)
                                                                                      tGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAg                                                                                     <  4:147828/126‑1 (MQ=255)
                                                                                      tGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAg                                                                                     >  3:147828/1‑126 (MQ=255)
                                                                                      tGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAAt                                                                >  8:299044/1‑147 (MQ=255)
                                                                                      tGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAAt                                                                <  7:299044/147‑1 (MQ=255)
                                                                                      tGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATgg                                                              >  4:109870/1‑149 (MQ=255)
                                                                                       ggAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCACTGAATAACGCGACTGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGCCGCCTTCCGCCCCGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGc                                                             >  8:129206/1‑149 (MQ=255)
                                                                                        gAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGc                                                             <  6:297713/148‑1 (MQ=255)
                                                                                         aaTTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCtt                                                           >  1:356076/1‑149 (MQ=255)
                                                                                         aaTTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCtt                                                           <  7:248330/149‑1 (MQ=255)
                                                                                          aTTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCtt                                                           >  5:126195/1‑148 (MQ=255)
                                                                                                 cggcggTGAGGGGGCTGGAAAATAAATATCCCTGAATAACGCGACAGTGTATCTTGCGTATAATCTCAAATTTCTCTTTGTTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGt                                                      <  5:118174/146‑1 (MQ=255)
                                                                                                 cggcggGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTaa                                                    <  4:270165/148‑1 (MQ=255)
                                                                                                  ggcggGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAg                                                   <  2:247638/148‑1 (MQ=255)
                                                                                                  gccggGTAGGGGGCGGGAAAAGCAATATCCCTCAATCACGCGACAGTGGATCTTGCGTAGCATCTCAAATTCCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAg                                                   <  4:9621/146‑1 (MQ=255)
                                                                                                     ggggAGGGGGCGGGAAAAGAAATTTCCCTGAATACCGCGAACGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCGGCCACGACGGTTAAATTGAGTCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCaa                                               <  4:259817/149‑1 (MQ=255)
                                                                                                               cGGGAAAAGAAATATCCCTGAATACCGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGttt                                     >  5:227984/1‑149 (MQ=255)
                                                                                                               cGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGttt                                     >  1:152862/1‑149 (MQ=255)
                                                                                                                gggAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGtttt                                    >  8:65523/1‑149 (MQ=255)
                                                                                                                    aaaGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCgg                                >  8:220658/1‑149 (MQ=255)
                                                                                                                      aGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCtt                                                           <  7:356317/120‑1 (MQ=255)
                                                                                                                      aGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCtt                                                           >  8:356317/1‑120 (MQ=255)
                                                                                                                      aGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTc                              <  3:132836/149‑1 (MQ=255)
                                                                                                                       gAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCa                             <  8:73798/149‑1 (MQ=255)
                                                                                                                             tCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAAc                       >  6:118083/1‑149 (MQ=255)
                                                                                                                               ccTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGa                     >  8:21048/1‑149 (MQ=255)
                                                                                                                                cTGTATAACGCGACGGTGGATCTTTCGTAGCATCTCAAAGTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGCGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGAt                    <  6:280289/149‑1 (MQ=255)
                                                                                                                                 tGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATc                   >  7:29107/1‑149 (MQ=255)
                                                                                                                                   aaTAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGa                 <  7:220658/149‑1 (MQ=255)
                                                                                                                                   aaTAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGa                 <  7:21048/149‑1 (MQ=255)
                                                                                                                                               gTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGATGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGt     >  8:273481/1‑149 (MQ=255)
                                                                                                                                                tGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGCCGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGt     >  2:92479/1‑148 (MQ=255)
                                                                                                                                                tGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGt     >  3:45592/1‑148 (MQ=255)
                                                                                                                                                  gATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCaa  >  4:47221/1‑149 (MQ=255)
                                                                                                                                                   aTCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGATGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtttt       >  3:302212/1‑143 (MQ=255)
                                                                                                                                                   aTCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGATGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtttt       <  4:302212/143‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAA  >  NC_000913/1563284‑1563578

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: