Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 874,623 Δ1 bp 5.5% coding (1645/2349 nt) yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913874,6230C. 5.5% 392.9 / 15.9 109coding (1645/2349 nt)yliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (52/51);  new base . (3/3);  total (55/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 113 total aligned reads displayed.

TGGCGATTTCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCA  >  NC_000913/874477‑874765
                                                                                                                                                  |                                                                                                                                               
tgGCGATTTCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCg                                                                                                                                               >  1:1114790/1‑149 (MQ=255)
      tttCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGTGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑Gtt                                                                                                                                        <  1:1223219/149‑1 (MQ=255)
        tCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGa                                                                                                                                      >  1:1634039/1‑149 (MQ=255)
           aGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGatat                                                                                                                                   <  1:970327/149‑1 (MQ=255)
            gTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATAtt                                                                                                                                  >  1:945190/1‑149 (MQ=255)
            gTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAA‑CGATCGGTTGat                                                                                                                                     <  3:450218/146‑1 (MQ=255)
            gTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAA‑CGATCGGTTGat                                                                                                                                     >  4:450218/1‑146 (MQ=255)
               tCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAac                                                                                                                               <  2:395716/149‑1 (MQ=255)
                  gTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGGAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGatat                                                                                                                                   >  3:1492670/1‑142 (MQ=255)
                  gTGACCATACGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACAcc                                                                                                                            >  4:1800409/1‑149 (MQ=255)
                        atatGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGcc                                                                                                                      >  3:1359290/1‑149 (MQ=255)
                         tatGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCa                                                                                                                     >  3:1824733/1‑149 (MQ=255)
                             ttCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCt                                                                                                                 >  4:896155/1‑149 (MQ=255)
                                tcAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTggg                                                                                                              >  2:1203823/1‑149 (MQ=255)
                                tcAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTggg                                                                                                              >  4:787269/1‑149 (MQ=255)
                                tcAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTggg                                                                                                              <  1:1081618/149‑1 (MQ=255)
                                tcAGAGTTTATCTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTCTCAACCGATC‑GTTGGTATTAACACCCGCGCCATTCTggg                                                                                                              >  3:1244752/1‑149 (MQ=255)
                                  agagTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGgcg                                                                                                            >  2:1354448/1‑149 (MQ=255)
                                  agagTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGgcg                                                                                                            >  2:281329/1‑149 (MQ=255)
                                      tttATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGa                                                                                                        >  3:724663/1‑149 (MQ=255)
                                       ttATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAg                                                                                                       <  3:787269/149‑1 (MQ=255)
                                            ttCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGtt                                                                                                  <  3:896155/149‑1 (MQ=255)
                                                gtgATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGc                                                                                              <  1:1203823/149‑1 (MQ=255)
                                                    tgaCGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGtt                                                                                          >  2:62050/1‑149 (MQ=255)
                                                    tgaCGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGtt                                                                                          <  1:1354448/149‑1 (MQ=255)
                                                      aCGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTGATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTgg                                                                                        <  4:1244752/149‑1 (MQ=255)
                                                       cgcgAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTggg                                                                                       <  1:281329/149‑1 (MQ=255)
                                                        gcgAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTaa                                                                                                                                <  1:381801/107‑1 (MQ=255)
                                                        gcgAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCgcgc                                                                                                                       <  2:1542799/116‑1 (MQ=255)
                                                        gcgAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCgcgc                                                                                                                       >  1:1542799/1‑116 (MQ=255)
                                                          gAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAac                                                                                                                               >  3:809836/1‑106 (MQ=255)
                                                          gAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAac                                                                                                                               <  4:809836/106‑1 (MQ=255)
                                                           aaaCTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTc                                                                                   <  2:1634039/149‑1 (MQ=255)
                                                               tCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCg                                                                               >  1:1804288/1‑149 (MQ=255)
                                                                 aaCGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTaa                                                                                                                                >  3:1118470/1‑98 (MQ=255)
                                                                 aaCGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTaa                                                                                                                                <  4:1118470/98‑1 (MQ=255)
                                                                 aaCGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAacac                                                                                                                             >  4:152205/1‑101 (MQ=255)
                                                                 aaCGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAacac                                                                                                                             <  3:152205/101‑1 (MQ=255)
                                                                 aaCGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGa                                                                             >  4:1232986/1‑149 (MQ=255)
                                                                 aaCGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGa                                                                             >  2:351216/1‑149 (MQ=255)
                                                                    gAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCg                                                                          >  1:1111653/1‑149 (MQ=255)
                                                                     aaaTAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgg                                                                         <  1:351216/149‑1 (MQ=255)
                                                                         agaagaGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAAc                                                                                                                                                 >  3:377563/1‑74 (MQ=255)
                                                                          gaagaGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCa                                                                                                                     >  2:343873/1‑100 (MQ=255)
                                                                          gaagaGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCa                                                                                                                     <  1:343873/100‑1 (MQ=255)
                                                                          gaagaGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTgg                                                                                        >  4:646567/1‑129 (MQ=255)
                                                                          gaagaGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTgg                                                                                        <  3:646567/129‑1 (MQ=255)
                                                                             gagaATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCAt                                                                 <  4:1824733/149‑1 (MQ=255)
                                                                               gaATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGata                                                                                                                                    <  3:1519725/80‑1 (MQ=255)
                                                                               gaATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGata                                                                                                                                    >  4:1519725/1‑80 (MQ=255)
                                                                               gaATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCAttt                                                               >  2:1026263/1‑149 (MQ=255)
                                                                                aaTATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTc                                                              >  4:247049/1‑149 (MQ=255)
                                                                                aaTATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTc                                                              >  3:1088863/1‑149 (MQ=255)
                                                                                   atCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGc                                                                                              >  2:983577/1‑114 (MQ=255)
                                                                                   atCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGc                                                                                              <  1:983577/114‑1 (MQ=255)
                                                                                    tCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTg                                                                                         >  2:912682/1‑118 (MQ=255)
                                                                                    tCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTg                                                                                         <  1:912682/118‑1 (MQ=255)
                                                                                     cTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGc                                                         <  3:247049/149‑1 (MQ=255)
                                                                                       ggATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTg                                                       <  3:1800409/149‑1 (MQ=255)
                                                                                         aTTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGaa                                                     >  1:538494/1‑149 (MQ=255)
                                                                                          ttGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAg                                                    >  4:844467/1‑149 (MQ=255)
                                                                                              gcgTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCa                                                >  2:1382604/1‑149 (MQ=255)
                                                                                               cgTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCAt                                               <  2:1804288/149‑1 (MQ=255)
                                                                                                gTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTtgt                                                                                                >  1:1353844/1‑99 (MQ=255)
                                                                                                gTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTtgt                                                                                                <  2:1353844/99‑1 (MQ=255)
                                                                                                  aaTATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTg                                                                                                                                       <  1:1744981/58‑1 (MQ=255)
                                                                                                  aaTATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTg                                                                                                                                       >  2:1744981/1‑58 (MQ=255)
                                                                                                  aaTATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTaa                                                                                                                                >  4:1397981/1‑65 (MQ=255)
                                                                                                  aaTATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTaa                                                                                                                                <  3:1397981/65‑1 (MQ=255)
                                                                                                     atTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgggg                                                                       <  3:1524483/119‑1 (MQ=255)
                                                                                                     atTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgggg                                                                       >  4:1524483/1‑119 (MQ=255)
                                                                                                     atTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgg                                                                         <  2:1830911/117‑1 (MQ=255)
                                                                                                     atTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgg                                                                         >  1:1830911/1‑117 (MQ=255)
                                                                                                      ttCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCAtt                                        >  1:401699/1‑149 (MQ=255)
                                                                                                        cGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCAttac                                      <  2:538494/149‑3 (MQ=255)
                                                                                                           cATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGa                                   >  1:177670/1‑149 (MQ=255)
                                                                                                            aTGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGaa                                  <  4:1088863/149‑1 (MQ=255)
                                                                                                            aTGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGaa                                  <  3:1232986/149‑1 (MQ=255)
                                                                                                               cAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCggg                                                                        >  2:375452/1‑108 (MQ=255)
                                                                                                               cAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCggg                                                                        <  1:375452/108‑1 (MQ=255)
                                                                                                                 aTGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGatat                             <  2:1111653/149‑1 (MQ=255)
                                                                                                                   ggAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCg                           <  1:62050/149‑1 (MQ=255)
                                                                                                                        ttGGCGAACTATTCCTCGTATATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGttt                      >  3:1117287/1‑149 (MQ=255)
                                                                                                                           gCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTa                     >  1:423133/1‑147 (MQ=255)
                                                                                                                           gCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTa                     <  2:423133/147‑1 (MQ=255)
                                                                                                                             gAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGatat                             <  2:318154/137‑1 (MQ=255)
                                                                                                                                 tATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCggg                                                                        >  2:1457546/1‑90 (MQ=255)
                                                                                                                                 tATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCAcc                                                           <  4:1450108/103‑1 (MQ=255)
                                                                                                                                     ccTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTgg         <  1:1382604/149‑1 (MQ=255)
                                                                                                                                      cTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGa              <  3:799898/143‑1 (MQ=255)
                                                                                                                                      cTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGa              >  4:799898/1‑143 (MQ=255)
                                                                                                                                          tCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCgg                          >  1:581847/1‑127 (MQ=255)
                                                                                                                                          tCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCgg                          <  2:581847/127‑1 (MQ=255)
                                                                                                                                           cTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTcg                                                             >  3:1790534/1‑90 (MQ=255)
                                                                                                                                           cTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTATGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCa                                                             <  4:1790534/91‑1 (MQ=255)
                                                                                                                                           atATCAA‑CGATCGGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGaa                                  <  4:258904/117‑1 (MQ=255)
                                                                                                                                            tATCAA‑CGATCGGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGa                                   >  1:188547/1‑116 (MQ=255)
                                                                                                                                            tATCAA‑CGATCGGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGa                                   <  2:188547/116‑1 (MQ=255)
                                                                                                                                                 aCCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGat                               >  2:584746/1‑115 (MQ=255)
                                                                                                                                                  ccGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCaa                 <  3:971650/128‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TGGCGATTTCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCA  >  NC_000913/874477‑874765

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