Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,563,341 T→A 5.8% K798I (AAA→ATA)  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,3410TA 5.8% 216.9 / 3.6 69K798I (AAA→ATA) dosPoxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator
Reads supporting (aligned to +/- strand):  ref base T (35/30);  new base A (2/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

TTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGG  >  NC_000913/1563197‑1563485
                                                                                                                                                |                                                                                                                                                
ttGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAg                                                                                                                                              >  1:187544/1‑149 (MQ=255)
  gtgCAGCGTAGTGCAGTTTTGGTGCAAGAGTGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCg                                                                                                                                            <  1:152796/148‑1 (MQ=255)
  gtgCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGc                                                                                                                                             <  2:187544/148‑1 (MQ=255)
      aGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTaa                                                                                                                                        <  4:326169/149‑1 (MQ=255)
                  ttttGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCAAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCAt                                                                                                                            >  4:164692/1‑149 (MQ=255)
                           aaGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCAGAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCTGCCTgg                                                                                                                   <  6:268168/149‑1 (MQ=255)
                           aaGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTaa                                                                                                                                        >  8:256499/1‑128 (MQ=255)
                           aaGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTaa                                                                                                                                        <  7:256499/128‑1 (MQ=255)
                                gggAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCTTGGAAttt                                                                                                              >  1:220134/1‑149 (MQ=255)
                                  gAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGGGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATttct                                                                                                            >  4:121973/1‑149 (MQ=255)
                                    aGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATttcttc                                                                                                          >  2:2980/1‑149 (MQ=255)
                                     gTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcg                                                                                                         <  4:193450/149‑1 (MQ=255)
                                         aggaaggaATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCTTGGAATTTCTTcggcgg                                                                                                     <  2:220134/149‑1 (MQ=255)
                                         aggaaggaATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTcggcgg                                                                                                     <  8:340942/149‑1 (MQ=255)
                                          ggaaggaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGTAATTTGTCAGATTTTCAGCGGTAAAACGCTGCTCATCCAGACTGGAATTTCTTCTGCggg                                                                                                    >  5:274621/1‑149 (MQ=255)
                                           gaaggaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCggg                                                                                                    >  7:197685/1‑148 (MQ=255)
                                              ggaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAgg                                                                                                <  2:443/149‑1 (MQ=255)
                                               gaaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAggg                                                                                               <  2:64677/149‑1 (MQ=255)
                                                aaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAgggg                                                                                              <  4:319877/149‑1 (MQ=255)
                                                aaTCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAgggg                                                                                              <  2:198798/149‑1 (MQ=255)
                                                      ccGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCggg                                                                                         >  5:44754/1‑148 (MQ=255)
                                                      ccGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCggg                                                                                         >  5:174630/1‑148 (MQ=255)
                                                              gTATCTTAACTCACAGTATTTAAGTGACAGTGGCACGTTAAATGAAAACCCGCGAGTGCTGGCGAGAGTAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGaaa                                                                                 <  6:274621/148‑1 (MQ=255)
                                                               tAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATttcttc                                                                                                          >  5:61570/1‑122 (MQ=255)
                                                               tAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATttcttc                                                                                                          <  6:61570/122‑1 (MQ=255)
                                                               tAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAata                                                                               >  4:110652/1‑149 (MQ=255)
                                                                  cTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTATCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATcc                                                                            >  4:228660/1‑149 (MQ=255)
                                                                  cTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTATCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATcc                                                                            >  4:224288/1‑149 (MQ=255)
                                                                     aaaTCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCt                                                                          >  1:306788/1‑148 (MQ=255)
                                                                     aaaTCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCt                                                                          <  8:126108/148‑1 (MQ=255)
                                                                      aaTCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTg                                                                         <  1:2980/148‑1 (MQ=255)
                                                                         cacaGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAAt                                                                      >  8:227046/1‑148 (MQ=255)
                                                                          acaGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATaa                                                                    >  8:282592/1‑149 (MQ=255)
                                                                           caGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCCGCCTGGAATTTCTTCGGCGGGGAGGGGGCGGGAAAAGAAATATCCATGAATAAc                                                                   >  3:62104/1‑149 (MQ=255)
                                                                                 ttAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGTGTAGGGGGCGGGAAAAGAAATATCCCTGACTAACGCGa                                                               >  7:169021/1‑147 (MQ=255)
                                                                                  tAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAg                                                            >  5:195032/1‑149 (MQ=255)
                                                                                   aaGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGt                                                           >  8:118298/1‑149 (MQ=255)
                                                                                     gTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTg                                                          >  2:133607/1‑148 (MQ=255)
                                                                                      tGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATACCGCGACAGTGGa                                                        >  3:215770/1‑149 (MQ=255)
                                                                                      tGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGa                                                        >  1:31210/1‑149 (MQ=255)
                                                                                      tGACAGTATCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGa                                                        >  2:211802/1‑149 (MQ=255)
                                                                                        aCAGTGGCACGTTAAATGAAAACCCGCGAGGGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGAt                                                       <  4:208562/148‑1 (MQ=255)
                                                                                            tgtCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGc                                                  <  7:118298/149‑1 (MQ=255)
                                                                                              tCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCg                                                 >  4:246610/1‑148 (MQ=255)
                                                                                              tCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCg                                                 >  5:45928/1‑148 (MQ=255)
                                                                                              tCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTATCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCg                                                 <  3:228660/148‑1 (MQ=255)
                                                                                              tCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTATCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCg                                                 <  3:224288/148‑1 (MQ=255)
                                                                                                aCGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAg                                              <  2:319180/149‑1 (MQ=255)
                                                                                                  gTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCa                                            >  2:250845/1‑149 (MQ=255)
                                                                                                      aaTGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAAt                                                                                                                <  5:34990/77‑1 (MQ=255)
                                                                                                      aaTGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAAt                                                                                                                >  6:34990/1‑77 (MQ=255)
                                                                                                       aTGAAAACCCGCGAGTGAGGGCCAGAGGAATTTGTCAGATTTTCAGCGGTAACCCGCTGCTCAACCAGCCTGGAATTTCAGCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCAtctc                                        >  2:238167/1‑148 (MQ=255)
                                                                                                           aaaCCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAAtt                                   <  2:306788/149‑1 (MQ=255)
                                                                                                            aaCCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTg                                  >  2:215615/1‑149 (MQ=255)
                                                                                                             aCCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGc                                 <  7:282592/149‑1 (MQ=255)
                                                                                                               ccGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCTGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGctc                               <  8:259097/149‑1 (MQ=255)
                                                                                                               ccGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGctc                               <  7:295249/149‑1 (MQ=255)
                                                                                                                   gAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTg                           <  6:44754/149‑1 (MQ=255)
                                                                                                                   gAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTg                           <  6:174630/149‑1 (MQ=255)
                                                                                                                   gAGTGCGGGCAAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTg                           <  3:164692/149‑1 (MQ=255)
                                                                                                                      tGCGCGCGACAGGAATTTCTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGtt                        >  2:312583/1‑149 (MQ=255)
                                                                                                                           gCGAGAGGAATTTGCCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCtt                             >  5:195050/1‑139 (MQ=255)
                                                                                                                           gCGAGAGGAATTTGCCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCtt                             <  6:195050/139‑1 (MQ=255)
                                                                                                                             gagaGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGc                 <  3:110652/149‑1 (MQ=255)
                                                                                                                                  gAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTcc            >  4:91140/1‑149 (MQ=255)
                                                                                                                                   aaTTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCg           <  7:227046/149‑1 (MQ=255)
                                                                                                                                         tCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCacga     >  8:336491/1‑149 (MQ=255)
                                                                                                                                           aGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCacgacg   >  1:305734/1‑149 (MQ=255)
                                                                                                                                            gATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACgg  >  1:115115/1‑149 (MQ=255)
                                                                                                                                              ttttCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCa                                       >  3:252934/1‑110 (MQ=255)
                                                                                                                                              ttttCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCa                                       <  4:252934/110‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TTGTGCAGCGTAGTGCAGTTTTGGTGCAAGAGGGGAAGTTAAGGAAGGAATCTCCCGGAATCGTAGCTGAAATCACAGTATTTAAGTGACAGTGTCACGTTAAATGAAAACCCGCGAGTGCGGGCGAGAGGAATTTGTCAGATTTTCAGCGGTAACACGCTGCTCATCCAGCCTGGAATTTCTTCGGCGGGTAGGGGGCGGGAAAAGAAATATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGG  >  NC_000913/1563197‑1563485

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: