Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 874,628:1 +G 5.5% coding (1650/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,6281.G 5.5% 494.2 / 22.4 109coding (1650/2349 nt)yliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base . (52/51);  new base G (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 = 6.52e-01

Only 100 of 112 total aligned reads displayed.

CGATTTCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAA  >  NC_000913/874480‑874776
                                                                                                                                                     |                                                                                                                                                    
cGATTTCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATc                                                                                                                                                       <  2:1114790/149‑1 (MQ=255)
   tttCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGTGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑Gtt                                                                                                                                                   <  1:1223219/149‑1 (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                                                                                                                                              <  4:1492670/142‑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)
                            ctcAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTgg                                                                                                                          <  2:945190/149‑1 (MQ=255)
                             tcAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTggg                                                                                                                         <  1:1081618/149‑1 (MQ=255)
                             tcAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTggg                                                                                                                         >  2:1203823/1‑149 (MQ=255)
                             tcAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTggg                                                                                                                         >  4:787269/1‑149 (MQ=255)
                               agagTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGgcg                                                                                                                       >  2:1354448/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‑GTTGATATTaa                                                                                                                                           >  2:381801/1‑107 (MQ=255)
                                                     gcgAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCgcgc                                                                                                                                  <  2:1542799/116‑1 (MQ=255)
                                                     gcgAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCgcgc                                                                                                                                  >  1:1542799/1‑116 (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)
                                                                       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                                                                         >  3:1088863/1‑149 (MQ=255)
                                                                             aaTATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTc                                                                         >  4:247049/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)
                                                                                  cTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGc                                                                    <  3:247049/149‑1 (MQ=255)
                                                                                      aTTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGaa                                                                >  1:538494/1‑149 (MQ=255)
                                                                                       ttGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAg                                                               >  4:844467/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                                                                                                                                                  >  2:1744981/1‑58 (MQ=255)
                                                                                               aaTATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTg                                                                                                                                                  <  1:1744981/58‑1 (MQ=255)
                                                                                               aaTATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTaa                                                                                                                                           <  3:1397981/65‑1 (MQ=255)
                                                                                               aaTATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTaa                                                                                                                                           >  4:1397981/1‑65 (MQ=255)
                                                                                                  atTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgggg                                                                                  <  3:1524483/119‑1 (MQ=255)
                                                                                                  atTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgggg                                                                                  >  4:1524483/1‑119 (MQ=255)
                                                                                                  atTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgg                                                                                    >  1:1830911/1‑117 (MQ=255)
                                                                                                  atTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCgg                                                                                    <  2:1830911/117‑1 (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)
                                                                                                                          gAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGatat                                        >  1:318154/1‑137 (MQ=255)
                                                                                                                           aaCTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAAt                           <  3:844467/149‑1 (MQ=255)
                                                                                                                              tATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCggg                                                                                   <  1:1457546/90‑1 (MQ=255)
                                                                                                                              tATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCggg                                                                                   >  2:1457546/1‑90 (MQ=255)
                                                                                                                              tATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCAcc                                                                      >  3:1450108/1‑103 (MQ=255)
                                                                                                                                  ccTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTgg                    <  1:1382604/149‑1 (MQ=255)
                                                                                                                                   cTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGa                         >  4:799898/1‑143 (MQ=255)
                                                                                                                                   cTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGa                         <  3:799898/143‑1 (MQ=255)
                                                                                                                                       tCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCgg                                     >  1:581847/1‑127 (MQ=255)
                                                                                                                                        cTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTcg                                                                        >  3:1790534/1‑90 (MQ=255)
                                                                                                                                        cTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTATGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCa                                                                        <  4:1790534/91‑1 (MQ=255)
                                                                                                                                        atATCAA‑CGATCGGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGaa                                             >  3:258904/2‑118 (MQ=255)
                                                                                                                                         tATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCa             <  2:177670/149‑1 (MQ=255)
                                                                                                                                         tATCAA‑CGATCGGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGa                                              >  1:188547/1‑116 (MQ=255)
                                                                                                                                         tATCAA‑CGATCGGTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGa                                              <  2:188547/116‑1 (MQ=255)
                                                                                                                                              aCCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGat                                          <  1:584746/115‑1 (MQ=255)
                                                                                                                                              aCCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGat                                          >  2:584746/1‑115 (MQ=255)
                                                                                                                                               ccGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCaa                            <  3:971650/128‑1 (MQ=255)
                                                                                                                                                    c‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTaaaa  <  4:1117287/149‑1 (MQ=255)
                                                                                                                                                    c‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTaaaa  <  4:724663/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGATTTCCAGTATCCGTGACCATATGTTCTCAGAGTTTATTTTCTGTGATGACGCGAAACTCAACGAAATAGAAGAGAATATCTGGATTGCGCGTAATATTCGCCATGCAATGGAAATTGGCGAACTATTCCTCGTCTATCAACCGATC‑GTTGATATTAACACCCGCGCCATTCTGGGCGCGGAGGCGTTGTGCCGTTGGGTGTCTGCGGAGCGGGGGATCATTTCACCGCTGAAGTTCATTACCATTGCTGAAGATATCGGGTTTATCAATGAGCTGGGTTATCAGATTATTAAAA  >  NC_000913/874480‑874776

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: ‑